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And so long as we have this itch 
for explanations, must we not 
always carry around with us this 
cumbersome but precious bag of 
clues called History? 
— G. Swift, Waterland (1983) 
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A Primer for Black Hole Quantum Physics 



Abstract: The mechanisms which give rise to Hawking radiation are 

revealed by analyzing in detail pair production in the presence of horizons. In 
preparation for the black hole problem, three preparatory problems are dwelt 
with at length: pair production in an external electric field, thermalization 
of a uniformly accelerated detector and accelerated mirrors. In the light of 
these examples, the black hole evaporation problem is then presented. 

The leitmotif is the singular behavior of modes on the horizon which 
gives rise to a steady rate of production. Special emphasis is put on how 
each produced particle contributes to the mean albeit arising from a partic- 
ular vacuum fluctuation. It is the mean which drives the semiclassical back 
reaction. This aspect is analyzed in more detail than heretofore and in par- 
ticular its drawbacks are emphasized. It is the semiclassical theory which 
gives rise to Hawking's famous equation for the loss of mass of the black 
hole due to evaporation dM/dt ~ — Black hole thermodynamics is 

derived from the evaporation process whereupon the reservoir character of 
the black hole is manifest. The relation to the thermodynamics of the eternal 
black hole through the Hartle-Hawking vacuum and the Killing identity are 
displayed. 

It is through the analysis of the fluctuations of the held conflgurations 
which give rise to a particular Hawking photon that the dubious character of 
the semiclassical theory is manifest. The present frontier of research revolves 
around this problem and is principally concerned with the fact that one calls 
upon energy scales that are greater than Planckian and the possibility of a 
non unitary evolution as well. These last subjects are presented in qualitative 
fashion only, so that this review stops at the threshold of quantum gravity. 
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Introduction 



This review is conceived as a pedagogical essay on black hole quantum 
physics. 

Black hole physics has a rich and varied history. Indeed a complete cov- 
erage of the subject would run parallel to a good part of the development 
of modern theoretical physics, embracing as it does general relativity, quan- 
tum field theory and thermodynamics. The present review, a primer that 
is designed to bring the reader to the threshold of current research, covers 
the middle ground. It is a review of that aspect of the theory that dates 
from Hawking's seminal work on black hole evaporation up to but not in- 
cluding quantum gravity. This latter chapter has been the subject of intense 
investigation over the past few years and is still at a very speculative stage, 
certainly not yet in a state to receive consecration in a review article. What 
has emerged from investigation of present quantum field theory, in the pres- 
ence of a black hole is the awareness that the problem cannot be confronted 
without at least some aspects that emerge from the (as yet unknown) quan- 
tum theory of gravity. It is our intention here to develop the theory to the 
point where the reader will have a clear idea of the nature of the unsolved 
problems given the present state of our ignorance. An optimist would say 
that the solution to these problems leads to the discovery of the quantum 
theory of gravity. But even a minimalist will admit that a good deal is to 
be learnt from black hole physics at its present stage of development. It is 
perhaps not too much to hope that a careful presentation in terms of present 
concepts can illuminate the way to the unknown. Such is our endeavor. 

There exist several situations in quantum field theory which give rise to 
phenomena which are similar to black hole evaporation. Our own experience 



7 



has been that it is very helpful to analyze those examples in detail. They 
furnish a body of information so as to constitute a useful theoretical labora- 
tory for the study of the more complex issues which arise in the black hole 
problem. To this end the initial chapters of this review are devoted to these 
laboratory exercises : pair production in a static electric field, accelerated 
systems which become spontaneously excited in Minkowski-space and accel- 
erated mirrors. Black hole evaporation is then presented in the light of these 
examples. 

The common leitmotif that runs through all of this will gradually emerge 
upon reading. The two central elements are : 

1) The existence of a horizon in each case, a separation of space-time into 
portions, which would have been in causal contact in other circum- 
stances, that are no longer so. Straight trajectories become curved so 
as to approach or recede from the horizon exponentially slowly as seen 
by an external observer. The result of this law of approach is a modi- 
fication of the spectral decomposition of the quantum matter fields in 
such manner as to give rise to a steady state of production. 

2) The mechanism of production is the steady conversion of vacuum fluc- 
tuations into particles. A large part of the analysis is devoted to the 
description of how a particular fluctuation gives rise to a particular 
production event. 

In this way, we reveal the mechanisms at work in black hole evaporation. 
These are understood in terms of usual quantum mechanics. For the black 
hole problem they turn out to be surprisingly complex in that they involve 
a subtle interplay between the conversion of vacuum fluctuations to physical 
particles and a restructuration of the vacuum itself. Each gives its contribu- 
tion to the energy momentum tensor that provides the source of gravitational 
feed-back. Some of this richness is already present in the semi-classical ap- 
proximation wherein the source of Einstein's equations is given by the mean 
energy-momentum tensor. But this approximation will certainly break down 
at later stages, obviously at the Planckian stage, but possibly before, even 
much before. What we have set out to do is display the structures at work in 
the hope that this approach will sharpen the focus on the types of problems 
which will have to be confronted at later stages. The hope is that, from 
thence, quantum gravity wifl emerge. But, as introduced in one BBC radio 
program, that was next week's news! 
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Considerable effort has been put into exhibiting the configurations of field 
energy which contribute to the production mechanism. This is the content 
of Section ll.4[ 12.5.31 12.61 and 13.51 Particularly relevant for black hole physics 
is Section I2.5.3[ Of necessity, the analysis is not without subtlety (though 
rigorous and mathematically straightforward). The reader, who wants to get 
into the midsts of the black hole problem immediatly, may choose to skip 
Sections 12.5.31 and 12.61 in a first reading. For his benefit Section 13.51 contains 
a summary of all the necessary information related to the configurations of 
the fluctuations. The reader may continue on from this point (and of course 
then return if he wishes to the demonstrations and more detailed conceptual 
discussion of Sections 12.5.31 and 12.61 

One word on presentation. Generally speaking, the type of physics dis- 
cussed is more conceptual than technical in character. As such the mathemat- 
ical formalism is usually quite simple, unencumbered by extensive algebraic 
manipulation. Nevertheless, from time to time, the algebra does become 
heavy. When we have judged that the details of the mathematics interrupts 
smooth reading, we begin the topics under discussion with a resume of the 
results and relegate details to a passage enclosed within square brackets, or 
to an appendix. 
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Chapitre 1 



Pair Production in a Static 
Electric Field 



1.1 Qualitative Survey 

It was shown early on by Heisenberg and collaborators [50] that the vacuum 
of quantum field theory in the presence of a static electric field is unsta- 
ble against the creation of charged pairs. Subsequently using techniques of 
functional integration in the context of the action formalism, Schwinger [86] 
showed in an almost one line proof that the overlap function between vacuum 
at the time the field was turned on with vacuum at subsequent times decays 
exponentially fast. 

After Hawking discovered black hole radiation [H] several authors [67] . 
[87] , [15] , [73] pointed out the analogy with this kind of pair production. Though 
the analogy is not strict the formalisms which are used to encode these two 
phenomena have many points in common thereby making the electric pair 
production problem a fruitful exercise. To introduce the problem we first 
give a brief sketch of the physics in terms of a given mode (solution of the 
wave equation). 

We shall work in one dimension with a charged scalar field. This contains 
enough physics to make the exercise applicable to the black hole problem. 
For definiteness we choose the gauge At = Ex, = (with the charge 
absorbed into the definition of E) whereupon the Klein-Gordon equation is 

□(^ = ^+dl - {dt - iExf] ifi = mV • (1-1) 
Since the electric field is static, in this time independent gauge the operator 
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dt commutes with the d'Alembertian operator □, and modes can be taken of 
the form ip^it^x) = e~*'^*/t^(a;) where fuj{x) satisfies the equation 

+ {uj + Exf] Ux) = m^Ux) (1.2) 

Upon dividing eq. (11.21) by the factor 2m, one comes upon a non relativis- 
tic Schrodinger equation in the presence of an upside down oscillator po- 
tential centred at the point Xc = —uj/E. The eigenvalue for the effective 
energy, is — m/2. Therefore the solution tunnels between two turning points 
{xc ± m/ E = Xc ± a~^) where a is the classical acceleration due to the field. 

Take the case E > 0. Then classically a particle is uniformly accelerated 
to the right, following the classical trajectory {x — Xc)^ ~ {i ~ tcY = fl"^- So 
if it comes in asymptotically from the right along its past horizon {x — Xc) = 
— {t ~ tc) it will turn around at x = Xc + and t = tc- In quantum 
mechanics this is translated into the following description. A wave packet 
centred for example around u = (i.e. Xc = 0) will come in from the right 
asymptotically at the speed of light, slow down and turn around at x = 
and fly off to the right approaching once more the speed of light along its 
future horizon. The incident packet is localized around the classical orbit 
x^ — (t — tc)^ = Near the point x = +a~^ there will be some amplitude 
to tunnel through to x = —a^^ and then continue from that point on to 
accelerate to the left. Since E > 0, the tunneled particle must therefore 
have been "mesmerized" into an antiparticle in the tunneling region. This 
situation can only be met by second quantizing ip (confronting the Klein 
paradox [56]). The "mesmerization" is simply reflecting the fact that there 
is a probability amplitude to create a pair in the tunneling region. Then 
the above scattering description must be amended to: the initial particle is 
scattered as classically, but its final flux after the scattering is increased by a 
factor to accommodate for the creation of an antiparticle on the other side. 
The flux of this latter is denoted by Thus we see that |a;p — |/3p = lis 
the statement of charge conservation. 

In Fig. O 
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Fig. 11.11 The function In \ip{t, x)\ is plotted for an incoming wave packet 
solution of eq. (1.1)- It is a minimal packet as described subsequently in the 
context of eq. (1.54). We use the logarithm because of the very small 
amplitude of the produced particle (^-'^'^l'^°- where we have taken m/a = 9). 
The particle (on the right ) is centered on its classical trajectory 

— = a~^. The wave function of the antiparticle behaves classically for 
t > and X < —a~^ and vanishes in the past. 

we have plotted the modulus of a wave packet. One sees both the classical 
path about which the packet is centered and the created pair. Note that 
the separation between the turning points ±a~^, the length of the region 
from which the pair emanates, coincides with the distance |Ax| necessary to 
have the electrostatic energy compensate the rest mass (i?|Aa;| = 2m), so 
as to make possible on mass shell propagation from these points outward. 
It is the realization of this possibility, without energy cost, that causes the 
vacuum instability. Note also that since the amplitude of probability for a 
vacuum fluctuation to have its particle-antiparticle components separated by 
a space-like interval, |Aa:|, is exp(— m|Ax|), it may be anticipated that the 
amplitude for production is ~ exp [— C(m/a)] where C is some constant 
of 0(1) to be calculated by the formalism, and found to be equal to vr. We 
have thus shown how a wave packet not only describes the classical orbit but 
gives an idea of the location of the produced pair as well as its probability 
amplitude. 



12 



1.2 Mode Analysis 

We shall start with the analysis of eq. (11.21) considered as the Schrodinger 
equation of a fictitious non-relativistic problem. The solutions are well known 
[99] . We will nevertheless present them by a method [73] which is suitable 
for subsequent use in that it displays the singularity of the solution along the 
horizons. It is this singular behavior which in every case (pair production 
in an electric field, accelerating observer and mirror, black hole) encodes the 
production phenomena that is peculiar to quantum physics in the presence 
of horizons. 

The function solution of eq. (II. 2p is the x representation of a ket |H), 
i.e. fuj{x) ={x\'B). Go over to dimensionless variables centered at the origin 
by introducing 

C = VE{x + u;/E) (1.3) 
so as to write (11.21) in the suggestive form 

^^^m) = -em) (1-4) 

where 

e = m^jlE = m/2a (1.5) 
[tt,^] = —i (vr = —id/d^ in ^-representation) 

Now introduce the new canonical variables [5], analogous to the standard 
annihilation and creation operators used to quantize the harmonic oscillator, 
defined by 



so that (II. 4p becomes 



[U,V] = -I 
m)= Jdum{u\E) 



(1.6) 
(1.7) 
(1.8) 



2 

We can go back to ^-representation through 



where we have adopted the notation U\u) = u\u), defining the u representa- 
tion whereupon V = idu- To find the kernel {S,\u) of this unitary transfor- 
mation one must solve the conditions 



V -u 


u) = 










du 




u) = 
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to yield 
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exp 



:i.9) 



where we have supplied the norm to make the transformation unitary. 
The general solution of (11. 7p is 



K{u) = {uiE) = A e{u) 



u 



ie-l/2 



+ B 



-u 



|ie-l/2 



-u 



;i.io) 



The constants A, B follow from initial conditions. Consider the classical 
orbits associated with eq. (11. 4p . described by the Hamiltonian problem H = 
uv. They are given by u = UQe~'^,v = v^e^ with UqVq = —e. At early times 
T the orbits located near u = — oo, v = 0^ (i.e. incoming from ^ = +oo) are 
described hj A = and the ones near u = +oo, v = 0~ hj B = 0. These two 
sets of solutions (labeled by e) are orthogonal, and complete. In the effective 
Schrodinger problem which we are discussing their norm (with respect to the 
measure du) is such that /j!"^ duX^K' = 6{e-e'). This fixes |Ap+|5p = 1. 

Inserting each of the above solutions and (II. 9p into (II. 8p gives an integral 
representation for the scattered mode Reference [99j, identifies these 

as integral representations of parabolic cylinder functions, more succinctly, 
Whittaker functions Dy[z) (with v = —ie — 1/2) and provides their prop- 
erties and connection formulae. For our purposes however it is instructive 
to continue to work as "quantum mechanics" (see ref. [73]) and show how 
to relate asymptotic amplitudes to on-shell quanta. The effort is worthwhile 
because of remarkable analogs with the modes that arise in black hole physics 
(the reason being that the singular character of the solution eq. (1 1.1 OP near 
the origin encodes the existence of a horizon). 

Take the solution (11.101) with B = 0. We then must analyze 



du 



23/4 



-u 



ie-l/2 



TT 



exp t[^y2 + V2^u + uy2] 



(1.11) 
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-eir/i i-iv/ 



23/4 



-e --e^'"/^r(l/2 



TT 



in the limit ^ 



±cxD. The saddle points of the integral lie at 



(1.12) 



For ^ ^ -oo, the saddle at = -^2^ + ^/(v^O 

may be used to evaluate 

the integral 11.121 by the saddle point method since it lies well within the 
limits of integration and one has therefore the reliable asymptotic estimate 
(W.K.B. approximation) 

1 



2vr|el 



: exp —I 



eV2-eln(y2|e|)-vr/4 



;i.i3) 



Since the solution that we are describing (i? = in eq. (11.101) ) corresponds 
to an incoming particle coming from the left, hence with Udassicai oo as 
^ —oo for the incident beam, eq. (11.131) is to be identified with the incident 
wave function (/). For ^ — > +oo, the saddle at u* = — + e/{\/2^) gives 
no contribution since it lies completely outside the limits of integration. 

The other saddles at = e/ {V^$,) cannot be used to evaluate the integral 
by saddle point integration since, for ^ — >• ±oo, they lie at the edge of the in- 
tegration domain. However we expect that, since classically, v +oo for the 
transmitted solution and v —>■ —oo for the reflected one, one might fruitfully 
exploit the f-representation. This is indeed the case. The f-representation 
of \E) is 



{v\E) 



+ 00 



{v\u) {u\E)du 



+00 g™" yie-l/2 



'2t^ V27r 



-du 



gi7r/4g-e7r/: 



Bigl[9{v) 



-is-l/2 



2n 



_|_g-i7r/4ge7r/2 



/27r 



r(l/2 + ze) 



;i.i4) 



= T 



9{v) 



- 



R 



-ie-l/2 



d{-v) 



V2n 

-ie-l/2_ 
^ - 



(1.15) 
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This result follows from standard analysis, most easily by deformation of the 
contour from the real axis to the positive (resp. negative) imaginary axis 
for positive (resp. negative) values of v. It is this difference that gives the 
different weights R, T of the reflected to transmitted amplitudes, the ratio 
of which is 

lfl=e^^ (1.16) 

This exponential ratio is generic to all problems having horizons since the 
behavior of the wave functions at the horizon gives rise to the singular form 
ti**^, hence analytic continuation as exhibited in ll.l5[ 

Since we have normed the unitarity condition + |Tp = 1 is 
fulfilled by eq. fll.lSp . as is easily checked thanks to the formula |r(l/2 + 
ie)\'^ = 7r/cosh(7r£:). 

The identification of the waves multiplied by R and T as the reflected 
and transmitted waves, thereby giving a precise expression to the physics of 
the tunneling process, follows from the asymptotic behavior of (■ClH). Indeed 



T 



+R 



m{v\E)dv 

+00 ^i[-e/2+V2(v~v^/2+7T/4] 



23/% 

gj[-5V2+v^5t'-fV2+7r/4] 



1 



23/4 

g£7r/4g-i7r/8 



TC 



le 



1-^l^\f2T^ coshvre 

+Z},,_i/2(v^e^'^/^0 



A.-i/2(v^e-='-/^0 



:i.l7) 



;i.l8) 



The phase of {i\v) ( = -^V2 + \f^iv - + n/A) is fixed by the relation 
(^1"^)= / du. The saddle points are now located at 



± V2e - £ 



:i.i9) 



The saddles that go to ±00 as ^ ^ ±00 are responsible for the transmitted 
and reflected solution respectively. For these the integration can be estimated 
by a saddle point approximation to give 

^ = T = , (1.20) 



5-^+00 



v/2ie 
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= R . (1.21) 

The other saddles {v ±0 as ^ — >■ ±00) are concerned with the inci- 
dent solution which has already been analyzed in w-representation. This is 
fortunate in that the width of these saddles extend outside the domains of 
integration in eq. (11.171) . We call attention to the complementary roles of the 
saddle points in u and i;-representations in isolating the asymptotic incident 
wave on one hand and the transmitted and reflected parts of the solution on 
the other. As previously mentioned, this is due to the fact that the locus of 
saddle points is along classical trajectories. 

1.3 Vacuum Instability 

The above solution and notation describes the effective Schrodinger wave 
function which solves (11.21) . We now analyze how it encodes vacuum insta- 
bility and pair production. As we have announced the transmitted wave must 
have opposite charge. This can seen by following the movement of physical 
wave packets, solution of the Klein-Gordon equation (II. ip . These are of the 
form 

duof{uo-Uo)e-'^^'~'''\,{y/E{x + uj/E)) . (1.22) 

Note the novelty of the construction wherein the eigenvalue, tu, appears in 
the argument of the function Xe(0 = (^1'=') since ^ = \fE {x + uj/E). It is 
this circumstance that allows for the correct charge and current assignments 
of each asymptotic part of the wave (11.131) . (I1.2ip . (11.201) : 

1 e-^^(^+'^/^)'/2|a; + ^/^|imV2i?-i/2 =x {x^ -00, t +00) 

--t^ /?e^^(-+-/i?)V2|^ + ^/^|-imV2i?-l/2 =7^(3;^_oo,t^_oo)(1.23) 

Te'^^''+^/^^'/\x + cj/E)-'^™'/2^"i/2 =r {x^ +00, t +00) 



The motion of each branch is given by the group velocity obtained by set- 
ting the derivative with respect to u of the phases equal to zero. Thus, X 
corresponds to motion to the left at late times, whereas TZ describes motion 
to the right at early times and T to the right, but at late times. The roles 
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of X and TZ have been swapped once we consider the dynamics with respect 
to the physical time appearing in the Klein-Gordon equation fll.ip instead of 
the fictuous dynamics associated to the solutions of the Schrodinger equation 

Charge assignments follow from the sign of the charge density 

jt = ^ ^*^S = K +2Ex)^^ = 2{u: + Ex)x*eXe (1-24) 

Thus X and TZ which live at large negative x have negative charge whereas 
T has positive charge. In summary, for E > {] the solution of eq. (11.101) with 
i? = represents an incoming anti-particle of amplitude i?, and outgoing 
reflected anti-particle of amplitude 1 and a transmitted particle of amplitude 
T. The unitarity relation + |Tp = 1 is then best recast into a form out of 
which one reads charge conservation rather than conservation of probability. 
Divide this relation by \R\^ and recast it into the form 

-|a|2 + |/5|2 = _i 

|«P = l/|i?|2 ; = |T|V|i?|2 (1.25) 

The right hand side then corresponds to incident flux of unit negative charge 
coming from x = — oo, is the outgoing transmitted flux of positive charge 
going to +00 and (lap — 1) is the increase of flux of the reflected wave 
necessary to implement charge conservation. 

For each value of u, in addition to the above mode, there is its parity 
conjugate, obtained from the transformation [x + uj/E) —>■ —{x + uj/E) and 
complex conjugation. Clearly it represents the charge conjugate since one 
need only run through its mirror image about the axis Xc{= —u/E). It 
corresponds to the solution of eq. (II. 4p proportional to 9{v). The modulus 
of a wave packet built of this latter mode is given in Fig. II. 1[ 

From these considerations, it follows that these basis functions for the 
quantized field, whose quanta are single particles in the past (formed from 
wave packets), lose this property in the future due the production. More 
specifically let us envisage the situation wherein E is switched on during a 
time interval < t < T, over a length L. Then the modes fall into two 
classes: 

• Class /: Those that scatter {j3 ^ 0) i.e. whose turning points lie within 
this space-time domain of area LT . 
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• Class //: Those which are only deflected, so as to follow the classical 
trajectory before they reach the would-be turning point (in the sense 
of wave packets). For these latter P = 0. 

Indeed, we shall see subsequently that, for modes in the first class, if a 
pair is formed, it will emerge from a region of area 0{a~^ x a~^) whereas 
in the second class no pair is formed to within "edge corrections" . By this, 
we mean for L » a^^ and T >> a~^, no pair is formed except near the 
edges. Thus to obtain asymptotic results of 0{L x T) for the production 
of pairs, it is a legitimate idealization to divide the modes into these two 
classes. For the remainder of this chapter, all considerations are devoted to 
those in the non trivial {(3 ^ 0) class /. For these modes, the future fate of a 
single particle mode in the past is to become a many particle mode, hence not 
convenient to describe the result of a counter experiment in the future. Thus, 
in addition to the modes previously presented - denoted by "in-modes" - we 
shall be obliged to introduce "out-modes" as well. These latter will represent 
single-particle quanta in the future of their turning points. 

It is in the framework of second quantized field theory that this discussion 
makes sense. Therefore we begin by quantizing the field operators in the in- 
modes. The complete set is given by 

^7Ji^^) = (4^)1/4^. - uj/E) (1.26) 

^Zit, x) = ^^^|i/4^ e— *x.(x + uj/E) (1.27) 

Note that for completeness, u should span the range [— oo, +oo] contrary 
to the case without E-field where the energy has to be taken positive only. 
The wave functions have been normed according to the Klein-Gordon scalar 
product 

Jdx^TJt,x)tV,cpi-*it,x) = -5{uj-uj') . (1.28) 

Here the label p designates particle and a anti-particle (described by the 
parity conjugate as previously discussed). The reason for the complex con- 
jugation in the l.h.s of eq. fll.27p is that the function is a solution of the 
charge conjugate field equation (here obtained by changing the sign of the 
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charge, i.e. E — > —E in the field equation). Then yj^"^ and y^™* obey the 
same equation and the second quantized field is written in the in-basis, 

Ht,x) = /ci^(a>;^(t,a:)+CVaX(^'^)) (1-29) 
with the usual commutation relations 

= 6{uj-uj') (1.30) 

The out-basis is obtained simply by observing that if one follows a single 
outgoing branch backwards in time it will trace out the same backwards 
history as that of one of the in-modes when it moves forward in time. These 
modes at fixed u are obtained by sending t into —t and complex conjugating 
so that the new functions remain solutions of the field equation fll.lj) . 

= ^Tj-t.^) (1-31) 

Furthermore, since for each u, the two sets of in-modes eqs fll.26| 1.271) are 
complete, these out-modes must be linear combinations of them. In fact they 
are given (after choices of phases) by 

r (1-32) 
- (1-33) 

where a and (3 are the coefficients introduced in (11.251 ll.25p . now fixed 
precisely as: 

(3 = T/R , a = e'^'l^/R (1.34) 

To understand this result it suffices to follow the histories of the y)™'s as 
wave packets [16]. In eq. (I1.32p . for example, yjp"^ will give rise to a refiected 
wave of amplitude a* on the right and a transmitted wave of amplitude (3* on 
the left, whereas (p^^*^ gives a transmitted wave of amplitude j3 on the right 
and refiected wave of amplitude a on the left. The contributions on the left 
therefore cancel and those on the right have total amplitude |ap — (= 1) 
as required for a particle out-mode. The assiduous reader will check that the 
right hand side of eqs (11. 32111. 33p are indeed the complex conjugates of the Xs 
functions defined in the in-modes (all multiplied by e*"^*). The easiest route 
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to obtain these relations is to use eqs fll.lllll.17p in the context of integral 
representations for 

XsiO = e"^ [Txli-O + (1-35) 

and the correct identification of these wave functions as used in second quan- 
tization. 

One develops $ in the out basis by inverting eqs (11.321 11.33P , 

in „ * „out I n* ntj 

^p,uj = « ^p,u;+P ^a, 

in* „ out* I n „out 

and substituting into eq ll.291 to give 

$ = / rfc.(ar <j + CV?) (1-37) 

where 

C*^ = + (1.38) 

These equations fll.38p define a Bogoljubov transformation [lOj wherein eq. f ll.25p 
ensures that a°"* and 6°"* obey the correct commutation relation. Historically 
it came up when Bogoljubov [12] noted that the free particle Bose-Einstein 
ground state was unstable against creation of pairs of equal and opposite 
momentum once interparticle interactions were introduced. The interested 
reader will find a brief account of Bogoljubov's considerations in Appendix lAl 
From f ll.38p it is clear that the in- vacuum |0, in) (i.e. the state annihilated 
by in-annihilation operators) will contain out-particles, ie. we set up this 
Heisenberg state at early times and count the mean number of out particles 
that are realized from it at some later time. 

(n^) = (0, in|a°"*ta-*|0, in) = \(3\^ . (1.39) 

One can express the in-vacuum state as the linear combination [55] of out 
states since the Bogoljubov transformation is unitary. 

|0,in) = N-^/^exp 



(1.36) 



a 



/ J UJ UJ 



|0, out) 



:i.4o) 
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This result is obtained by setting |0,in) = f{a°^* \ ^]fut t^jg^out) and impos- 
ing that = a^"|0,in) = (aa'"^*-/36°"* t)y|o,out) = (a[a°"*, t/)|o,out). 
This in turn imphes that a{df /da°^* ^) — /96°"* ^ = 0, hence that / has the 
form written in ll.40[ Note that the creation operators a°"* \b^^ appear 
only as a product. Hence to each produced particle in mode p, a;, there cor- 
responds one and only one antiparticle in mode a, uj. From what has been 
described these particles are born in pairs of opposite values of their con- 
served quantum numbers such as charge and energy . Thus the energy of the 
states containing pairs spontaneously created is equal to the energy of the 
in-vacuum. This is readily seen by expressing the hamiltonian of the field in 
terms of the a°"*, 6°"* operators. One finds 

/oo 
duuial^'^a-f - (1.41) 
-oo 

To compute the normalization factor N in equation I1.4UI we have 

(0,in|0,in) 

N 

From this it follows that 

I (0, out|0, in) p = 1 = n ^ = exp - ^ ln(l + ^ (1.43) 

thereby delivering the probability to find no pairs at future times. 

To complete the calculation we must now give a meaning to the X^o; i-^- 
we must count the number of modes in the non trivial class (those which 
turn around in the interval < t < T and < x < L). 

The density of orthogonal modes of frequency uj in the interval < t < T 
is ^du. To calculate / du, we note that the scattering centers of non-trivial 
modes are at the points x = —uo/E where < x < L. Therefore the 
total number of non trivial modes is {T /2tt){EL). The values of uj which 
contribute to / duj can also be obtained from classical mechanics. Indeed 
the classical equation of motion is md'^t/dT'^ = E dx/dr which integrates to 
rndt/dr = Ex + uj in the region < x < L. The turning point is given 
by dt/dr = 1 (i.e. hj uj = —E{x — a^^) ). Therefore the modes that 
have turning points in the region < a; < L have energies in the interval 
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a\ 



:i.42) 



—EL < < 0, so that for L » a ^ the number of modes in class / is 

Of course the same result can be obtained in 

the gauge = —Ft, At = 0. 

From this discussion a physical picture of production emerges wherein a 
certain class of vacuum fluctuations (class /) in the past contains the po- 
tentiality of making pairs. This comes about because these fluctuations fall 
into resonance with a state which contains pairs of out-quanta. The pair 
production is possible because of the degeneracy of states of zero energy (see 
eq. (11.411) ). In each case that we review in this article the same physics is 
repeated. Only the details of the mechanism of conversion from vacuum to 
physically propagating states changes from problem to problem. Moreover 
the selection of those fluctuations that are predestined to become physical 
particles is such as to give rise to a steady rate of production, the common 
ground being the constancy of acceleration. 

The Schwinger formula [86] now follows from (11.431) 

|(0,out|0,in)|' =exp-:^ln(l + e-'"'"/^) (1.44) 

where we have used (I1.16P = |T/i?p. Actually Schwinger worked out the 
case of 3 -|- 1 dimensions. This is obtained by replacing by the transverse 
mass squared (= -|- kj_) and integrating over k±{J (fk±/27i) in the expo- 
nential so as to replace ETL/2'k by E'^TV / At^"^) . For the interested reader 
Appendix [B] contains a modified version of Schwinger's derivation obtained 
by functional integration. The mode by mode analysis given here is obvi- 
ously far more detailed in revealing physical mechanisms and hence of more 
pedagogical value for the black hole problem. 

The physical interpretation of eq. (11.441) is illuminating when one recalls 
that the mean number of quanta produced in the mode uj is = {nj) (c.f. 
eq. (11.391) . The argument of the exponential in eq. (I1.44p is thus Y.^j^'^i)- + 
(n^^)), as in a partition function. 

It is interesting to remark that the population ratio of produced pairs of 
two charged fields one of mass m and the other of mass m + Am is given, in 

^It is amusing that it is precisely the same consideration that affords a very simple proof 
of the chiral anomaly of Fermi fields in 2 dimensions in the presence of a magnetic field H . 
In its Minkowski version, replace H hy E and chirality (75 = 7172) by velocity (= 7071)- 
Then the change in velocity of the vacuum due to E {= J dxtpjoliip) is (l/27r) / Edxdt. 
See for example [SS] 
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the case Am « m and m? / E « 1, by 




27r Am/a 



(1.45) 



{n)m 



= e 



i.e. the ratio of the mean number of created pairs with neighbouring masses 
is Bolzmannian with temperature a/27r. That the physics of accelerated 
systems is associated with this temperature is the subject of the next chapter. 
From 11.451 and Section [21 we can say that particles are born in equilibrium. 

1.4 Pair Production as the Source of Back 
Reaction 

Following [16] let us examine more closely the physical situation that arises 
when the E field is turned off after a finite time lapse {E ^ 0; in the interval 
< t < T), wherein the (Heisenberg) state is the vacuum |0,in). This 
state, at time T, the content of which is then measured by a counter system 
(i.e. a measuring device sensitive to out-quanta), is expressed in eq. fll.40p . 
It appears as a linear superposition of different outcomes according to the 
number and nature (mode number) of produced pairs. For example eq. f l 1.441) 
gives the probability to find no pairs and eq. fll.39p gives the mean number 
of pairs in a mode. These pairs produce a non vanishing expectation value 
for the electric current at intermediate time t: 



where the second equality results from the insertion of a complete set of states 
in the out-basis. The coefficients (= (m, out|in)) are the probability 
amplitudes to find the system in the state |m, out) at time T. For a system 
with many pairs, the state |m, out), is very complicated due to interactions 
among these pairs. The whole development of the previous sections has been 
devoted to the production of non interacting pairs and we shall continue to 
work in this approximation valid for sufficiently small times T or large masses 
m. More precisely we require 



(J^(t,x)) 



(0,in|J^(t,a;)|0,in) 





m 



ET 



—nm' 



'/^[a-^+T] « 1 
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This inequality ensures negligible overlap among the pairs produced in time 
T, hence negligible interactioiu. Equation fl 1.471) comes about by multiplying 
the density of pairs (= (ET/27r)e-™'/^) by their mean separation (= a ^ + 
T) since they are produced at separation whereupon the members of the 
pair fly apart at the speed of light. 

The opposite limit of large density, where mean field approximation is 
valid has been the subject of detailed investigation [22] where the back re- 
action has been solved. Unfortunately many body effects which give rise to 
plasma oscillations are an important element of this development. Hence it 
is not a relevant laboratory for the study of black hole radiation. Therefore 
we shall continue with the approximation of non interacting pairs, which 
turns out to be more relevant to the study of black hole evaporation. When 
eq. (11.471) is valid, the relative fluctuations of density of pairs are large, so 
the mean can no longer be used. It is no longer appropriate for calculating 
the effect of back reaction on the production. For this reason it becomes of 
interest to take apart eq. (11.461) term by term. 

For the case of negligible interaction the normalized states |m, out) are 
direct products over pair states : 

\„OUt\'LOUt] 
- . /.I /.I 



Since is quadratic in the field operator $, it contains creation and an- 
nihilation operators in a double sum over modes. It is easy to check that 
the diagonality of the Bogoljubov transformation (11.381) not mixing different 
values of cj, reduces the double sum to a single one. Thus the in- vacuum 
expectation value (J^) in eq. (11.461) is expressible as a single sum over ui. It 
is given by 



which may be found by either hard work from eq. (11.401) or using the identity 
for the Feynman propagator: 




(1.48) 



(J^(t,x)) 



(0, out| J^(t, a:;)|0, in) 
(0, out|0, in) 






(O,out|<l>(t,x)<l>t(t',x')|0,in) 




^Recall that the pairs are one-dimcnsional dipoles 
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= (0, m|$(t, x)$t(t', x')|0, in) + ^ ^^^(t, x)^;Z(^', x') 

(1.50) 

easily obtained by expressing in the in-basis, and expanding the out- 
vacuum state (O,out| up to the term quadratic in the in-annihilation opera- 
tors acting on the in-vacuum : 

(O,out| = A^~^/2(0,in| exp 

The first term on the r.h.s. of eq. (11 .49^ is the current contained in the 
state |0, in) which arises when no pairs are produced. Hence it should vanish. 
It does. There are two proofs. The first one is related to subtraction problems 
in general. This method will be presented in Section 13.31 (see ref. [16] for an 
explicit treatment in the case of a background electric field) . The other is to 
prove that upon expanding into modes, each term is separately zero. We 
leave this as an exercise to the reader (see also ref. [16j). This latter proof 
is completely satisfactory since the summation over the alternating series is 
absolutely convergent due to the finite domain of space-time L x T in which 
the field is non vanishing. 

The second term on the r.h.s. of eq. (11.491) is the mean current, expressed 
as a weighted sum over the contribution of the current carried by pairs arising 
from the mode u. The weight is the probability that at least one 

pair be produced in this mode (since l/|ap is the probability that none be 
produced). The c-number in parenthesis which multiplies this weight will be 
called a weak- value. It is related to what one measures in an experiment in 
which the presence of a pair is detected, thereby making contact with the 
weak measurement theory of Aharonov et al. [3]. In our situation, the post- 
selection (i.e. the detection of a localized pair) is physically performed by 
registering the click in a counter. More generally, whether or not the pairs 
are treated as independent, one may write (J^) as (see eq. (I1.46P 

m 

where Pm = is the probability to find the system in the state |m,out). 

In other words the expectation value of is a weighted sum of non diagonal 
matrix elements (weak values) in exactly the same way as usual conditional 
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E 



in itn 
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(m, out| J^(t, a;)|0, in) 

Am 



;i.52) 
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probabilities. The physical relevance of these matrix elements as well as the 
concepts of post selection and weak value are discussed in more detail at the 
end of this section and in Appendix O 

For the nonce, we shall enquire into the properties of these matrix el- 
ements using well localized wave packets To this end a Gaussian envelop 
turns out to be most convenient for the following reason. First consider the 
asymptotic regions. For simplicity we take a mode which is centered on the 
space-time origin. Then, see eqs ( 11.24111.261) . we obtain : 



t- 



2 



^ ^-iEx^/2^-ix+ty/2^ 



lim - / rfwe-^'/^a^e-*"* l^^+iEix+u./E)y2 ^ ^*^-^Eix+./Ef/2 



where 



1±1 
cr2 E 



:i.54) 



and inessential phases factors have not been taken into account. The asymp- 

r 1-1/2 

totic widths are given by Re(l/S|) . If Ima 7^ the two outgoing 
branches in (11.531) will have different widths. Charge symmetry then dictates 
the choice Imcr = whereupon once sees that the width (= [[E"^ + a^) / E'^a'^]^^'^] 
is minimized by the choice o"^ = E. In what follows, wave packets of such 
width will be qualified as minimal. They are not only optimally localized, 
but they constitute a good approximation to a complete orthogonal set. This 
follows from the fact that their width is ~ E~^^'^ about the classical orbit so 
that the number of such (asymptotically) non overlapping packets is 0{ELT) 
as is required for the correct count of non trivial modes. With a little tin- 
kering on their size and shape they could be shaped into a rigorous complete 
orthogonal set [SS], but for our purposes the minimal packet is sufficient to 
show how the back reaction field emerges from a single pair. 

The subsequent analysis is facilitated by using the integral representation 
(II. lip which from (II. 9p is a Gaussian transform of {u\x)- Hence with ^ given 
by (ll.3p the Gaussian integral over uj is trivial and one finds once more that 
the packet is nothing but another Whittaker function. With cr^ = E, one 
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obtains 



-oo 

^ eW4r(l/2 + ^^) g^(.^+.^+.2x.t)i./4^^^^_^ r^W4 /^(^ _ ,^)] 

(1.55) 

and ^/'p" is given by the same superposition in terms of x*e see eq. fll.26p . 
The algebraic details are given in ref . [16] . We have plotted in Fig. 11.11 the 
(logaritm) of the modulus of ip^^^- 

The current (11.491) has been written as a sum over modes, this can be 
rewritten as a sum over packets. The contribution to this sum from ^/^^ and 

|a|V;(t,x) = ^i;;-*it,x)tV,i;rit,x) (1.56) 

(The factor |ap on the l.h.s. is introduced here for subsequent physical 
interpretation, see eq. (I1.58P ). We emphasize that this current is a product 
of two localized wave packets. It is a contribution to the non- vanishing second 
term of eq. ( 11.49P when the set of modes is recast into a set of packets. 

In Fig. ll.2l we have plotted the real and the imaginary part of Note m 
particular that this product vanishes in the remote past since we start from 
vacuum at early times. Indeed each factor {ip^^^ and ip^J^ respectively) of this 
product is a minimal packet centered on the space-time origin and they do 
not overlap before they reach their common turning point. Being in- modes, 
one comes in from the left only and the other from the right only. This is at 
is should be, the current produced by the pair should not exist before it is 
produced. 
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(a) 



(b) 



Fig. 11.21 The real (Gg. a) and imaginary (Gg. h) parts of Jgf are plotted for 
a minimal wave packet where we have taken m/a = 9. In the figure zero is 
grey, positive is white and negative is black. Within the production region 
At X Ax = 2a~^ X 2a~^ the real and imaginary parts are comparable and 
oscillate. Outside the production region, the quanta are on mass shell and 
propagate classically: the imaginary part vanishes and the real part takes 
its classical value. 



Moreover, in the far future, from the inverse of eqs. fll.32[ I1.33P taken 
together with the vanishing overlap of '?/'°"*and (where is the same 
wave packet as but made of out-modes) we find 



hm 

t— >+oo 
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This once more is at it should be. The future current is the sum of currents 
carried by the particle and anti-particle separately according to one's classical 
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expectation. Remark on the importance of the denominator in eq. fll.491) (or 
more generally in eq. (11.521) ). wherein {m,out\ is given by (l,?/'""*! the out 
state which contains only the pair described by and V'a"*- Hence it is 
given by 

(l,?/'°"*|0,in) al3*^^ ^ ' ' ^ ' > m v > 

It is the denominator a(3* that garantees that the current carried by the 
detected out pair is unity. The multiplicative factor |ap in eq. (11.571) is there 
to account for the fact that more than one pair may be produced in the 
packet. It is equal to the "induced emission" factor (n+ 1) where (n) = 
is the mean density number of pairs produced. It is not present in eq. (11.581) 
since in that case only one out pair characterizes the out-state. 

Having understood that asymptopia goes according to rights, we note 
further features seen in Fig. 11.21 

1) The region of production lies within a circle of radius 2/a, 

2) Within this region the contribution to is complex and it oscillates. 

This is the quantum region. The oscillations occur because for x and t <C 
the modes oscillate with frequencies of 0{m). This is seen most clearly in 
the gauge where At = 0, = —Et. Then the packets are built out of 
X-e(t + k/E)exp(ikx) and thes^ oscillate near the origin with frequency 
y/m"^ + k^. The amplitudes of these oscillations is very strong near the origin 
(like expE{a^ — — x^)) and then fade out as the particles settle down to 
get on to their mass shell. We are seeing how the transients work themselves 
out so as to arrive at the completion of a quantum event. 

Note however that in the sum for (J^) eq. (11.491) . the imaginary part 
vanishes identically and furthermore, since E is constant in the box (J^(a;)) 
is homogeneous, so that the oscillations drop out as well. Thus one should 
question the physical relevance of these complex oscillations. We therefore 
conclude this section with a few remarks on this count, including a short 
discussion of the back reaction problem. 

The simplest way to put them in evidence in a physical amplitude is by 
considering the change in the probability to find the ip pair upon slightly 

Which are the analytic continuation from to —m? of the functions introduced in 
eq. dLll). 
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modifying the electric field, i.e. by replacing E hj E + SE{x,t). Since 
the change in the action is given hj 6S = J dtdxJf^6A^, the change of the 
probability to find the ip pair is given, in first order in 6E{x, t) by 



Pe+5e = \{l,r''y\0,m)\^ 

= |(l,^°"*|(l + 255)|0,in)|2 



PeIi + 2 J dtdx6A^'{t, x)lm 





J, 


0,in) 







,in) 



;i.59) 



(We recall that Pe = | (1, '?/^°"*|0, in) p = \l3/a\'^/N where is given by 
eq. fll.42p ). Hence it is the imaginary part of the weak current (given in 
eq. fll.SSp ) which furnishes the change of the probability when one compares 
two neighboring iJ-external fields. But given that each created pair carries 
an electric field these weak values will also govern the modification of multi- 
pair production due to current-current interaction in the usual perturbative 
approach. One way to proceed is to include the interaction term in the hamil- 
tonian | / / dxdx' Jt{x)v{x — x')Jt{x') where v is the Coulomb potential. It 
can be shown that once one pair is produced, the electric field which acts on 
fluctuations to give rise to subsequent pairs is of the form E + j dx' jf{x') 
where jf{x') is the matrix element of the charge density due to the pair 
which has been created. The self interaction of a pair as it is produced is 
given by a more complicated loop correction. The statement of the back 
reaction problem is thus given in terms of a self consistent highly non trivial 
problem. One will then be led to treat the complete wave function of the 
system as a function of the configuration of both the matter and electric 
field. This will then give a generalization of the usual semi-classical theory 
[d^F^'^ = (inlJ'^lin), wherein F^" on the l.h.s. is classical, but modified due 
to the mean of J*^], to a full quantum Heisenberg equation among operators 
d^F^"" = r. 

The hard job then remains. How do these individual production acts 
affect subsequent production. This is the true feed-back problem and it has 
not been solved. Current research is now under way to find an approximate 
solution which is better than the mean field approximation (see for instance 
ref. |26j) wherein fluctuations and correlations among pairs are taken into 
account. 

Another way to understand the status of the matrix elements eq. (11.581) 
and the relevance of the complex oscillations is by appealing to the formalism 
developed by Aharanov et al [3]. In this formalism, it is shown how those 
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matrix elements give the result of a "weak measurement" . This development 
contains two elements that we briefly expose. The interested reader will 
consult Appendix [C] as well. 

• A weak detector, a quantum device in weak interaction with the system, 
where weak means that one can correctly work in flrst order in the 
interaction. The detector wave function is spread out so that flrst 
order perturbation theory is valid. Then one picks up not only the 
detector's displacement (in position) as in the von Neuwmann model 
of measurement [HS] , but one flnds also that the imaginary part of the 
matrix element ( here of J^) imparts momentum to the detector. For 
the case we discuss here the detector could be for example a test charge 
whose path is deflected by the external E fleld, modifled by the back 
reaction field (induced by the weak value of the current) due to the 
selected pair, i.e. the modification of the E field is given by Gauss' 
law: AE{x,t) = dx'jf{x',t), where jf{x',t) is the weak value of 
the current. Its statistical interpretation in the context of measurement 
theory is given below. 

• A post selection, the specification of the final state of the system, here 
the state |m, out). In the context of our wave packet development with 
no interaction among the pairs a possible out state is the one post- 
selected by the click of a localized counter. One selects, thus, that part 
of the wave function |0,in) which causes the counter to click. Physically 
this post selection consists of registering all results of the weak detector 
(as one usually does to establish mean values, in the manner prescribed 
by the Copenhagen interpretation of quantum mechanics) and then to 
keep only the results when the particular localized click is registered 
(as one does in constructing conditional probabilities). 
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Chapitre 2 

Accelerating Systems 



2.1 The Accelerated Detector 

A system in constant acceleration through empty Minkowski space, coupled 
to a field whose state is quantum vacuum in Minkowski space, will heat up 
and its internal degrees of freedom will become thermally distributed with a 
temperature given by 

T = l3-^ = a/2'K . (2.1) 

More succinctly : an accelerated detector perceives the vacuum as a thermal 
bath. This remarkable observation is due to Unruh [91], who thereby gave 
substance to an equally remarkable observation by Fulling [35] who showed 
that quantization of a field in Rindler coordinates [STj is inequivalent to the 
usual Minkowski quantization. 

The definition of Rindler coordinates p, r in 1 + 1 dimensions is 

t = p sinhar , x = p coshar 

p > , -oo < r < +00 (2.2) 

hence applicable to the quadrant t > 0, x > 0, hereafter referred to as the 
right quadrant (R). These coordinates are naturally associated to a uniformly 
accelerated system with acceleration a in the sense that its trajectory is 
p = = const; its proper time is r. In addition the Minkowski metric 
expressed in coordinates p, t 

ds^ = -dt^ + dx'^ = -p'^a^dr^ + dp^ (2.3) 
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is static, translations in r (boosts) being Lorentz transformations. 




Fig. 12.11 The coordinate system (p, t), the trajectory of a uniformly 
accelerated detector, the four quadrants of Minkowski space and the 
horizons H± (U=0,V=0) which separate them. 

The remarkable result, now stated with more precision is that Minkowski 
vacuum is described in R (or alternatively in the left quadrant (L) wherein 
t = —p sinhar, x = —p coshar, see Fig J2.1[ ) by a thermal bath in the 
quantization scheme based on the coordinates p, r defined by eq. f l2.2p . The 
connection between the remarks of Unruh and Fulling is through the proper 
local dynamics of the accelerator when it is coupled to the radiation field 
so as to undergo transitions. Conventional emission and absorption is then 
given in terms of Rindler quanta 

The easiest way to see this phenomenon is by looking at rates of absorp- 
tion and emission of an accelerated two level detector. One may for example 
take a two level ion whose center of mass is described by the wave function in 
the W.K.B. approximation to the minimal wave packet of Section The 
type of idealization envisioned is to let the mass of the ion tend to infinity 
at fixed acceleration, a, i.e. M —>■ oo, E —>■ cxd, E/M a. The packet then 
describes a 6 function along a classical orbit whose center is used to define 
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the origin of space-time 



\ t = a~^sinhar 



X = a 



^coshar i'^-'^) 



(Acceleration is constant in the sense that a^a^ = where a'^ = d^x^ /dr^ = 
a^x^ with r the proper time on the path of the accelerator, i.e. u^u^ = 
{dx^ I dT){dx fj, / dr) = —1.) The two levels are separated by the mass difference 
AM, taken to be finite. 

We take the example of a massless scalar field (here chosen hermitian) in 
1 + 1 dimension, coupled to the two level accelerating detector [H] , [2B] , [SI] • 
In the interaction representation, the interacting hamiltonian of the system 
is 



dk 



9 



e (T+ + e (T_ 



(2.5) 



where a± are the operators which send the system from ground (excited) 
state to excited (ground) state and where = u = \k\ ,ki = k and al (a^) 
are the creation (destruction) operators of field quanta of momentum k. 

Standard golden rule physics gives in lowest order perturbation theory 
the following formula for the rates expressed in terms of the proper time r, 
of absorption (-) and emission (+) 

/ + 00 
dAr e^'^^^'^^WiAr - le) (2.6) 
-oo 

where W is the Wightman function between two points on the orbit separated 
by proper time At. For the reader unfamiliar with this formula we sketch a 
brief derivation 

Take the case of absorption, then we begin at r = with the atom in its 
ground state and the radiation field in vacuum. The amplitude to find it in 
the excited state is the matrix element of —i Jq Hint{T')dT' . The probability 
is obtained by squaring the amplitudes and by integrating over all k and is 
therefore equal to 

/ rdT2 fdr, e-^*^(---) [(OmI^ (x^(r2)) ^ {x^{n)) |0m)] (2.7) 
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where |0m) is Minkowski vacuum. For emission change the phase e~'^^('^2-n) 
to e+*^^^(^2--ri)^ fjj^g quantity in brackets is called the Wightman function 
W^(t2, ti) and the positivity of the Minkowski frequencies (u) encoded in its 
decomposition into annihilation and creation operators impose an integration 
prescription in the complex plane given by ie in eq. (12. 6p 

The Wightman function for the case of a massless field in 1 + 1 dimension 
is — (47r)~-^ In [(At — ze)2 — Aa;^)] where for constant acceleration eq. (12. 4p . 
we have 

(At — ieY — Ax^ = (sinhar2 — sinhari — ie)^ — (coshar2 — coshari) 

= 4a"^sinh^a(Ar -2e)/2 (2. 

Therefore for this C3iSG, clS well as for the case where the path is inert ial 



const), the integrand in eq. (12.71) is a function of T2 — ti onljQ. Take 
T in eq. (12. 7p such that r >> Max(AM~^, a^^). [The usual condition to 
establish the golden rule is r >> AM^^ (for spontaneous emission) and 
T » (3 (for absorption of photons from a thermal bath). In the case of 
uniform acceleration this implies r >> a~^] One may change variables to 
(r2 + ri)/2 and T2 — ti = Ar. Since the integral comes from the finite region 
around T2 — ti = 0(AM~^), the integral over the latter may have its limits 
extended to ±cxd. (For (r2 + ti)/2 < AM a small error is made which does 
not contribute to the rate.) Then integration over {t2 + ri)/2 gives r thereby 
yielding a rate formula per unit propertime. The rate is given by eq. (12. 6p . 
The reader unfamiliar with this formalism will check out the usual golden 
rule in terms of density of states for the inertial case. Also he will relate 
eq. (12. 6p to the imaginary part of the self energy. 
Substituting for W in eq. (12. 6p we have 



sinh^ f a ie 

2 



±ig^— / dr -, — ^ ^ (2.9) 

47ri-oo AM sinh f§r-ze' ^ ' 



where we have integrated by parts accompanied by the usual prescription of 
setting infinite oscillating functions to zero. We remark that in 3 + 1 dimen- 
sions the same type of formula obtains with W given by a^sinh~^(ar/2 — ie). 



^ The origin of this miracle is that both kinds of trajectories are orbits of the Lorentz 
group. In the accelerated case r translations are boosts. 
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From eq. fl2.9p . it is seen that for the minus sign (case of absorption), 
complexifying r and extending the domain of integration to a closed contour 
in the lower half complex plane picks up poles at ar = 2Tmi {n = —1, —2, . . .). 
For the plus sign (emission) the contributing poles lie in the upper half plane, 
at ar = +ie and ar = 27ini {n = 1,2, . . .) whence the ratio 

p Y^oo -f3AMn 
^- _ 2^n=l ^ -/3AM 



e-^^^"^ (2.10) 



with j3 = 2Tx/a. Since R- is proportional to {n) and to {n) + 1 one 
obtains, following Einstein's famous argument, (n) = (e^^*^ — 1)"^ as in a 
thermal bath. Explicit evaluation yields 

= ^il-e~^^'r' (2.11) 

Note that these rates coincide exactly with the rates given in an inertial 
thermal bath for the model considered. This occurs for massless fields in 2 
and 4 dimensions only. 

Equation (12.101) . in fact, results from a general property of the propagator 
W appearing in eq. (12.61) related to the periodicity of the orbit eq. (12. 4p when 
T is imaginary, hence transforming the orbit to a circle. To see how this is 
related to thermal properties we first present the thermal Wightman function 
for a general system wherein H may include interaction of the field with 
itself 



W^thermal — tr 

= tr 
= tr 



e-^^0(ti)0(t2)J (2.12) 
e-'^^0(At)0(O)] (2.13) 
e-^('^-'^*V(0)e-'^^V(0)l (2.14) 



where eqs. (l2.13]) and (I2.14p follow from 0(t) = e*^*0(O)e~*^*. Because of the 
positivity of the spectrum of H (i.e. En — Eq > 0) , WthermaK^T) is defined 
by eq. (I2.14p in the strip in the complex plane 

-/5 + e < ImAr < -e (2.15) 

In this strip one has the identity 

W^thermal(i) = W^thermal(-«/5 " t) (2.16) 
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obtained by using the cyclic invariance of the trace. Under general conditions 
Wthermai IS analytic in the strip. [For a general discussion of thermal Green's 
functions see ref. [37]]. 

Analyticity and eq. fl2.16p then imply eq. (12.101) . This is derived by in- 
tegrating §dz e^'^^^Wthermaiiz) = ovcr the contour surrounding the strip 
f l2.15p . and we assume that there is no contribution from the ends at Re^ = 
±00. The integral along Im^; = — e is then —R+ and along Imz = —(3 + e is 
e^^^R- thereby recovering eq. fl2.10p . 

The point of all this is to note that W for the accelerating system enjoys 
the property fl2.16p since it is a function of (Ax^y = — 4a~^sinh^aAr/2 
and this is true not only for a free field but in general. Thus the ratio 
of the accelerating detector rates of absorption and emission is given by 
eq. f l2.10p simply in consequence of the imaginary periodicity of the orbit. 
We emphasize that it is not necessary that W{At) be equal to an inertial 
thermal propagator eq. fl2.12p and indeed for free field theory one may check 
that apart from d = 2 and = 4 it is not since the density of Rindler modes 
of energy AM differs from the inertial one. However eq. 02.160 and therefore 
eq. f l2.1Up always obtains. 

2.2 Quantization in Rindler Coordinates 

We now show how it is possible to interpret these results in terms of annihila- 
tion and creation of Rindler quanta (Rindlerons) in situ. Rindler coordinates 
are defined by eq. (12.20 in the first quadrant and it is in this quadrant where 
we situate the accelerator. We make this point somewhat more explicit. 
The transformation eq. (12. 2p is analogous to euclidean polar coordinates, 
with the role of the angle being played by ar (hence Imr has period 27r/a). 
However whereas euclidean space is completely covered by polar coordinates, 
Minkowski space is only partially covered in one quadrant: x > 0, x > \t\ 
designated in Fig. 12. II by R. 

The d'Alembertian —d^+d^ is [—p~'^a~'^d^+p^^dppdp] so that for massless 
particles the modes in these coordinates are solutions of 

' 1 



(2.17) 
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In R, a complete set of eigenmodes of idr, solutions of eq. fl2.17p . is 

A > 



taken together with their complex conjugates. Here 



(2.18) 



u 

V 



TTa~^lnap (2.19) 



The modes are normalized according to the Klein Gordon norm which for u 
modes can be written as 

/ + 00 <— > 
duip*y^jiidu ipx,R = 5iX-\') (2.21) 

and similarily for v. We note that even though u, v are defined in R only by 
eq. (12.191) . by trivial extension (px.Fj{u) is also defined in P since u is finite on 
the past horizon V = 0. Similarly (p\fi{v) is defined as well in F as well as 
R. 

In L, the corresponding complete set of modes of positive Rindler fre- 
quency is 

VxA-^l) = e{U)e+'^^^/y/^ = 0(f/)(at/)-^^/VV4^ , 

(2.22) 

taken together with their complex conjugates where in L we have 

t = -psinhar 
X = —pcoshar 

and Ul, Vl in L are also given in terms of r, p by eq. fl2.19p . Hence in L 
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Once more, by extension fx^Liui) is valid in L and F and (px^iyi) is valid 
in L and P. Note that in L, dt/dr < whereas in R, dt/dr > 0. Therefore 
the modes fl2.22p have Rindler frequency opposite to the usual one. The 
convenience of the convention (12.231) is that a given boost is represented by 
the same displacement in r in both L and R. In addition the mapping of R into 
L by reflection through x = 0, t = is very simply described by the analytic 
continuation r — >■ r + in/a, u —>■ u + tTr/a = Ul, v —>■ v + in/a = which 
maps (O into (12:231) and fl2:20D into (12:21) respectively. If one performs 
this inversion twice one obtains that R is invariant under r — >■ r + 2m/a 
which is the essential ingredient used to obtain the thermal properties in the 
discussion following eq. (12.111) . 

These sets of modes are the Rindler versions of the Minkowski light-like 
modes 



TTUJ 



= e-'^^/V^ (2.25) 

and their conjugates. 

The advantage of the use of these light-like variables is that the right 
movers (functions of u or U) and left movers (functions of v or V) do not 
mix under Bogoljubov transformations (e.g. (P\^r{u) is a linear combination 
of ^Lu{U) only). For the remainder of this section we shall work with the right 
movers only. The left follow suit. 

We have two representations of the right moving part of the field 

POO 

ct>{U) = / duj[a^UU)+^.^] 
Jo 

roo 

= ei-U) / dX [cx^RifxMiU] + h.c] 
Jo 

POO 

+e{U) rfA[cA,m,L(f/)+h.c.] (2.26) 
Jo 

To find the Bogoljubov transformation in R we express ^uj{U) as linear com- 
bination of (P\^r{U) 

POO |- . 

e{-U)UU)= d\ \atvx,R{U) + Pt^lAU)] (2.27) 



«L = / dU^*{U)idu UU) 







40 



du (fXniu) idu ^^{U{u)) 

1 fO 



2tt y/uj 



dU{-aUy 



iX/a—l^—iivU 



27ra V uj 



-T{-i\/a)(- 



T\/2a 



UJ 



r+oo ^ 

and = / du (px,R{u){-idu)^u;{U{u)) 

J — CO 



27ra V 00 



-T{i\/a) - 



UJ 



iX/a 



-nX/2a 



(2.28) 



In terms of operators eq. (12.261) gives the Bogoljubov transformation 



CX,R 



du 



(2.29) 



From the exphcit values (I2.28p . or more simply from their integral represen- 
tations, or more formally from the canonical commutation relations of Ca,r 
and c\ ^ one checks out the completeness of 



duj 



duj 



R* oR 



6{\ - X') 
. 



(2.30) 



Of course, the set ipx^R is complete only in R. For complete completeness its 
complement v^a.l, as expressed in eq. (12.261) . is required. Then the operators 
a^^ may be expressed as linear combinations of ca,_r and ca,l. So the totality 
of Bogoljubov transformations is eq. (12.291) . its analog in L (wherein one has 
the simple rule a^^ = af* and = Pxu)y inverse given by 



dX 



a 



cx,R + ailcx,L - l3tc\R - PxA,L 



R J 



X(jj 



(2.31) 



To calculate the Rindler content of Minkowski vacuum we use eq. (12.291) 
to obtain the mean number of Rindler particles in Minkowski vacuum 



(Oa/IcI^ca'kIOm) 
n{X) 



du^PtP^r^ = 5(A - X'MX) 



(2.32) 
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where we have used 12.251 and |r(ix)p = vr/xsinhvrx and n{X) is the Planck 
distribution. One finds also 

roo 

{OmIcxrCx'rIOm) = / diuatP^r^ = (2.33) 
Jo 

In this way one understands the thermal propagator W of eq. 02.61) . Indeed 
since W has been calculated only in R, the accelerator perceives Minkowski 
vacuum in R (to which it is confined) as a thermal bath (of Rindlerons since 
in the proper frame of the accelerator it is they which are in resonance with 
the excitations of the accelerator). 



2.3 Unruh Modes 

There is an interesting and powerful technique of handling the Bogoljubov 
transformation 12.281 which is due to Unruh [HI] . It couples together Rindler 
modes which lie in R and L respectively (</?A,i? and v^a.l) so as to yield a 
density matrix description for what happens in one quadrant when one traces 
over the states in the other. This technique makes contact with the physics 
of pair production in a constant electric field and finds important use in the 
black hole problem. It is implemented by inverting 12.271 (using 12.301) to give 

^UU) = / du; [a^lUU) - P^lCiU)] (2.34) 

Though the inversion has been carried out in R, eq. fl2.34p can be extended 
into L (since the ^^(f/) are defined there as well). It will be checked out 
below that this is perfectly consistent since we will find automatically that 
in this continuation the right hand side of eq. fl2.34p vanishes in L. We rewrite 
eq. (12.341) in the form 

^x,R = axVxiU) - PxV-xiU) (2.35) 

where 

1 1'^ 

0xiU) = — dua^lUU) 



ax JO 
1 



VUU) = - I du CiU) (2.36) 
Pa -'0 
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and 

1 

Px = , ^ (2.37) 

so that the modes (p\ are normed in the usual way, i.e. \_'^dlJ(p\{JJ) idu 
0x'{U) = S{\ - A'). Substituting eq. (12:281) into eq. fl236D gives 

Oa 27ra V 47r Jo uo \aj 

{ axe{-U){-aUy^/''/V4^ + /3xe{U){aUy^/''/V^ A>0 
j /3|A|^(-f/)(-af/)^^/VV4^ + a|A|^(f/)(af/)^^/VV4^ A<0 

I aA<^A,i? + /5a<^I,l a > 

1 AA|<^|A|,iJ + a|A|<^|A|,L A < 

We shall call the modes, 0x, Unruh modes. As announced the linear combi- 
nation fl2.35p does vanish in L. 

Unruh modes enjoy the following properties: 

1) They are eigenfunctions of iaUdu with eigenvalue A in both R and L. 

2) (fx are manifestly positive Minkowski frequency modes for both signs 
of A (c.f. eqs f l2.36p ). Together with their conjugates they form a complete 
orthogonal set just as plane waves. This is proven trivially by direct compu- 
tation. 

3) They are linear combinations of the Rindler modes ip^ and ip^ given 
by the Bogoljubov linear combination fl2.38p (we remind the reader once 
more that the mode W^^''9{U),X > is a negative frequency mode in L 
(eq. (^M)- 

An independent and interesting derivation of eq. fl2.38p is obtained by ap- 
peal to analyticity in the lower half U plane. This is because the Minkowski 
modes eq. (12.250 are analytically extendable in the lower half plane for 
positive frequencies only. This expresses the stability of the ground state 
(Minkowski vacuum) of the theory. Eigenfunctions of idu considered as func- 
tions of U obey the differential equation 

- ^U^uXx = -XX (2.39) 
a 
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having solutions 



XX = A 9{-U) 



a 



+ B e(u) 



(2.40) 



(Note the similarity to eq. (11.101) . This is no accident. Any problem with 
an exponential approach to a horizon in the classical theory will be reflected 
in such a singular differential equation when expressed in terms of global 
coordinates.) 

To determine A and B in eq. (I2.40p we require that x\ be a positive 
frequency Minkowski mode i.e. an Unruh function. Therefore set U equal 
to \U\e'^ with -TT < < so as to continue the function f/*-^/" analytically 
in the lower half complex plane to go from L to R. The boundaries at 6* = 
— TT, then give the ratio B/A = e^'^^l"- = e~^^/'^ (compare with eq. fll.l6p . 
Normalization then gives x\ = 'Px- This technique has been fruitfully used 
by Unruh [91] and Damour and Ruffini [27] in the black hole problem. 

The Unruh modes can then be synthetically written as 



with the ie encoding how the function should be continued in the complex 
plane {U — ie = U for U > and U — it = |?7|e~*'^ for U < 0). This expression 
for (fx has the added luster that for small but flnite e the high frequency 
behavior of the modes in the vicinity oi U = has been regularized [7T] . 
This will be seen to have important consequences when evaluating energy 
densities near and on the horizon, see Section 12. 6[ 

The quantized fleld can be decomposed in Unruh modes everywhere 
according to 



The bogoljuobov transformation among Rindler and Unruh annihilation and 
creation operators is 




(2.41) 




(2.42) 



dx = a\Cx,R - (^\Cx,L 



, A > 



(2.43) 



where we have used eq. fl2.38p . 
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In this way Minkowski vacuum can be viewed as a linear combination 
of pairs of L and R Rindlerons. Defining \ORindier) as the direct product 
|0_r)|Ol) of Rindler vacuum in the left and right quadrant (ca,l|Ol) = and 
ca,_r|0_r) = 0), we have following eq. fll.40p 



|OAf) = ^-'/'exp 
where Z is given by (see eq. (11.421) ) 



.A=0 «A 



A t t 



10 



Rindler 



) 



(2.44) 



ni"Ai 



exp 



POO 

/ dX ln(l + Pi) 
Jo 



L 1-°° 
2^ 







= exp 


12 /3_ 



(2.45) 



which is the partition function of a massless gaz in 1+1 dimensions in a 
volume of size L. Thus for an operator Or localized in R one then has 



(Or) 



Z (Offi„dZer|eXp 
2nA 



p. 



{nx} 



X 



ax. 



V — cx,rCx,l Or exp 

.t^Q "A J 

{{nx}\On\{nx}) 



P. 



.£"A 



^X,R^X,L 



I ^ Rindler^ 



{nx} 

tre-^^«Oij 



(2.46) 



where Hr is the Rindler hamiltonian generator of r translations restricted to 
R. So single quadrant operators have their means given by a thermal density 
matrix. Note that the Rindler hamiltonian (the generator of boosts) is equal 
to —ialldu and is given in terms of the operators cx,r, cx,l by 



H 



Rindler 



Hp — Ht 



/+00 p + oo 

du Tuu - e{ + U) / dUL 

r+oo 

/ dX\{cl jiCx,R - c\ lCx,l) 

JO 



(2.47) 



and possess therefore the same structure as the hamiltonian in the E'-field 
given in eq. (11.411) . Thus the pairs of Rindlerons in eq. (I2.44p have zero 
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Rindler energy. This degeneracy allows for the creation of pairs of Rindler 
quanta. 

Of course everything that has been derived for u modes applies equally 
for V modes where the role of P and F are interchanged hence the past and 
future horizons H± (see fig. 12.11) . Also one can perform the analysis for 
fermions and include the effects of a mass and of higher dimensionality (see 
for instance [BH], [53]). 

It is more interesting (and relevant for blackholology as well) to inquire 
into the physical condition of the radiation due to its interaction with the 
accelerator. We shall therefore in the next section delve more intimately into 
the details of the transition amplitude. 



2.4 Spontaneous Emission of Photons by an 
Accelerated Detector 

We here dissect the rate formula eq. (12.111) by displaying the transitions 
of the atom as resonance phenomena with Doppler shifted Minkowski pho- 
tons [74j. This will introduce in a natural way the conversion of a particular 
Minkowski vacuum fluctuation into an on mass shell quantum. Similar res- 
onance phenomena constitute the dynamical origin of particle emission from 
a non inertial mirror (Section 12.51) and of black hole radiation. 

We shall see that in this formulation the rate formula comes about by a 
somewhat different mechanism from that of usual golden rule analysis, i.e. in 
terms of density of states. It is rather a consequence of the steady process in 
which photons are brought into resonance through the ever-changing Doppler 
shift occasioned by the accelerator, much like the analysis of Section 11.31 
(paragraph following equation (I1.43P ). wherein different u values brought 
different vacuum fluctuations into resonance with real pairs giving rise to a 
time sequence of produced pairs. 

We focus for definiteness on right movers. In lowest order in g, the am- 
plitude for a transition in time r is 

A^{u;,t) = -t{±\{OMK r dT'H,r^t{r')\OM)\T) 

Jo 

\/4:TlUJ Jo 

where the plus (minus) subscript of A corresponds to spontaneous deexcita- 
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tion (excitation) and the ket |±) means excited (ground) state. Begin with 
the former, spontaneous deexcitation. The integrand presents a saddle point 
at T* given by 

AM = cje-""*^'^) . (2.49) 

If u is such that t*{uj){= Inu/AM) lies well within the domain of inte- 
gration in eq. 02.481) (i.e. the gaussian width = (AMa)~^/^ does not overlap 
the limits of integration), and furthermore if higher derivatives of the ex- 
ponential are relatively small (i.e. (AM/a) >> 1), then the saddle point 
estimate is valid and one finds 



The crucial inequality for legitimization of eq. (12.501) is that u should be 
bounded by 

0<a~Mn-^<r (2.51) 

AM - ^ ^ 

but we shall shortly loosen up on the condition (AM/a) >> 1. 

Once more there is, as for the electric division into class I, the 

frequencies uj which resonate, i.e. satisfy eq. fl2.5ip and class II modes which 
do not resonate, i.e. lie outside the range (12.511) and for which A^{uj, r) ~ 0. 
And once again as in Section 11.31 this must be qualified by " apart from edge 
effects" . Clearly r has to be sufficiently large for these asymptotic estimates 
to make sense, as in all "golden rule" type estimates. 

For spontaneous excitation the saddle point condition has a minus sign 
on the r.h.s. of eq. (I2.49P so that the saddle point lies at 

Rer*(cu) = a~Hnuj/AM 

Imr*(u;) = vr/a (2.52) 

whereupon one has 



= -A;(cj,r)e-^^^^/' (2.53) 

with (3 = 211 /a. Equation (12.531) is valid only for class I modes; for class II, 
A-{uj, r) ~ as well. Squaring these formulae one recovers the thermal ratio 
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of rates gven in eq. fl2.10p . Integrating uj over the bounds set by eq. 02.511) 
and dividing by r, one recovers half of R± given by eq. (12.111) provided one 
hmits oneself to the contribution from the leading poles in the integrand 
(12.91) [ r = — ze for emission and 27rz/a for absorption, consistent with the 
condition AM/a >> 1]. The other half is due to the left movers. Note how 
the rate formula comes about, e.g. 



The integral over a; is / din a; hence an integral over saddle times Rer*(a;). 
Neglecting edge effects this is equal to r. The physical interpretation is 
clear. At each time t*[io\ the Minkowski frequency lo enters into resonance 
because of the changing Doppler shift occasioned by the acceleration (cuRes = 

It is interesting to try to interpret the complex saddle encountered in 
eq. (I2.52P (the case of absorption). In the complex r plane the contour 
has to be deformed from Imr = so as to go through Imr = vr/a (as 
in the construction of Unruh modes-paragraph after eq. (12.391) . In both 
cases the voyage in the complex plane encodes the positivity of Minkowski 
frequencies). Reference to eq. (I2.20p then indicates that there has been a 
voyage from \J to — f/ (or x —>■ —x, t — > — t), hence from a point in the 
quadrant R to its inversion in the third quadrant L where r runs backwards 
(i.e. dr/dt < 0) hence where Rindler energy conservation is satisfied since 
the frequency uj appears as carrying negative Rindler energy (since it is given 
by iaUdu see eq. (I2.47p . In Sections 12.61 and 12.6.31 it will be shown by more 
detailed considerations that in fact the Minkowski photon which is emitted 
when the atom is excited arises from a vacuum fluctuation one part of which 
exists in the quadrants L and F, the other part lives in P and R. This latter 
has crossed the past horizon so as to be absorbed by the atom and the former 
continues out to infinity in quadrant F, i.e. it arises in L and radiates from 
there into F. It is this effective emission act in L which is encoded in the the 
rough Born approximation saddle point integral of saddle point eq. (I2.52p . 
The reader at this point can pick up some flavor for the true state of affairs 
by peeking ahead at Fig |2.5[ 




2AMa 




(2.54) 
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One can offer oneself the luxury of completing this type of analysis to get 
the full expressions for R± (i.e. by relaxing the condition AM/ a » 1) by 
noting that for any value of AM /a the contribution to the integrand from 
a given a; is dominated by the region (gaussian width) around t*{uj). So 
provided one maintains eq. fl2.5ip one can extend the r' integration over the 
whole real axis. In this limit the transition amplitudes eqs (12.481) and (12.531) 
become 



A- = -Wsj^Pluu (2.55) 

thereby giving a dynamical content to the Bogoljubov coefficients (I2.28p . 
Squaring and integrating over uj within the bounds limited by eq. (12.511) 
leads to the exact rates given in eqs (12.531) multiplied by the interval r. This 
is how the Golden rule comes about in this approach, successive resonances 
with the Doppler shifted frequencies of the proper atomic frequency (here 
= AM). 

For the skeptical reader we now present a more rigorous analysis of this 
argument as it plays an important role in the black hole problem as well. At 
the same time we shall clearly exhibit the difference between the class I and 
class II modes. 

[ The amplitude to emit a photon of frequency uj in the interval (— r/2, r/2) 

is 

A^^^^r) = -j^(-) / dx e-'^x'^''/^-' (2.56) 



^/A^^UJ \(lJ JLja-le-"^/2 
—in /,,,\-iAM/a p~AMTT/2a 

= ^(-) 7(^AM/a,.c.a-e-/^) 

--f{iAM/a,iuja-^e-''^/^)] (2.57) 

where the interval is taken symmetric around r = for mathematical con- 
venience; the physics is unmodified by this translation in r and where 7 is 
the incomplete gamma function [99] which for large and small values of its 
argument takes the form 

7(z/x,zx) - r(z/i) - ze^'^^/V^e-^^a;-^ (2.58) 



'y{ifi,ix) - -ie-'^^l^x'^ (2.59) 
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Before analyzing equation fl2.57p we discuss the unphysical infra-red di- 
vergence which arises as u ^ 0. This divergence exists in the inertial case 
as well. Indeed as lo —>■ one has yl+ = 2ig sin( AMr/2) / y/ Attlj independent 
of the acceleration. It is unphysical because in one dimension the coupling 
of the massless field to the atom becomes strong as u —>■ 0. Therefore the 
perturbation theory which has been used is inadequate. No doubt a resum- 
mation of all terms is possible to give the true infra red physics (like the 
Bloch Nordziek [TT] theory in QED). But this is irrelevant to our purpose, 
since these infra red photons do not contribute to the rate. 

One way to get rid of the problem is to subtract off the inertial amplitude. 
We shall do the computation in another manner so as to obtain an infra red 
finite answer. This is accomplished by adiabatically switching on and off the 
coupling to the field by a function /(r). Then the amplitude for u —>■ will 
decrease as a function of the period T of switch on and off. To see this we 
remark that 

1) for uj 0, is given by the Fourier transform of / {A^{uj ~ 0, r) ~ 

/dre-^^^V(r)/v^). 

2) If the interval r during which / is constant is equal to r = 27rfc/AM 
with k an integer then Aj^{uj ~ 0, r) is independent of k. 

3) Therefore we can take k = 0. 

4) If / can be differentiated n times, the Fourier transform of / decreases 
for large AM as AM"". 

5) Since the only dimensional parameter is T, the Fourier transform of / 
must be of the form T-^+^AM"". 

The correct procedure to calculate the rate of excitation consists in first 
taking the adiabatic limit T —>■ oo (with however the condition T << r ) 
and only then performing the integral over uj. We shall show that with this 
order of operations the concept of resonant frequency described in the main 
text appears. 

To perform thus the adiabatic limit it proves rather more convenient to 
introduce a time average of the amplitude 

A+{uj, t)=J dr' g{T')A+{uj, r') (2.60) 

where g{T') is a n times differentiable function, centered on r' = r, with 
width T, and normalized such that Jdr gij) = 1. One verifies by permuting 
the integral over r in eq. (12.601) and the integral over x in the definition of 
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(eq. (12.561) that the averaging of is equivalent to the introduction of 
a switch function /(r) with the same regularity as g. 

We now permute the series expansions eqs (12.581) and (12.591) with the 
averaging in eq. (I2.60p to obtain that is given by a similar expression to 
eq. (12.571) with 7 replaced by 7 where (see eq. (I2.59P ) 

^{iAM/a^iLua'^e^"^/^) ~ T{iAM/a) + 

±aT/2 

— ^00 

a 

/ g-7rAA//2a \ 

+0 --r-—H] (2.61) 



/ ±iAMT/2 \ 

^iiAM/a^iuoa-^e^''^'^) ~ qI -.AM/2a ^ , , (2.62) 



We are now in a position to rederive the results obtained in the main 
text. Three cases are to be considered according to the values of uj: 

1) oj < ae~"^/^. Then both 7 functions in eq. (12.571) are given by eq. (I2.62p . 
hence their difference yields 

A+{lo,t) ~ 0[ae-^^^'^'^''sui{AMT/2)/iT''-^AM'')]. 

2) ae^''^/^ < UJ < ae°^/^. Then the first 7 function in eq. (I2.57P is given by 
eq. (12.6 ip and the second by eq. (12.621 hence their difference is T{iAM / a) + 

{uj, t) is given by eq. (I2.55P as announced. 

3) > a"^/^. Then both gamma functions are given by eq. (I2.6ip and 
y4+ vanishes once more. 

In Fig J2.2l we have plotted a numerical calculation of |7(iAM/a, ico'a"^e"'^/^) 
^{iAM/a, iuja~^e~"''^/'^)\'^ as a function of r and Intu. The plateau of this func- 
tion when — ar/2 < \nuj/a < ar/2 is clearly apparent. The oscillations on 
the plateau are due to the fact that the adiabatic limit has not been taken 
in this figure. 
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ln(co/a) 



Fig. 12.21 The square of the amphtude A{t, u) to have made a transition in 
time T by emitting a photon of Minkowski frequency uo. plotted as a funtion 
of T and Inuj. The plateau that arrises when the resonance condition 
— ar/2 < Inuj/AM < ar/2 is satisfied is clearly seen. 



In this way we liave confirmed tliat eq. (12.481) arises from tlie saddle point 
region to yield the result ( 12.55^ . Hence that the successive resonances do 
build up to yield a rate.] 

2.5 The Accelerating Mirror 
2.5.1 General Description 

A highly instructive chapter in the physics of accelerating systems is that of 
the accelerating mirror [36], |29j . [TO] in that the analogy to the production 
of Hawking radiation is remarkably close. 

We have shown that a uniformly accelerating system with internal degrees 
of freedom displays thermal properties due to the exponentially changing 
Doppler shift (eq. (12.491) ) which relates the inertial frequencies to its local 
resonant frequency. Similarly a non inertial mirror scatters modes with the 
changing Doppler shift associated with its trajectory, hence giving rise to 
physical particles. If the mirror trajectory is such that the Doppler shift is 
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identical to eq. fl2.49l) the particles emitted by the mirror will be thermally 
distributed as well. We write this Doppler shift as 



ue-''^ = k (2.63) 

where the energy difference AM is now replaced by the frequency of the pro- 
duced particle k, and u is, as before, the frequency of the Minkowski vacuum 
fluctuation in resonance with the produced quantum which is reflected from 
the point U on the mirror. In the next paragraph the order of magnitude of 
k is taken to be of 0{a/2Tr) whereas u varies strongly. We now determine 
the mirror trajectory which leads to this Doppler shift. 
The general solution of n4> = d"^ /dudv4> = is 

(t) = F{y) + G{U) (2.64) 

The reflection condition is that vanish on the mirror hence the solutions 
take the form 

(t) = F{y)-F{ym{U)) (2.65) 
where the mirror trajectory is expressed as 

V = V^iU) (2.66) 

To determine the Doppler shift associated with this trajectory we take 
an incoming Minkowski mode: e~*'^^/v47rci;. The reflected wave is then 
_^-iuiVm{u) I ^J^^^ This mode should be decomposed into the inertial outgo- 
ing modes e~^^^ / y/Ank so as to determine its particle content. The scattering 
amplitude (the Bogoljubov coefficient a^^k) is the overlap of the modes 

/pikU ^ „—iuiVm.{U) 
dU^==id u r— (2.67) 

The resonance condition given by the stationary phase of the integrand is 

U^J^ (2.68) 

Hence in order to recover eq. fl2.63p we must take as trajectory for the mirror 

Kn(f/) = --e-"^^ (2.69) 
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where we have set to zero an irrelevant integration constant. This trajectory 
tends exponentially fast towards the asymptote = which is the last 
reflected ray. It plays the role of a horizon in this problem (see Fig 12.31) . It 
is by construction that this mirror trajectory has the same expression as an 
inertial trajectory expressed in Rindler coordinates m, f as it approaches the 
future horizon u = oo (see reference [TU] for a mapping of the one problem 
into the other by a conformal transformation). 




Fig. 12.31 The mirror trajectory Vm{U) = —a ""^ which gives rise to a 
steady thermal Hux. 



The mirror therefore follows a very different trajectory from that of the uni- 
formly accelerated detector. 

In order to describe the particles emitted by the mirror it is necessary to 
work in the second quantized context. As in the electric fleld problem one 
introduces two bases. The initial one is given by 

if'j; = ^^(e"'"^ - e-*"^"'^"'''') uj>0 (2.70) 

and corresponds on the past null inflnity surface T~ {U = — oo) to the usual 
Minkowski basis. The in-vacuum is the state annihilated by the destruction 
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operators associated with these modes: 

0= / rfo; (a>r + h.c.) (2.71) 
Jo 

The final basis corresponds to the usual Minkowski basis on the future null 
infinity surface {V = +00), hence given by 

^outL ^ _^^^\^y\^k/ag^_y^ _ ^-ikU^ ^ ^ q ^2.72) 

V47r/c 

The V part of ^9°"*^ is identical to the left Rindler modes expressed in 
Minkowski coordinates (see eq. f l2.22p ). Similarly the U part of ip^^ is the 
same function as a Minkowski mode when it is expressed in Rindler coor- 
dinates. Hence the mathematics as well as the physics of Sections 12.11 to 
12.41 apply to obtain the number and type of particles emitted from the mir- 
ror. Indeed the alpha Bogoljubov coefficient (12.671) and the beta Bogoljubov 
coefficient 

dU^==td u _ (2.73) 

are identical to the Bogoljubov coefficients obtained in eq. (12.281) . In par- 
ticular the ratio IPujk/o^ujk] = e~'^^/°' obtains and implies a constant rate of 
particle production in a thermal spectrum of produced particles at temper- 
ature T = a/27r. 

The mechanism wherein one finds a constant rate of emission is the same 
as in Section 12.41 It arises from the successive resonances of in-modes of 
varying energy uj with the emitted photons of energy k. To a photon in a 
wave packet of energy k emitted around the value U = Uq corresponds one 
and only one \^-mode. Its frequency is = fce"^° as in equation eq. (I2.49p . 
And the number of photons of energy k emitted during a certain AU lapse 
is given by the integral of IPujk]"^ = {l/uj)n(k) ( where n{k) = {e^^'^ — 
over the resonant frequencies uj. This integral yields AUn{k), i.e. a thermal 
fiux times the time lapse, as in eq. (12.541) . 

We now calculate the energy momentum carried by these quanta. In the 
next subsection the mean energy will be obtained by two different techniques. 
In the third subsection we shall describe the fluctuations of the energy den- 
sity: we shall obtain the energy momentum correlated to the observation of 
an outgoing particle around a particular value of U. We wish to emphasize at 
this point that we are now seriously trespassing into the domain of the black 
hole problem. The next two sections contain a great deal of the essential 
physics of black hole evaporation. 
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2.5.2 The Mean Energy Momentum Tensor 

To start, from eq. (12 .701) we have 

(0m|Tv'y|0j„)|(7=_oo = 



(2.74) 



The surface {U = — oo) is where Minkowski vacuum is laid down. It 
is the Cauchy surface from which emanates the modes (I2.70p . This defines 
the Heisenberg state |0j„) which is annihilated by the operators a^. From 
eq. fl2.70l) the modes are pure ^-like on X~ (in the sense of broad wave 
packets) and are Minkowski modes. Thus eq. 02.741) is true in the usual 
sense of normal ordering. Furthermore since duTyv = for massless fields in 
Minkowski space it follows that eq. (12.741) remains true for all points {U, V) 
which lie to the right of the mirror's trajectory fl2.66p : 



{Oin\Tvv\Oin)\v>V„ 



0. 



(2.75) 



More interesting is (ruu)in- We calculate it in two different ways: 

1) by mode analysis 

2) more synthetically through use of Green's functions. 
In the first method we have 

{^in\Tuu\^in) — {^mit\Tuu\^cm.t) = 



dk duvT'duvl 



^outL* 



(2.76) 



wherein we have implemented the prescription of normal ordering by making 
a subtraction of the value of (Tuu) in Minkowski vacuum. Indeed by defini- 
tion the U part of the out modes eq. (12.721) is Minkowski in character. We 
use the Bogoljubov coefficients (12.671) and (12.731) to express in terms of 



.^outL 



dk akuj^l""^ + Pkuj^k 



out* 



Thus obtaining 



dk / dk' 

Jo 



{Oin\Tuu\^in) — {Oout\Tuu\^out) — 



(2.77) 
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dk duifT'^'duiplf''* 
The unitary relations (12.301) simplify the result to the form 



(2.78) 



dk / dk' 2Re 
Jo V\Jo 

We now recall eqs. fl2:32l [233D 

^-^{k+k')/(2a)Y(^ikia)T{-ik' /a)Vkl^ 



(2.79) 



dupkiuPl 



6{k-k')n{k) 



where 



n{k) 

One also verifies that 



^^/'^|r(iA;/a)pA; 
(27ra) 



(dc^/cj)(a/cj)'('=-'=')/" 
(2.80) 

(2.81) 



/•oo 

/ duiP, 
Jo 



Mk-k')/(2a)Y^^k/a)T{ik'/a)y/kk' 



{2TTay 



(rfcj/cj)(a/cj)'('=+'^')/" 
(2.82) 

is proportional to S{k + k') hence does not contribute to eq. (12.791) . The 
vanishing of this interference term is the expression in the particular case we 
are considering of the general theorem proven in eq. (I2.46P that expectation 
values of operators restricted to one quadrant are given by their average in a 
thermal density matrix (see also eqs (12.321 [2733]) ). This theorem is applicable 
here since the U part of each mode arises from the part of the mode to the 
left of the horizon {V < 0). 
The final answer is thus 



roc 

{Oin\Tuu\Oin) - {Oout\Tuu\Oout) = dk 2 n{k) du^T du^l 

Jo 



1 

2^ Jo 



dk k n{k) 



outL , ^outL* 
2 



n / a 
12 \2^ 



(2.83) 
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which is a thermal flux in one dimension with T = a/2'K as announced. 

The other technique to calculate the flux relies only on the trajectory of 
the mirror eq. fl2.69p [30] . We now pass it in review. 

In usual global null coordinates the metric of Minkowski space reads as 

ds^ = -dUdV (2.84) 

and the modes used to quantize about the usual Minkowski vacuum are e~*'^^ 
and e"*"^^. The U part of the reflected wave in eq. fl2.70p is in terms of modes 
Q-ii^fiu) f{U) = — a~^e~"'^. If we adopt f{U) in place of U as coordinate 
we have 

ds^ = -^dfdV (2.85) 

The U part of the Green's function in usual Minkowski vacuum is —-^ln\U — 
U'\ whereas in the modes e'^^^f^^^ it is -^ln|/(f/) - f{U')\. Thus the 
difference in {Tjju) = {du(j)du(j)) between the two is 

A{Tuu) = -^\imdudu'[\n\f{U)-f{U')\-\n\U-U'\] 
1 



12n 



f/^dlf-V^ (2.86) 



where /' = df/dU and where the second line is obtained by expanding / to 
third order in U — U'. In the present case f{U) = Vm{U) which yields 



{Oin\Tuu\0^n) - {Oout\Tuu\Oout) = = TT^T' {T = a/2n) (2.87) 

as required. The second line of eq. fl2.86p is a remarkably elegant formula 
which relates the average energy momentum to the coordinate transformation 
f{U) only. It does not refer to the presence of the mirror. It suffices to refer 
to the non inertial coordinates used in the expression fl2.85p of the metric. 
We shall see in Section (13.30 that there is a natural generalization to curved 
space which is of important use in the black hole problem. 

Note that the complete in-in Green's function involves terms in V, V and 
mixed terms U, V due to the linear combinations which appear in eqs (12.701) 
and (I2.72p . The term in V,V' gives (Tyv) unchanged from Minkowski vac- 
uum, eq. (I2.75p . The U,V term gives U,V energy-energy correlations but 
no contribution to (Tf^u) itself since the trace (4 {Tuv) = m?{(l?)) vanishes 
identically for a massless field in Minkowski space (see Section (13. 3p and ref. 

Iini)). 
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2.5.3 The Fluctuations around the Mean 



In preparation for black hole physics, we now ask a more detailed question, 
concerning fluctuations. We wish to display the field configuration in vac- 
uum that is responsible for the emission of a particular photon whose wave 
function is localized around U = Uq. [ Recall that the pure state Minkowski 
vacuum is a linear superposition of field configurations (it is the product 
of the ground states of the field oscillators: |0m) = Ile"''^"^'^'-" ^" in rep- 
resentation)]. In particular we shall describe the energy distributions which 
characterize this fluctuation (we consider this quantity since in the black hole 
problem the energy momentum is the source of the back reaction). . At the 
same time we will also answer the complementary question, to wit: what 
is the energy distribution when no photon is observed around Uq (since we 
know the average). The calculation is straightforward but somewhat lengthy. 
We therefore sketch below the main points and then go on to the proofs. 

The Heisenberg state |0j„) is a linear superposition of out states. It is 
a straightforward exercise to describe the contribution of each final state to 
the energy on X"*" (in order to do this in a local manner it is necessary to 
use localized wave packets). In addition to the positive energy density on 
T"*" associated to the production of a particle at U = Uq, there is correlated 
to it a "partner" which is a bump of field that propagates along the mirror 
(hence with > 0). Before reflection there is also a bump in the region 
V < 0. This is the "ancestor" of the produced photon. The total energy 
carried by this pair of bumps vanishes as behooves a vacuum fluctuation. It 
does so in rather subtle fashion. It is positive definite for \^ > whereas 
in the region V < 0, there is a positive energy bump, mirror image of the 
former as well as an oscillating broader piece which is negative. The sum of 
all these contributions vanishes on I~ and therefore for all U until the wave 
packet starts reflecting. 

If two (or in general n) photons are produced around Uq then the ancestors 
carry twice (n times) the energy if one photon is produced. If no photons are 
produced there is also a vacuum fluctuation which is proportional to minus 
the energy if one photon is produced. The coefficient is such that upon 
averaging the energy over the production of zero, one, two ... photons one 
recovers the mean, (zero for (Tyy), see eq. f l2.75p and the thermal flux for 
Tuu (see eq. (^M))- 

To see all of this we first display the pair. This stands in strong analogy 
to Section 12.31 We then consider the non-diagonal matrix elements of T^^ 
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associated to the produced photon at Uq. 

To display the pair we introduce the analog of the Unruh modes (Sec- 
tion 12.31) since they diagonalize the Bogoljubov coefficients. An out mode 
(eq. fl2.72p ) is proportional to 9{—V). By adding to it a piece proportional 
to 9{+V) one can obtain a purely positive frequency mode on T~ . Writing 
only the V part of the modes we have 

= (ak\aVf^''9i-V) + Pk{aVy'/''9i+V)) k > (2.88) 

The U part is given by the reflection condition (12.651) . Analyticity in the 
lower half of the complex V plane (i.e. positive frequency lo in eq. fl2.7UI) ) 
fixes the ratio (see eqs (12.391) et seq.) 

^ = e-"^/" (2.89) 

and the normalization of (fk with the Klein Gordon scalar product fixes 
Q^i ~ /^fc = 1- To obtain a complete orthogonal set of positive frequency 
modes one must include the modes for /c < 

= (P\k\\aV\''/''9{-V) + a\k\{aVy'/''9{+V)) k < (2.90) 

/Aivlkl ^ ' 



From eqs (12.881) and (I2.90p it is seen that the set of "Rindler" type modes 
given by eq. (12.721) and the mode given by 

^^tR ^ Q^J^Y^_^^^YY^k|a ^2.91) 

constitute a complete orthonormal set. It is to be noted that the modes (12.911) 
do not have a f/ part since they don't reflect (when the mirror follows 
eq. (12.691) forever). 

The quantum number k in eqs (12.881) to (I2.9ip is not the usual Minkowski 
energy. Rather it is the eigenvalue of the boost operator iaVdy whereas the 
energy lo is the eigenvalue of idy. (We shall call it therefore "Rindler energy" 
see eq. (I2.47p .) However upon reflection the time dependent Doppler shift 
(I2.63P gives k the meaning of Minkowski energy for out modes. It is then the 
eigenvalue of idij. 



60 



Equations fl2.88|2.90l) therefore are to be read as the Bogoljubov trans- 
formation ( written in terms of eqs 



„, i^outL I Q i-outR* ~1 

~ TLL.\^'^\utR U>0 (2.92) 

^-k = Pk^k + Ctk^k J 

or in terms of operators: 

p. — ^, „outL o outR^ ^ 
O-fc — —Pk<lk + Oik^k ) 

The creation of particles of energy k by the moving mirror gives physical 
content to the modes fl2.9ip as the partners of the created particles since the 
Minkowski vacuum on X~ (|0j„)) can be expressed (see eq. f l2.44p ) as 

]^ 0}^ outL'\ outR-\ 

\0^^)= Hexp^"^ Kut) (2.94) 

V ^ fc>o 

where lO^.^) = \OoutL) ® KutR) with *^|Oo„tL) = and af ^^lOo^t^) = 0. In 
this way one sees that to each produced v^^"*^ particle corresponds a partner 
^pOutR liying on the other side of the horizon {V > 0) with the opposite Rindler 
energy. Upon tracing over the states of different i?-quanta in eq. fl2.94p one 
obtains a thermal density matrix, with temperature T = a/27r, defined on the 
subspace of L-quanta states. This is exactly what was seen upon computing 
the flux on J+ in eqs f l2T6] ^ 12:831) . 

As stated, to each particle created on there corresponds a bump of 
something on the other side {V > 0) which is correlated to it. To understand 
that this correlated bump is really present, consider a charged field and a 
measurement that reveals the production after reflection of a quantum of 
positive charge. Then the correlated bump necessarily has unit negative 
charge. 

To exhibit the correlations between the produced particle and its partner 
configuration we consider 

a packet localized around the line Uq, having mean energy k^. This packet 
issues from a vacuum fluctuation which is propagating in the V direction and 
which (from the reflection condition) is centered around V = —a~^e~°'^° . 
From the perfect symmetry {V —>■ —V) between v'^"*'^ and v'^"*'^, the "part- 
ner fluctuation" is centered around V = +a"^e~"^". This particular configu- 
ration of the field gives one of the contributions to (Oj„|T^j,|Oj„). Our object 
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therefore is to decompose this in-vacuum expectation value into its compo- 
nent parts, these latter constituting a complete set of post-selected photons, 
i.e. out states. Thus we introduce a complete set of localized out-states (wave 
packets), considering the tensor products of all possible states of arbitrary 
numbers of right and left quanta, we write (as in eq. fll.46p ) 

(2.95) 

Let us consider that piece of this expression wherein the post selected states 
are of the form 

(l^'rffc c,ar*^t) \OoutL)\{nlf''}) (2.96) 

i.e. where we have specified that the outL state factor contains one particle 
in the packet 

^= dk clvT''' (2.97) 

JO 

(with the normalization / dk\ck\'^ = 1) but the number or type of i?-quanta 
in K^ifc"*^}) not prescribed. This partial specification is appropriate in the 
present situation wherein only the out L quanta are realized on-shell [63] . 
However we shall soon prove that, due to the correlations in the pure state 
|0j„) the configuration in R is automatically specified as well. 
It is convenient to define the projection operator 

poo roo 

n = W®/ dkcka^'''^0,^tL){OoutL\ dkclar'"- (2.98) 
Jo Jo 

which projects onto states of the form (12.961) since loutR is the identity oper- 
ator in the > region. One may then rewrite the decomposition fl2.95p in 
such manner as to isolate the contribution from 11: 

{Oin\T^u\Oin) = {Oin\UT^u\Oin) + {Oin\{I-Il)T^,\Oin) (2.99) 

We rewrite the first term on the right hand side of eq. (12.991) as 



(o,„|nT^,|o,„) 



(0 



m I n I Oin 



(2.100) 
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{T,u)., = ' r! r^Tr! (2-101) 



SO as to express it as the product of the probabihty of being in the state fl2.96l) 
times the weak value of T^j, in that state (see eq. fll.52p ). In the sequel it is 
this non diagonal matrix element (weak value) 

(o.„|n|oj 

which we shall analyze, since from the above analysis this is the value of T^jy 
which corresponds to the final state n|Oj„) which we now construct explicitly. 

We first check that the specification of the out L quantum uniquely fixes 
the configuration in R to be its partner. To this end we change basis from 
the waves fl^'^^ to a complete set of wave packets labeled by i. The matrix 
of this change of basis is the unitary matrix 7^^. We shall take the wave 
packet i = to be that specified in eq. (12.961) . i.e. 70^ = Cfc. We now rewrite 
the argument of the exponential in eq. fl2.94p so as to isolate the creation 
operators of the wave packets labeled by i: 

f dk ^aZ'^'^al^''^^ = fdk fdk'y2 —aT^'^^lk^ik'alf''^ (2. 102) 

J ak J J ^ ak 

Hence eq. (12.941) becomes 

outL] f ^* I 



Om) = -T= exp i J dk'ck'ttj!^^^'^ J dk cl 




111 Oh 

dk '^ik'ttu, 



outLj 



rffc 7*fc— ar*""^") \Oout) . (2.103) 
ak J 



Note that we have arranged this construction so as to put into evidence the 
combination that creates the observed wave packet (/ci/c'cfc/a^"*'^^). This con- 
struction shows clearly the asymmetry between the particle and its partner 
wave functions induced by the presence of (3k/ ak- This will be crucial in what 
follows. 

Since by construction all the states created by the operators / (i/c'7jfc/a^"*^^ 
(i 7^ 0) are orthogonal to the states involving ^ dk' Ckia'^^^\ the state n|Oj„) 
is easily found to be 

]= [ dk'ck'alf^^^OoutR) [dk cl — 
'Z J J ak 

In this way we see that in the projection 11 of eq. (12.981) onto |0j„) there is 
an implied specification of the partner. This EPR ^9j effect results from the 
global structure of the Heisenberg state |0m). 



n|0,„) = ^ / dk'cyaZ'''^\^outR) I dk c*^ar*''>,„,i) . (2.104) 
\/Z J J ak 
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All is now ready for the evaluation of eq. fl2.10ip . In order to simplify 
the notation we shall calculate rather than T^j,. The latter is 

obtained by taking derivatives with respect to x, x' and then the coincidence 
limit. By expressing the out operators which appear on the right hand side 
of eq. fl2.104p in terms of in operators and writing (f)[x)(j){x') in terms of the 
in basis a straightforward calculation yields 

(O,„|n0(a;)0(x')|O,„) 



(Oi„|n|o,„) 

(Oo.t|0(x)0(x')|O,„) 



+ 



+ ix) ^ (x'] 



S^dk \cunh/auY 

(2.105) 

To derive eq. fl2.105p we use eq. (12.930 to obtain the sequence of equalities: 

1 1 

(Oo«t|apaf00|O,„) = (Oo^tl— a_p — (og + /5af"^)00|Oi„) 



-6{p - g)(Oo«t|00|Oi„) + {Ooutl d^pdg(f)(j)\Oin) 

(2.106) 

Expanding (j) in and making the packet construction indicated in eqs (12.1031) 
and (I2.104P yields the two terms in eq. (12.1050 when the denominator (Oj„|n|Ojn) 
is taken into account. This denominator gives the probability to find the state 
eq. (12.1040 on X"*" . It is given by 



(0,„|n|0,„) = \{0out\0in)n dk\ck\\f3k/ak)^ 

Jo 

1 

= ^ dk \ckW(3u/auf. (2.107) 
Z Jo 

(as for the expression for Pe given after eq. (I1.59P ). 

The first term in eq. (12.1051) is background. It is evaluated by expressing 
i){x) in terms of owt-modes and (pix') in terms of m-modes and one obtains 

{OouMx)(P{x')\Oin) 



{Oin\<P{x)<j){x')\Oin) 

dk^[0l{x)0Ux') + {x ^ x')]{2.108) 

Oik 
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The first term gives the mean value of T^^j calculated in the previous 
subsection. 

The second term is equal to —2 /q°° (3l{\ip'j^\'^ + \ip^\'^)dk . Taking derivatives 
with respect to U and the coincidence limit, gives — /o°° kn{k)dk which is 
the negative Rindler energy density in Rindler vacuum. Indeed, in Rindler 
vacuum (paragraph after eq. (12.431) ). using Rindler coordinates (I2.20p . one 
has 

{ORindler\Tuu\0 Rindler) = -(7r/12)T^ (2.109) 

This can be interpreted as the removal of the thermal distribution of Rindler 
quanta present in Minkowski vacuum. In the present case it corresponds to 
the removal of the energy of all the produced quanta. 
Putting the two contributions together, we thus have 





iTuul 


Oin) {Ooutl 




du(i)\Oin) 


{Ooutl 


Oin) {Ooutl 


lOin) 



= (2.110) 



{Ooutl 


I'^vvl 


lOin) _ {Oout\dv(f> 


dv(p\Oin) VT 1 1 


{Ooutl 


Om) {Ooutl 


|0,„) 12 (27ra)2 {aVy 



since dU/dV = 1/aV see eq. (I2.69p . The absence of outgoing flux in the 
in-out matrix element (I2.110p is natural since it corresponds to a state with 
no out-particle produced. This specification implies that before reflection 
the vacuum fluctuations leading to the production of out-particles be ab- 
sent, hence to the negative Rindler energy (I2.11ip . We note that we have 
neglected to treat properly the singularity at = 0. This shall be analyzed 
subsequently. 

We now consider the second term of eq. (I2.105P (hereafter denoted {T^y)^) 
which contains the contribution of the selected particle iIj eq. (12.970 . Its 
contribution to {Tuu)w is 

''^""^* - ' Wdkwm^) 

To calculate the energy of the particle we shall take a gaussian packet = 
g-«fcc/og-A (fc-fco) /2 ^Y^^QYe the phase factor e~^^^° locates the produced particle 
around U = Uq, its energy being approximatively /cq- Then one verifies by 
saddle point integration that {Tujj)^ is also located around U = Uq and 
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carries also the energy 

As in the electric field, once the selected particle is on mass shell the ip part 
of the weak value behaves classically. 

We now consider the contribution of ip to {Tvv)w Before refiection {U < 
Uq), there is a piece both for < (which upon refiection becomes {Tuu)^) 
and for > (the partner contribution): 





1 {J.^dk' Cl,dyvlf^*) 


J^dk 


\ckm/al) 



2e{+v)- 



S^dk \c,\\Pl/al) 

(2.114) 



After refiection, for U > Uo, only the 6{+V) piece remains since the 6{—V) 
is refiected and gives {Tuu)i)- To grasp the energy content of the particle and 
partner let us first check that each "Rindler" piece (proportional to 0{+V) 
and 6{—V) respectively) carries the same "Rindler" energy (i.e. iaVdy) 
approximatively equal to ko- To this end we introduce v = a~^ln \aV\ and 
analyse the "Rindler" fiux Ty^. One has: J^^dv (Tj,„) = J^°°dV aV{Tyy) 
with the ± corresponding to the 0{±V) pieces respectively. In both cases 
one finds 



^ dv (T,„) = ^ dk k\ck\\Pl/al)j I ^ dk |cfcp(/5,Va^)j ^ 

(2.115) 

For the Q{^—V^ piece this simply follows from eq. (12.1131) since v = U upon 
refiection. For the 9{+V) piece it follows from the orthogonality rules enjoyed 
by the modes. 

Thus {Tyv)w is simply obtained by adding the "Rindler energy" density of 
the 'j/'-Rindler pair to the background, i.e. the first term of eq. fl2.105p which 
we have explained gives the Rindler vacuum see eq. (12.1111) . 

However (Tyy)^ is quite asymmetric and its analysis is more subtle. This 
asymmetry is already present in eq. (I2.104p where the partner wave function 
is given explicitly. The source of the asymmetry in the formalism comes from 
the use of packets which is necessary to exhibit correlations in space-time. 
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Indeed the detection of a particle by a counter requires a description in terms 
of a broad packet for this outgoing particle. Once this is done the partner's 
packet becomes complicated: the convolution of the particle Fourrier com- 
ponents with the energy dependence of the ratio see eq. fl2.104p . One 
loses therefore the simplicity of the Bogoljubov transformation between Un- 
ruh and Rindler modes (eq. fl2.92p ). Hence the two terms on the right hand 
side of eq. fl2.114p are not symmetric with respect to = 0. 

Continuing with the configuration before reflection we emphasize that 



for U < Uq. To see this, recall that J^^dV TyvlOin) = on X~ in virtue 
of normal ordering (i.e. subtraction of the zero point energy of each mode, 
here the energy carried by each orthogonal packet). Thus J^^dV Twl^in) 
vanishes mode by mode on I~ . 

Equation (12.1171) implies that the energy distribution in eq. (I2.11ip has 
a positive singular contribution at = which exactly compensates the 
negative energy density for V ^ [71], [63]. This singularity will come up 
once again in Sections 12.61 and 1 3.51 where it will play a critical role in ensuring 
the consistency of the theory. 

Let us now decompose J^^dV {Tyy)^ into the contribution from the 
particle jl^dV (Tyy)^ and from the partner j^"^ i^vv)^- 

Starting with the latter one sees that the integral is positive since the 
integrand is (as is seen by inspection of the second term in the right hand 
side of eq. (12.1140 ). One estimates ^^°°dV {Tvv)i) — koe°''^° (since we know 
that Jf^'^dV aV{Tvv)i} — and it may be checked by stationary phase that 
the main contribution comes from Vq — a^^e~°'^°). This relation between 
"Rindler" and Minkowski energy is exactly that given by the Doppler shift 
equation (I2.68p . 

These properties imply that J^^dV {Tvv)%i) is negative and ~ —kQc"-^". 
But we have seen that the "Rindler" energy J^^dV \aV\{Tvv)ii! — is pos- 
itive. How is that possible? The answer is that 0{—V){Tyv)^ is not positive 
definite and contains strong negative oscillations for small values of V. When 
evaluating J^^dV {—aV){Tvv)i> (or upon reflection, dU {Tuu)%p) these 
oscillations give negligible contributions. But in the Minkowski energy the 












(2.116) 



(2.117) 
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oscillations near V = play a dominant role ( because of the weight factor 
\V\~^ in going from Rindler to Minkowski energy ) in such fashion as to 
make the Minkowski energy of particle plus partner vanish. We shall discuss 
similar oscillatory effects in quantitative detail in the next section. 

In conclusion the post selection of an outgoing photon of frequency kg 
entails a distribution of Tyv which exists on both sides of the horizon = 0. 
The partner piece, localized around the line Vq = +a~^e~'^^° propagates for- 
ever and has Minkowski energy equal to the Doppler shifted value given by 
/coe"^°. The particle piece propagates out toU = Uo and is then reflected. Be- 
fore reflection it carries energy equal and opposite to that of the partner. This 
energy has a positive piece equal to /coe"^° centered around V = — a~^e~"'^° 
and a broader oscillating piece which is net negative in such fashion that 
the sum of all contributions vanishes. After reflection this piece has gotten 
converted into Tjju. The energy carried in the reflected wave is almost all 
in the center positive piece (see eq. (12.1121) ). The oscillating negative piece 
is diffuse and negligible. For U > Uq the partner continues alone carrying 
a net positive energy always equal to fcoe'*^". However it is not a quantum 
in the usual sense. It does not manifest itself on the average since the field 
configurations, on the average, are still Minkowski (by causality). To pick 
up the effect of the partner requires an EPR correlation type experiment. 
What is important is that the fluctuations in play when an outgoing parti- 
cle is produced have energies which blow up exponentially and which hugs 
the horizon at exponentially small distances. It is this circumstance which 
constitutes a major hiatus in the more realistic case of black hole collapse. 

Had we post selected the absence of an outgoing photon we would come 
upon an "anti-partner" whose energy is negative being the weighted sum of 
the energies of all the negative energy corresponding to the absence of of 
1, 2, . . . , n, . . . photons. As previously calculated (eq. (12.1111) this is precisely 
the energy of the Rindler vacuum (i.e. the absence of the average thermal 
energy ). The sum of all weak values is of course the net average as given in 
Section fl2X2D . 

Perhaps a more physical way to exhibit the correlations between the emit- 
ted U photons and the V partners is to decelerate the mirror after a while. 
As seen in Fig 12.41 
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Fig. 12.41 The trajectory of a mirror which decelerates after a while. The 
classical trajectories (stationnary phase ) of a pair of particles are indicated 
by wavy lines. It is apparent that the deceleration of the mirror allows the 
partners to be realized on shell. 

ipOutR ^jjj then be reflected at late times, transformed into a real quantum 
and can therefore be observed, in particular in coincidence with v^^"*^. One 
could choose a trajectory which starts from rest V — U = const and come 
back to rest after having followed the trajectory eq. fl2.69p for a while. The 
main point is that after the mirror becomes inertial again the system goes 
back to vacuum plus the radiation that goes out to infinity, always a pure 
state. This was pointed out by Carlitz and Willey in ^21j who suggested 
possible applicability to the black hole problem (this was also mentioned in 
[71] and discussed in [100] ). 
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2.6 The Energy Emitted by the Accelerated 
Detector 



2.6.1 Introduction and Qualitative Description 

Though not directly related to black hole physics we include this section in 
our review, not only because it is of interest in itself, given the rather stormy 
history of the problem, but also because one comes upon concepts that arise 
in the black hole problem as well. 

We shall present in detail the energy emitted by the accelerated detector. 
The first subsection is qualitative and shows how the paradoxes which have 
been raised (see for instance and in historical order [M], [H], (SD], [M], 
in, can be rationalized. In the second subsection the energy emitted 
is computed in perturbation theory and particular attention is paid to how 
it is correlated to the final state of the atom. In the third subsection the 
same decomposition of final states is used to display the vacuum fluctuations 
which induce the transitions giving rise to these final states. 

It has been seen in Section 12.41 that in perturbation theory the leading 
order in Hint corresponds to photon emission in both cases, excitation and 
deexcitation, as it should be since one perturbs Minkowski vacuum. Never- 
theless as Grove pointed out [H], when equilibrium is reached there is no net 
change of the state of the radiation in the quadrant, i?, of the accelerator ex- 
cept for transient effects. We refer to this as Grove's theorem. His argument 
is the following. The accelerator feels the effect of a thermal bath. So first 
consider the inertial two level system in thermal equilibrium. The principal 
ingredients which guarantee the absence of net energy flux to or from the 
atom is the time independence of the Hamiltonian and the stationarity of 
the state of the atom (in the thermodynamic sense), so that each photon 
which is absorbed is re-emitted with the same energy, i.e. energy conser- 
vation results from time translational symmetry. Equilibrium is maintained 
through the implementation of the Einstein conditions. 

This argument is immediately applicable to the accelerator since the fact 
that a = constant implies that his physics is translationally symmetric in 
his proper time (ie. invariance under boosts). Since Minkowski vacuum is 
also an eigenstate of the boost operator, the implication is that the total 
eigenvalue of d/dr is conserved. The dynamical realization is the time aver- 
aged conservation of the energy of the totality of Rindler quanta which are 



70 



absorbed and emitted, in strict analogy to thermal equilibrium. This implies 
no net energy flux. 

The above considerations result in a dilemma since, as we said, both 
excitation and deexcitation of the atom leads to emission of a Minkowski 
photon (see for instance the amplitude eq. (12.481) wherein the first order 
interaction of the atom with the field always results in the creation of a 
photon of frequency uj). The resolution of the dilemma will be shown to lie 
in a global treatment of the radiation field which also takes into account the 
transients due to switching on and off the detector. These transients play a 
role similar to the oscillations encountered near V = in equation (12.1161) 
and (12.1171) which ensure the global vanishing of the energy. 

The above general considerations have been verified in an exact model, 
that of an accelerating harmonic oscillator coupled linearly to the radiation 
field |80], [92], [M]. It is noteworthy that in this work Heisenberg equations 
of motion have been integrated to give a long time steady state solution 
wherein initial conditions become irrelevant. This is what makes the analogy 
to thermal equilibrium possible. Rather than describe the exact oscillator 
system we shall continue with the two level atom in perturbation theory 
(since this is more relevant for the understanding of some corresponding 
problems which come up in the black hole problem). 

A detailed picture of the steady state emerges from the following con- 
sideration. Focus on the ground state of the accelerator which excites by 
absorbing a Rindleron coming in from its left. Then the field configurations 
to its right is depleted of this Rindleron. Since this Rindleron carried positive 
energy, its removal can be described as the emission of negative energy to the 
right. In equilibrium there is also to be considered the process of deexcitation 
corresponding to the emission of positive energy to the right. The Einstein 
relation eq. (12.101) guarantees that the two cancel. 

A nice way to express the physics is to appeal to exact eigenstates of 
the photon field, scattering states. Their energy (Minkowski energy for the 
inertial atom, Rindler energy for the accelerating one) is the same as that of 
the free states and their number is also the same since the scattering matrix 
is unitary. Therefore the average energy of the radiation field of any ( mixed) 
state is unperturbed by the scatterer. 

We now discuss qualitatively the transient behavior in both the Rindler 
and Minkowski representation of the radiation field by calculating the en- 
ergy density emitted in terms of the mean energy momentum tensor, {T^J). 
Take, for example, right movers. The relevant energy density is {Tjjjj) = 
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m/duf). 

By Grove's theorem, one should have (Tu„), the energy density measured 
by a co-accelerator, equal to zero in the absence of transients. (Once more we 
normalize the energy so that the expectation values of Tjju and Tyv vanish 
in Minkowski vacuum, therefore of T„„ and T„„ as well). Now consider the 
transients, with switch on (off) time at (rj) modeled by some function 
/(r) which vanishes outside the interval (Tj,rj) and is equal to 1 inside the 
interval except for a time At during which / passes smoothly from to 
1. We choose Tf — Ti >> Max(AM~^, a~"^) which is the time necessary 
to establish the Golden rule. In addition we take Tf — » At. The 
compensation mechanism based, as it is, on translational invariance in time 
is then no longer operative. Insofar as T^u is concerned this will introduce 
minor effects, but these become dramatic for Tuu owing to the exponential 
character of the Doppler shift measured by the inertial observer. Indeed the 
total Minkowski energy emitted is 

/O r+CO (J^y^ 

dU < Tuu >= du< Tuu > 3777 (2-118) 
-oo J — oo dU 

where du/dU = e"". The first integral in eq. fl2.118p is limited to the domain 
U < because the accelerator lives in R thereby confining right movers (i.e. 
those which contribute to Tuu in the integrand) to the quadrants P and R. 
To pick up the total energy emitted, one must integrate along the surface 
V = Vq with Vq > a~^e"-'^f so that all emitted right movers cross the surface 
(see Fig. (12.51) ). From what has been said the integrand vanishes except at 
the endpoints Ui (= r,) and uj {= tj) whereupon eq 12.1181 integrates to a 
form 

Em = Cfe''^i -de''^^ (2.119) 

where Cj and Cf depend on the exact form of the switching function as it 
turns on and off the interaction. 

Most surprisingly eq. fl2.119p can be approximately written as an integral 
over the naive rates of absorption and of emission of Rindler photons taking 
into account that each transition, be it excitation or deexcitation, is accom- 
panied by the emission of a Doppler shifted Minkowski photon, according to 
the resonance condition uj{t) = AMe°^ (eq; ( 12.491) ) i.e. 

Em = HdT f{T){R+p^+R^p+)AMe''^ 

= {R+p. + i?_p+)AM(e""-f - e""0 + (C".e""-f - C^e^'O (2-120) 
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Here R± are rates of excitation (deexcitation) and p± are probabilities of 
being in excited (ground) states (Einstein's condition eq. ( 12.101) is R+p- = 
R-P+ =constant). The constants C^' and Cj depend on the form of the func- 
tion / near the endpoints. So the integral eq. fl2.120p is of the same form as 
eq. (12.1191) . In the subsequent development we shall see how it is that the 
energy density < Tuu > which appears in eq. (12.1181) can be decomposed 
into a positive steady piece and an interference term. The positive piece has 
an integral of the form eq. (I2.120p and should be interpreted in a similar way. 
The interference term carries no net energy, however it plays an essential role 
in ensuring that causality be respected. In the interval (tj, Tf) far from the 
transients the interference term is negative and exactly compensates the pos- 
itive piece so as to recover Grove's theorem in the region where it should hold, 
this being characterized in good approximation by translational symmetry. 
This detailed method proceeds event by event so as to provide a description 
of the radiation field in all four quadrants and its correlations to the detector. 

[This state of affairs wherein transients have a global content which de- 
pends on the whole history also occurs in the problem of the classical elec- 
tromagnetic field emitted by a uniformly accelerated charge in 3 dimensional 
Minkowski space. Here also one can argue, in a way strongly reminiscent of 
Grove's result, that no radiation should be emitted. To wit, the equivalence 
principle asserts that a uniformly accelerated charge is equivalent to a static 
charge in a uniform gravitational field. In a static situation no energy should 
be emitted (this is confirmed by everyday experience at the earth's surface) 
hence none should be emitted in the accelerated case. On the other hand the 
charge should radiate at a constant rate given by Larmor's formula 



where ^ is the energy radiated per unit proper time of the particle, u"^ is 
the square of the acceleration of the particle. 

This dilemma was resolved by Boulware |14j who showed that both results 
are compatible. The crucial remark is that the equivalence principle is valid 
only "inside" the right quadrant (where the static Rindler coordinates eq. 
(12.21) are valid) and indeed inside R the field is the Coulomb field of the 
charge with no radiation part. However Maxwell's equations imply that 
along the past horizon U = there is a delta like singularity in the field. 
This singularity can be interpreted as the infinitely blue shifted transient 




(2.121) 
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which occurred when the charge was set into acceleration. This singularity 
along the past horizon is unobservable by a coaccelerator, but in the future 
quadrant it propagates and gives rise to a flux of radiation at exactly the rate 
eq. (12.1211) . Hence once more transients in the accelerated frame acquire a 
global content for the Minkowski observer.] 

2.6.2 The energy emitted at 0{g^) 

We first consider the energy emitted during spontaneous excitation of the 
atom. The coupled state of atom and field is 

|^_) = Tc-^/^-'^-'IOm)!-) 

= |0m)|-) - ig rdT'e'^^'-'<P{T')\OM)\+) 

JTi 

JTi JTi 

(2.122) 

where T is the time ordering operator. |+) (|— )) refer to excited (ground) 
state of the atom respectively. The mean energy momentum is (considering 
once more only right movers and therefore Tjju) 

(Tuu)^ = {^-\Tuu\i^-) (2.123) 

The mean Tm; eq. (12.1231) is dissected into its constituent parts by consider- 
ing the final state of the atom ^94j, ^J, [63j (i.e. by making a post selection 
similar to eq. (12.991) where the final state of the radiation field was specified). 
This is realized by inserting into eq. (I2.123P the projectors 11+ (n_ = 1 — 11+) 
onto the excited (ground) state of the atom: 

{Tuu)- = {^p.\U+Tuu\i^-) + i^-lU^Tuuli^-) (2.124) 

We have 

JTi JTi 

(2.125) 

Equation (12.1251) is the energy emitted due to excitation of the atom as it 
would be calculated in lowest order perturbation theory ( 0{g) in the wave 
function, 0{g'^) in the energy). But, to 0{g^), one requires the correction to 
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the wave function as well, thereby leading to interference terms. These are 
in the n_ term of eq. 12.1241 corresponding to emission and reabsorption in 
the wave function. This term is 



{ij.\U_Tuu\i^-) = -2/Re(0 



M 



{T2) / dne 

Jt2 



It is convenient to reexpress eq. (12.1261) as 

{ij.\U_Tuu\i^-) = -g'Re[D{U) + J^{U)] 



(2.126) 
(2.127) 



where 
D{U) 



rfr2e-*'^^^^-0(r2) / ' dne'''''-^ct>{n)Tuu\^M) 



7 rf 
dT2 / dTie(T2 - ri)e" 



*^*'"¥(r2)e*^*^^^0(ri),Tf 



uu 



\0m) 



(2.128) 



Complications from time ordering are no longer present in D{U). The term 
T{U) plays no important role in the physics. Indeed a detailed analysis shows 
that it enjoys the following properties: 

1) It is smaller than D{U) by a factor l/AM(r/ — r^) except in the 
transitory regime where it is comparable to D. 

2) It vanishes in the non causal domain U > 0. 

3) The integral / dUT{U) vanishes (ie. it does not contribute to Em = 
JdUTuu). 

4) The integral / du{dU / duYJ^{U{u)) vanishes also (ie. It does not con- 
tribute to the Rindler energy = J duTuu)- 

All these properties are also valid when one introduces a switch on and 
off function /(r) as in equation eq. fl2.142p below. Hence from now on we 
shall drop the term T{U). 

To execute the calculation, and so check out Grove's theorem at this order, 
we need to evaluate {il)^\U.^Tuu\il)_) and D{U) in the limit Tf — ti — 00. 
Begin with eq. (12.1251) and expand 0(r2) and (t>iji) iii Unruh modes (eq. 
(I2.43P ). These are the most convenient because their vacuum is Minkowski. 
Then 0(r2) (and 0(ti)) contain the creation of an Unruh mode in R. Since 
AM > 0, the part of (j){T2) which contributes to the integral (i.e. resonates 
with the factor e*^^''^'^^^ jg ^^j-^at which corresponds to the annihilation of a 
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Rindler mode in R, hence carrying a factor Px in (fx (eq. fl2.43l) ). The result 
is thus (compare with eq. fl2.55p ) 

J c?re^^^^Xr)|OM) = /9AMv/^/AMaL,|0M). (2.129) 

This Unruh creation operator must be contracted with the corresponding 
Unruh annihilation operator appearing in to give 



^ {Oi-U)Pi,, + 9{U)PlMalM) (2.130) 



2 {aUy 

The calculation of eq. (12.1261) proceeds along similar lines and yields 



^2 I 

^^jjy pAMOiAM (2.131) 

The important minus sign is the same as comes up in the verification of 
unitarity wherein to 0{g^) the interference term cancels the direct term in 

1 = {ip\i)) = {ipo + iSiplipo - iSip) 

= {i/jom + {s^pm +2im {i/jom = {h^^o) (2.132) 

This minus sign is in fact essential to satisfy causality which requires that 
all expectation values of vanish for U > 0. The photons emitted cannot 
affect the region U > 0. Indeed there is a rigorous theorem stating that 
{ip-\Tuu\ip_) = for U > which we now prove 

{Tuu{U>0))- = (-|(OAf|Te*/'^"^'"'T[/t/(f/ > 0)Te-'/'^"^^"*|OM)|-) 

= {-\{Om\Tuu{U > 0)Te*/'^"^»"*Te-^/'^"^^"*|OM)|-) 

= (-|(OM|Tc/t/(f/>0)|OM)|-) 

= (2.133) 

where we have used the commutativity of 0(r) appearing in Hint (i-e. on the 
accelerating trajectory) with 4>{U) for f/ > appearing in Tm;. The same 
proof applies immediately if the initial state is \ip+) = |0m)|+)- [This proof 
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may be also generalized to massive fields and interacting fields for Tuu{t^x) 
with (t, x) in L since then it is space like separated from 0(r) but in general 
we will have {Tuuit,x)) ^ for (t,x) in F.] 

The sum of eqs. (12.1301) and 12.1311 gives a negative result to the average 
energy density 

(Tuu)- = -jf^lMj^ (2-134) 

in accord with the expectation that absorption of a Rindleron diminishes the 
energy density. More precisely the total reduction of Rindler energy in the 
interval (Tj,r/) at 0{g'^) due to absorption of right movers is 

{Tuu)-du = J {Tuu)^{-^fdu = -—plMi^f - Ti) 

= -^AMiTf-T,) (2.135) 

where we have used eq. (12. lip . Dividing eq. (I2.135P by the probability of 
excitation P+ = i?+(r/ — r^) yields the energy emitted in U modes if the atom 
is found excited 

/ {T^u)-du = -AM/2 (2.136) 

Thus one obtains a steady absorption of energy exactly as in the usual 
golden rule. The factor 1/2 in eq. (12.1350 arises because we have taken 
into account right movers only. We have set the integral du = Tf — t^. 
This very reasonable result can be obtained rigorously by going to the limit 
Tj — Tj — >■ cxD in a more controlled manner (see [63]). 

The same procedure as that discussed in the paragraph preceding equa- 
tion eq. (12.1301) applies to the energy emitted if the initial state is One 
finds 

g\, e{-u ) 
{auy 



= '-o^lu-hr^ (2-137) 



whence 

R 



J {^+\T^^\^+)du = —AM{Tf - n) (2.138) 

We can now check out Grove's theorem for the equilibrium situation. In 
the inside region (tj < r < tj) one has at thermal equilibrium 



= = e-'^^*^ = PIm/c^Im (2.139) 



77 



where p_ is the probabihty to find the atom in the ground (excited) 
state. This guarantees no net fiux: 



V+{Tuu)++V-{Tuu) 







(2.140) 



It is important to remark that the separate contributions to (Tuu)- given 
by eq. fl2.130p and eq. (12.1311) are each non causal. Each has 0{+U) contribu- 
tions which cancel in the sum. These non causal contributions occur because 
in the calculation of eq. (12.1251) and eq. (12.1261) we have introduced the pro- 
jectors n+ and n_ (the 11+ contribution was already anticipated in the simple 
saddle point approximation (Section 12.41) . They encode the correlations in 
Minkowski vacuum between L and R Rindlerons (recall eq. (12.441) ). Thus in 
the description wherein an atom excites (deexcites) by means of Rindleron 
absorption (emission) its transitions correlate to the presence (absence) of 
the corresponding space like separated Rindler quantum in the other quad- 
rant {U > 0). This is summarized in Fig. (12.51) which is presented at the end 
of this section. 

Instrumental in the explicit realization of Grove's theorem is the negativ- 
ity of {Tuu)~ and hence of {Tjju)_ in the region where translational symmetry 
is valiqj (recall we have made the computation in the limit rj — — oo). 
But one has dU {Tjjjj)- > since to order g'^ the contribution to it from 
eq. (12.1311) must in fact vanish on the basis of a rigorous theorem 



After this manuscript was completed, further research on this subject was done. It 
has now been proven [72] that for a detector of finite mass M (but nevertheless with 
M/a >> 1) the interference term becomes negligible after a few transitions, thereby 
reinstating the validity of the naive Born approximation for the energy emitted. This 
occurs because, when recoil is taken into account, the detector shifts its orbit from p = 

to an orbit characterized by a new horizon (ie. the center of the hyperbola which 
describes the detector's trajectory shifts). Thus one loses the translational invariance in t 
(boost invariance) and Grove's theorem is no longer applicable. More formally, the term 
eq. (j2.13ip exists because the atom which has emitted a photon and then reabsorbed it 
interferes with the atom which has not made any transitions. When the atom recoils 
these two amplitudes no longer interfere destructively. Recoil induces decoherence. This 
decoherence occurs after a logarithmically short time r ~ aln(A//AA/). This very short 
time is a manifestation of the exponentially large frequencies which resonate with the 
accelerated atom at early times (see eq. (|2.49[) ). The same exponential rise of frequencies 
in the Hawking radiation is a source of anguish when the gravitational back-reaction of 
these frequencies is considered. More on this in Section [3.71 




(2.141) 
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This is because J_^dU Tuu is the total energy operator carried by right 
movers; hence it annihilates the vacuum. (Recall we are normalizing the 
vacuum energy to zero). Furthermore the integral J^^dU {■ip-\Il+Tuu\'4'-) 
is strictly positive (it is the expectation value of a positive definite operator: 
the hamiltonian, see eq. fl2.125p ). 

Thus the sum of the two is positive, as stated. Therefore there is positive 
energy density which does not appear in eq. (I2.134p and which can only come 
from transient behavior at the end points. 

In order to study the details of the energy distribution one must intro- 
duce explicitly a switching function, for instance by taking the interaction 
Hamiltonian to be of the form 

H.Ur) = /(r)e-*^*^V_0(r) + h.c. (2.142) 

where /(r) controls the switching on and off of the interaction. In order 
for /(r) not to induce spurious switch on and off effects it should have a 
sufficiently long plateau Ar that the golden rule can establish itself (Ar >> 
AM~^ and Ar >> a'^) (see discussion after eq. (12.81) ). Furtherore in ^3] it 
is shown that if /(r) obeys the condition 

J cire'^l"l/(r) < oo ^ J dtf^rit)) < oo (2.143) 

then the energy fiuxes are regular and no singularities appear on the horizons. 

When /(r) does not obey eq. (I2.143P the transients do not appear in 
finite Rindler time. Rather they are found on the horizons f/ = and 1^ = 
where they give rise to singular energy fluxes. This situation is analyzed by 
regulating the Bogoljubov coefficients eq. fl2.28p which amounts to replacing 
U hj U — ie as in eq. fl2.4ip . When this is done one finds that the two terms 
eq. and eq. fl^TTM]) take the form 

(^_|n+T^H^_) = ^/3iM«iM(f/-^e)-^^^-^/'^-i([/ + .e)^^*^/'^-i 



^2f^lMjj^, (c^IaAU) + PlMOi-U)) (2.144) 



and 



{^.\U.Tuu\i^-) = -;f^/5iM«L/Re[(f/ + Z6)-^^-^/'^-l(f/ + Z6)^^*^/"-l 

2a^ L 



9 q2 ^2 



2a? 



{U + ief 



(2.145) 
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Thus for \U\ > e, eq. (12.1341) remains valid whereas for \U\ < e all terms 
are positive. The total energy emitted is obtained by integrating over U (by 
contour integration). The contribution from eq. fl2.145p vanishes as required 
whereas eq. (12.1441) gives apositive contribution equal to {g'^/2a^){n/e)P\^ja\j^ 
ie. total positive energy is radiated, one part of which is concentrated along 
the horizon and is positive, and the rest is in the detectors quadrant. In this 
way both causality and the positivity of energy of Minkowski excitations are 
respected. 

A similar analysis is possible for (ip^lTuulip^). In the equilibrium situa- 
tion the terms from {ip_\Tuu\ip_) and (ip+lTuulip^) combine in such a way 
that they give zero everywhere except for a singular positive energy flux on 
the horizon —e<U<e. 

This completes the formal proof of the qualitative discussion presented 
in the preceding subsection. We thus have shown that as in the accelerated 
charge problem considered by Boulware this singular flux can be interpreted 
as infinitely blueshifted transients which occurred when the atom was set 
into acceleration at r = a~^lnae. When eq. (I2.143P is satisfied complicated 
expressions arise wherein the role of e is played by Ar~^. As these functions 
are not very interesting to display here, we reserve a more interesting spot 
in Section [2.6.31 to put them on exhibition (see Fig. (12.51) ). 

It is now clear how eq. (I2.120p makes sense. The total contribution to 
I~^dU {4'^\Ti/u\ip_) is from the 11+ contribution only and hence can be 
expressed as coming almost entirely from the wrong quadrant. Thus for 

<< 1 one may write the n_|_ contribution as an integral over ul (see 
eq. (12.240 ) thereby obtaining the absorption part of eq. (I2.120p . The emission 
part (ie. the contribution of |^/'+)) is obtained in similar manner to yield 
the sum eq. (I2.120p . Causality is verified when one analyzes the complete 
distribution of energy, taking into account the n_ interference term. The 
mean energy density radiated is then found only in the transients. 

2.6.3 The Vacuum Fluctuations Correlated to the Ex- 
citations of the Atom 

In the previous subsection, we have analyzed the mean energy radiated by the 
atom and we saw how it can be decomposed into two contributions. These 
correspond to the energy emitted by the atom when it is found excited or 
not excited at t = cx). We now address the question: what configurations 
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of energy-momentum were present in Minkowski vacuum which give rise to 
spontaneous excitation of the two level atom? They certainly carry zero total 
energy but locally should have a positive Rindler energy density to excite the 
atom. 

We shall use the weak- value formalism since it is shown in Appendix [C] 
that it provides a framework to investigate on X+ or on the nature of the 
field configurations correlated to the excitation of the atom. 

Before the atom interacts with the radiation, we have that in the mean 

(-|(0Af|Tf/[7(J-)|0M)|-) = (2.146) 

where Tuu{T~) is the Heisenberg operator, Tuu, evaluated on the surface 
X~, ie. on a surface which is temporally situated before the time Tj when the 
atom begins to interact with the field. 

In this mean appear two cancelling contributions according to the class 
of final states considered: contributions for which the atom excites in the 
period (tj, t/) where both r, > t and Tf > t and that for which the atom does 
not excite in this same time interval. Thus we must project |0a/)|— ), the 
Schrodinger state at times r < Tj, into the various outcomes that are realized 
at later times. We carry this out by inserting at t = +00 the projectors 11+ 
and n_ (see eq. (EH^H)) 

= {-\{0m\Tuu{I-)\0m)\-) = (-|(OM|Te^^^'''''''"'Te-'^^''''''"'Tt;,;|OM)|-) 
= (-K0M|Te'^^7"^^-*(n+ + n_)Te"'-^^7''^^"-Tt;c;|0M)|-) (2.147) 

Let us examine the contribution 

(-1 (Oa/ |Te' ''"''"^'n+Te"' ""^""'""TuuIOm) \ -) (2. 148) 

where 11+ is the Heisenberg operator which projects onto the excited state at 
T = Tf. This contribution is anticipated to be positive since it is the contri- 
bution to the mean {Om\Tuu{I~)\Om) that results in excitation. The correct 
normalization which appears in the weak value formalism (see eqs. (11.521) 
and fl2:T00D ) is to rewrite eq. 02.1481) as P+{TuuiI~))w+ where P+ is the 
probability for excitation in the interval {Ti,Tf): 

p+ = (^-|n+|v^_) 

(T^^(J-))^+ = -L(^_|n+T^^(J-)|^_) (2.149) 
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In this rewriting {TuuiX'))^^ is what has been identified in Section [T^ and 
2.51 and Appendix O with the energy of the field configuration which gives 
rise to excitation of the atom. 
To order g"^ we find 



{Tuu{I-))r 



{-\{0m\ {l + ^9I;/ e^^^-0(r)cr+) n+ {h.c.)Tuu\OM)\-) 



{-\{0m\ (1 + zgj;/dT e^AA/r0(^)^^j n+ (h.c.) |0m)|-) 



(2.150) 



As in eq. (12.1251) the g"^ term in the wave function does not contribute to 
matrix elements when 11+ is inserted. 

It is instructive once more to consider the resonant piece (ie. t/ — —>■ oo) 
of {Tuu(I~))w+ even though we know from the previous subsection that 
the transients are not correctly described in this approximation. Thus we 
replace the double integral in Tuu by a simple integral and take all Rindler 
frequencies to be equal to the resonant frequency AM and obtain 

Pam 

AM 1 , , , 

Where we have used eq. ( I2.38p . Notice that {Tuu{T~))w+ is non vanishing 
both for f/ > and U < 0. It can be interpreted in R by appealing to the 
isomorphism with the thermal bath. Since we have post-selected that the 
atom will get excited necessarily there was a Rindler quanta in R which will 
excite the atom. The weak value therefore contains the Rindler energy of 
this particle. The factor a\j^.^ = n{AM) + 1 rather than 1 takes into account 
that in a thermal bath their may be more than one quantum. The two level 
atom is sensitive to this since it responds to the mean number of quanta. 

In L, the partner of the Rindleron in R appears with the same Rindler 
energy AM because Minkowski vacuum is filled with correlated Rindlerons 
in R and L (eq. (I2.44p ) of zero total Rindler energy eq. (12.471) . But since 
we are in Minkowski vacuum and the interaction has not yet occurred, we 
necessarily have 

/+00 
dU{Tuu{I-))n.+ = (2.152) 
-oo 
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because the operator dUTjjij annihilates Minkowski vacuum. Hence 
{Tuu{X~))w^ carries no integrated energy as behooves a vacuum fluctuation 
(see eqs. (EHUD and flMTTj) ). 

Upon integrating eq. (12.1511) a contradiction seems to arise as in the 
previous section. The transients have been incorrectly taken into account. A 
correct handling of the ultraviolet Minkowski frequencies (eq. (12.411) ) solves 
this problem and in the limit — — > oo, the weak value is the distribution 

{Tuuil-)).^ = ^^I(^7^«AM (2.153) 

with a regularized energy content on the horizon U = which restores the 
property eq. (12.1521) (verified by closing the contour in the complex plane 
whereupon the double pole at [/ = —ie has zero residue). 

We now consider the model eq. (12.1421) so as to have smooth Tuu^s (func- 
tions rather than distributions). The manner in which eq. ( 12.152^ is realized 
is a subtle interplay of several effects which we now summarize. 

We first obtain that for > the weak value is real and positive since 

/T ((T- TT ^ nn - (Oa/I ^drHUr) Jdr' H,^,{t')Tuu\Om) 
{Iuu[[l , > U);.+ - ^Q^i jdrH,^^{r) Jdr'HUr'MOM) 

{Om\ JdTHint{r)Tuu Jdr' H,ntir')\OM\ 
(OmI SdrH^r) jdT'Hi^t{T')\QM) ^ ' ' 

which is manifestly real. We have used the causality development of eq. (12.1331) 
to make the necessary commutation. Moreover eq. (12.1541) is positive since 
it is the expectation value of Tjju in a one particle state. 

Since the integral over all U vanishes, the integral over f/ < must yield 
a negative result which exactly compensates the integral over f/ > 0. On the 
other hand the Rindler energy J^^du T^ui^'] f/ < 0) is positive (although 
the integrand contains small end effects) since it describes the energy of the 
Rindleron which shall be absorbed by the atom. The net negative Minkowski 
energy on the right {{U < 0) once more finds its origin in the transients which 
are small in the Rindler description but large in the Minkowski description 
due to the enormous Doppler shift du/dU = e""^. Furthermore for U < 0, 
the weak value is not real. Nevertheless the complexity appears only in the 
transients (in eq. (I2.150p at U = 0). This imaginary part also integrates 
to zero by virtue of eq. (I2.152p . In Section 13.51 we shall dwell somewhat 
further on the physical content of the imaginary part of (T^j^(X~))^+ (see 

eq. mm)- 
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The reader may check that the energy fluctuations correlated to the tran- 
sitions of the atom posess all the properties of the fluctuations correlated to 
the production of a specific photon by a mirror as studied in the previous 
section. However whereas in the mirror the post selection was an abstract 
construction, we have here shown that it is realized operationally by coupling 
an external system to the radiation. 

Putting together the analysis of this subsection and the preceding one we 
have obtained a description in terms of energy density of the field configura- 
tions in Minkowski vacuum which make an accelerated two level atom excite 
(eq. fl2.150p ) and the resulting field configuration after excitation (eq. fl2.130p ). 
This is depicted shematicaly in Fig. (12. 5b ) for an atom which begins in the 
ground state around time Tj and excites in the interval between and Tf. If 
no excitation occurs one has the complement of eq. fl2.150p and fl2.130p . This 
is depicted in Fig. (12. 5b ). The mean energy momentum is recovered when 
summing over excitation and no excitation and is depicted in Fig. (12.5b ). 

For simplicity of drawing it is the Rindler energy density T^u which is 
drawn in every case. This is because the rapid fluctuations of the Minkowski 
energy in the vicinity of the horizon makes it virtualy impossible to represent 
in a comprehensible diagram. In all three figures the trajectory of the atom 
is drawn as a hyperbola between the points and Tf. Energy fluxes for 
U-modes are represented by dotted lines. These flow from the past energy 
configuration (drawn in the SO corner of the diagram) to the future energy 
configuration (drawn in the NE corner). 
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Fig. 12.51 Birds 'eve view of the energy of the field configurations both in the 
past and in the future of a two level atom which begins in the ground state 
around time Ti and which either does (Fig. a.) or does not (Fig. h) excite in 
the interval between Ti and tj. Fig. c represents the mean effect. 



In Fig (a) the atom has become excited as is witnessed by the absorption 
of the positive Rindler peak in the middle of the configuration for U < 0. In 
consequence this peak has disappeared in the future and there only remains 
the partner in the region U > 0. (In fact the energy in the future for U < 
is not quite zero but is order e~^^^'^ taken to be e~^^*^ ~ 10^^ in this case 
so that this is in fact too small to draw). 

Fig (b) is the complementary case in which the atom has not become excited. 
In order to make the drawing visible the scale of the axes has been changed 
by a factor ~ 10~^. 
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Fig. 2.5 b 



The mean effect (fig (c)) therefore having mean energy zero in the past gives 
rise to a negative Rindler energy in the future (illustrated by the central 
negative dip in the future). As emphasized in the text there are transient 
effects which carry positive energy and which contrive to render the mean 
emitted Minkowski energy positive. It is further to be noted that causality is 
respected in that the partner contributions {U > 0) have cancelled between 
fig (a) and (b). 
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Fig. 2.5 c 
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Chapitre 3 

Black Hole Evaporation 



3.1 Kinematics 

The reader is referred (for instance) to references [IQ] , [1] , [2] , [21] , [60] for an 
introduction to classical black hole physics. For the conceptual issues raised 
by quantum mechanics in the presence of horizons it suffices to work with 
the Schwarzschild black hole. Not that the Kerr hole does not give rise to 
interesting effects, but its complications appear out of context in the present 
review. By Schwarzschild black hole we include, and indeed mostly discuss, 
the incipient black hole wherein the star's matter has fallen into a region 
which is asymptotically close to its Schwarzschild radius, its future event 
horizon at rs = 2M. Throughout we shall take the Planck mass equal to 
unity. Thus r is measured in Planckian distances and M in Planck masses. 
For a star the size of the sun where M = 1.1 10^'' proton masses = 0.9 10^^ 
Planck masses, we have rs = 1.8 10'^^ Planck distances. 
Outside the star the metric is Schwarzschild 



(1 



2M 



)de + (1 



2M 



(3.1) 



r 



r 



(1 



2M 



)i-dt^ + dr*^)+r^dn^ 



r 



(1 



2M 



)du dv + r'^dVt^ 



(3.2) 



r 



where r* is the tortoise coordinate 



r - 2M 



and 



dn'^ = d9^ + sin^ 9 dip'^ 



(3.4) 



Thus radial light rays follow u or f = constant where 

l]=tTr*. (3.5) 

In addition to the Schwarzschild coordinates displayed in eqs (13.11 13.21) 
we shall refer to both Kruskal [57], [ID], and advanced Eddington Finkelstein 
[32] , [31] , [lQ| coordinates. In the first Schwarzschild quadrant (R), in which 
r > 2M, Kruskal coordinates are defined by 

V = AMe-l^^ . ^ 

U = -4Me-"/^^ ^ ^ 

The relation between the Kruskal U, V and the Schwarzschild u, v is thus 
identical to the relation, eq. (12.201) . between Minkowski and Rindler light-like 
coordinates. This isomorphism will play a crucial role in the understanding 
of the various properties of the Hawking radiation. In Kruskal coordinates 
the metric reads 

r 

where r is given imphcitly by UV = {AMY{1~ r /2M)e''/^^^ , and radial light 
cones are surfaces of constant U or V . 

Were the complete space given by the analytic extension of the geometry 
whose metric is eq. (13. 7p there would be four quadrants (see Fig. (13. ip ) sep- 
arated by horizons f/ = and V = Q the first of which, R, is coordinatized 
by the Schwarzschild coordinates m, v of eq. (13. 5p . The other three are co- 
ordinatized by Schwarzschild local coordinates in a manner analogous to the 
coordinatization Minkowski space into four Rindler quadrants. 
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Fig. 13.11 Penrose diagram of the maximal analytic extension of 
Schwarzschild space. The horizons, the singularities as well as r = const 
and t = const lines have been represented. 



For the collapsing black hole the outside space is confined to two of these 
quadrants R and F (see Fig. fl3.2p ). 

Whereas Kruskal coordinates can be used to describe both of them, the 
equations which relate t and r to U and V through eq. (13.31) to eq. (13.61) 
are good only for the quadrant R. In the F quadrant, one may introduce up 
given hj U = 4Me"^/^^. The relation between V and v is still given by 
eq. (13.61) since v is finite on the future horizon U = 0. Then, with t,r given 
in terms of Mi?, f as in eq. (13. 5p the metric is once more eq. (13.11) . Note that 
in R, t and r are time-like and space-like variables respectively whereas in F, 
where r — 2M < 0, t is space-like and r time like. Finally, we introduce the 
advanced Eddington-Finkelstein set v, r that covers both R and F and will 
be found convenient when we study the back reaction. In regions exterior to 
the star one has 

ds'^ = -{l-—)dv^ + 2dvdr + r'^dn'^ (3.8) 
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Fig. 13.21 aThe exterior of the star (in white) is a recopy of the relevant 
part of the complete Penrose diagram Fig. (3. 1 ) whereas the grey part 
representing the inside of the star is obtained by extending ingoing and 
outgoing light cones through the two regions. The locus of the apexes of 
these cones being the line r = 0. As in all Penrose diagrams radial light 
rays are represented by straight lines at 45° degrees In te to the vertical, b 
The same Penrose diagram as in a redrawn by availing ourselves of further 
reparametrisation freedom of radial null rays inherent to the construction in 
such manner that the line r = is drawn as a straight line. This is the 
Penrose diagram which is usualy found in standard texts. The points 
labeled Pj, Pc, Pq, I and O refer to the discussion of Hawking radiation in 
Appendix D 

In the relevant quadrants the domain of variation of the variables are: 
• Schwarzschild in R 

— OO < U < +00 , — OO < V < +00 

also - OO < t < +00 , 2M < r < oo (or - oo < r* < +oo) (3.9) 
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• Kruskal in R 



- oo < U < , 0<V<oo 



(3.10) 



• Kruskal in F 



0<f/<oo, 0<V<ooandt/V<l 



(3.11) 



• Eddington Finkelstein in R and F 



< r < oo , — oo < f < oo 



(3.12) 



So much for the geometry of Schwarzschild space, i.e. the exterior of the 
star. We must now describe the interior of the star. A convenient ideahzation, 
first used by Unruh [9T]|, is the model of a collapsing shell wherein the interior 
is empty, hence described by flat space. We shall pursue this case explicitly 
and show how, in Appendix [D] the essential result of the analysis emerges 
from the general case. 

Inside the space is Minkowski so that 



where V = T =p r. We have set the spatial coordinate r to be the same in 
both regions so that the area of spheres (SO (3) orbits) is Anr"^ everywhere 
and for all times. 

In Fig. (13. 2p we have drawn the Penrose diagram which shows the salient 
features of the collapsing geometry. In this diagram we do not take into 
account the loss of mass due to the Hawking radiation, hence it is only useful 
to describe the early stages of evaporation. The heavy line traces out the 
trajectory of the surface of the star Rst{v). We parametrize the trajectory 
in terms of the Eddington Finkelstein time v rather than the Schwarzschild 
time t, so as to cover its whole history. At v = vq the shell reaches its horizon 
H, the light like surface given by r = 2M. One may extend this light like 
surface H into the interior of the star as indicated. 

At still later times the shell reaches the space-like singularity r = 0. This 
is a 3-surface not to be confused with the 1-dimensional time-like line r = 
which is the axis of SO (3) symmetry of the whole collapsing geometry. This 
is brought out in Fig. (13.31) . This line of symmetry is the locus of the vertices 
of light cones which trace out the paths of spherical waves in the journey 



ds^ = -dT^ + dr^ + r^dn^ = -dUdV + r'^d^f 



(3.13) 
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from to X+ (light like past and future infinity). In Fig. (13. 2p these cones 
are represented by lines reflected off r = 0. This picture is blown up into 3 
dimensions in Fig. (13. 3p wherein spheres are represented by circles. 




I axis r=Q 



Fig. 13.31 The Penrose diagram Fig. (3. 2 b) for a collapsing star blown up 
to 2 + 1 dimensions. We have represented the horizon, the singularity, the 
surface of the star, a typical light cone. The dotted line represents a 
r = const line. 

We have also represented the horizon light cone, H, its extension into the star 
and its backward history as a light cone that terminates on X" . A spherical 
flash of light on this cone is the last flash that can be emitted from X~ so 
as to arrive at X"*". It will be shown that Hawking radiation is concerned 
with the combined Doppler and gravitational red shifts which spherical light 
cones, traveling just before H, experience in their trip from X^ to X+. 

Therefore, in order to compute those combined effects, we need to match 
the internal coordinate system U, V to the external Schwarzschild set u, v. 
We now turn to this task. 

Since u =const. defines a future half light cone, it can be extended back 
in time and enter the star with apex at r = 0; similarly for v. Thus the light 
cones labeled outside the star by the Schwarzschild coordinates u, v can be 
used to coordinatize points inside the star. The same future half light cone 
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is described either by a constant value of or u (and similarly for the past 
half null cones in terms of V and v). In consequence W is a function of u only 
and V a function of v only. 

The radial coordinate r has the same meaning inside and outside the star. 
Thus on the shell the displacement in r is the same in both systems: 

2dr^dV-dU^{l-'^){dv-du) (3.14) 

Rst 

where differentials are along the trajectory. A further condition is that in- 
tervals of proper time on the shell be the same in both systems. Thus 

dVdU = (1 - '^)dvdu (3.15) 
Rst 

These equations together with a trajectory for the star's surface Rgtiv) 
are sufficient to solve for U{u) and V{v). We rewrite them as 

dV dU du _ 2M du. , /o ie\ 

dv du ^ Rst 

When the shell is far from its horizon {Rgt » 2M) U{u) and V{v) are 
slowly varying functions the exact form of which is irrelevant for our purpose. 
On the other hand when the shell approaches its horizon asymptotically 
these functions acquire a universal behavior characterized by M only. As 
it is precisely this domain of the variables which is relevant to describe the 
steady state Hawking evaporation we shall consider this case. 

To characterize the shell's trajectory near the horizon is very simple. 
Suppose Rst crosses the horizon at some Eddington Finkelstein advanced 
time f with finite acceleration (one verifies that this is indeed the case when 
the shell follows a geodesic) . Then to first order we have 

2M - Rst{v) = k{v - Vq) (A;>0) (3.18) 

Here v{= t-\-r*) is equal to t-\-R*^ on the surface so this same equation gives 

Rst{t) 

RJt) - 2M ^ , , , , 

- = Rltit) +t-vo (3.19) 
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Near the horizon R*^ ~ 2M + 2M\n{Rst/2M — 1) so we can solve iteratively 
to give 

Rstit) -2M = 2Me(^°-2^'^-*)/2A^ + O(e-*/*0 = Ae~'/^'' + O(e-*/*0 (3-20) 

where A is a positive constant. Moreover to the same approximation we may 
replace t by + fo)/2 so as to yield 

Rst{u) -2M = A'e""/^^ + 0(e-"/'*^) (3.21) 

where A' is another positive constant. The relation between u and v on the 
surface is thus 

v-vo = -^e-"/^^ + 0(e-"/2^) (3.22) 

rh 

whence asymptotically 

dv. _ A' _ Rstju) - 2M 

J^* " AMk ~ AMk ^ ' 

Inserting eq. (I3.23P into eq. (13.161) and eq. (13.171) one sees that the first 
term on the r.h.s. of eq. ( I3.16P is negligible with respect to du/dv. Further 
we have in the limit v —>■ vq, dV/dv\yg = A where A is some positive constant. 
It then follows that 

^ = Be-^/'^' + 0(e-^/2M) ^ ^ , 

du AM ^ ' ^ ' 

where B is another positive constant. These constants will be discussed in 
Appendix [Dl see also ^33j . Their value is irrelevant for the calculation of the 
asymptotic Hawking radiation. The upshot is that near the horizon 

V-4M = A(f-fo) 

U = 5(-4Me-"/'*0 +0(e~"/2Af^ ^325) 

Very important is the fact that U tends exponentially fast to Kruskal U 
(defined in eq. (13. 6p ) and that V does not but rather is lineary related to the 
Schwarzschild v. Note also that 

^ = (5M')(i?.i-2M) (3.26) 
du 
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i.e. is proportional to goo, rather than the usual ^/goo■ A direct calculation 
reveals that the factor Rgt — 2M is a composite of the static red shift and a 
Doppler shift due to the retreating surface (see Appendix [D|) . 

In eqs fl3.24[ I3.25P we have set V = 4M and W = at the point where 
the star crosses the horizon H. The light cone H which shall generate the 
horizon is given by W = after it reaches r = and therefore by V = 
before (since 2r = V — W and r = at the apex of H). 

Fig. (13. 4p gives a sketch of the exterior R and interior I of the star in 
the orthogonal coordinate grid u, v which have their usual meaning in R as 
functions of t, r eq. (13. 5p and where in I we have 

2r = V -U = Xv + —e-''/^^ + constant 

AM 

2T = V + U = \v - —e~''l'^^' + constant' (3.27) 




Fig. 13.41 The geometry of a collapsing star in u, v coordinates 



One should compare this drawing with the Penrose diagram Fig. (13. ID and 
compare the trajectory of light rays in the two drawings. One should also 
compare the shell's trajectory eq. ( 13.22^ with the mirror trajectory eq. ( 12.69^ 
and remark that the Doppler shifts of eq. (I3.24p and eq. (I2.68P are identical. 
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Therefore, we shall find also a steady thermal fiux (Hawking radiation) in 
the collapsing situation. 



3.2 Hawking Radiation 

As in the preceding chapters, quantization proceeds through the construc- 
tion of modes, here solutions of = 0. In the simple shell model the 
d'Alembertian is a simple operator in the inside and outside domains. Begin 
with the outside (Schwarzschild space) wherein from eq. (13.11) 



nip 



\f9 



(r - 2M) 



rip 

(3.28) 



Here L is the angular momentum operator. Since it commutes with the 
d'Alembertian we pass over immediately to states of fixed I, i.e. modes are 
written as ipi = {^/ATTr)~^ifi(r)Yi^(9,ip) where 



Q-j-2 Qj,*2 



V,(r) 







(3.29) 



with 



Vi{r) 



'1 



2M 2M + l 



-)( 



+ 



Similarly in the inside region one has 



52 92 /(/ + 1)\ 



(3.30) 



(3.31) 



Vi{r) plays the role of a centrifugal barrier, present even for s- waves, where 
it gives a positive potential energy bump at r = 8M/3. It has been shown 
by numerical calculation that about 90% of Hawking radiation is in s- waves 
[SHI, [83]. 

Therefore, in this section, we shall restrict ourselves to s-waves. As in 
Minkowski space, s-wave modes vanish at the origin (pi=oir = 0) = (since 
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ip = if/r) and define a one dimensional problem. Hawking radiation is con- 
cerned with the outgoing reflected part. Because of the existence of Vi=o{r) 
in eq. (13.291) . low energy modes will be reflected back and the problem be- 
comes a usual quantum problem of finding the transmission coefficient. Once 
more for conceptual purposes this complication is irrelevant and we shall put 
Vi=o{r) = 0. Nevertheless, at the appropriate place we shall mention the 
modifications that ensue when its effects are included. In fact, in the last 
analysis though this barrier is a nuisance for mathematics it will prove to be 
of benefit since one will not have to confront the infra-red problem associated 
with myriads of soft Hawking photons. Only frequencies A > 0(1/2M) are 
passed. So the scrupulous reader may waive his qualms in the knowledge 
that what follows is quite rigorous for sufficently large A. 

With the above mentioned simplifications the s-waves obey 

in R (3.32) 
in I (3.33) 







\ d d 






) dudv'^ 






d d 


grp2 


Qj.2 1 





From eq. (13.321) and eq. (13.331) the solutions are the sum of an ingoing and 
an outgoing piece 

= X{v) + E{u) inR 
^ = xiV)+m in I ^^-^^^ 

where continuity on the star's surface impose^ that x(V) = X(v{V)) and 
^(U) = E{u(U)). Moreover, since (p must vanish atr = V — W = 0, one finds 
that everywhere (p takes the form 

^ = x{V)-xm (3.35) 

Before discussing in detail the form of the modes, we note that eq. (13.351) 
already contains information about the division of modes into the producing 
(class I) and non producing classes (class II) as in chapters 1 and 2. In 
the distant past, before collapse begins, we set up vacuum on X~; this is in 
fact Minkowski vacuum since the space is flat on X~. The modes propagate 



^In a realistic four dimensional model there will be some scattering in the star (mixing 
U and V modes) at low frequencies. But at the exponentially large frequencies encountered 
in Hawking radiation the star is completely transparent 
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from X~; hence they are of the form e"*'^'" (once more we are thinking of a 
broad wave packet). They diminish in radius, and either penetrate or do not 
penetrate the star's surface in there inward journey according to the value 
of V (See Fig. fl3.2l) ). Hawking radiation comes about from modes in the class 
that penetrates into the star, but not all of them. Only those that reflect 
in a finite Schwarzschild time {v < vh in Fig. (13.41) ) thereby picking up a 
dependence on u which gives rise to the radiation. 

The asymptotic requirement that there be vacuum on fixes that the 
incoming modes be of the form y?^ = Xlj{v) = e~'^'^'" / y/47ruj . We shall be 
concerned with the behavior of the modes in the vicinity of the light cone 
H that generates the horizon when U{u) reaches its asymptotic behavior 
(eq. (13.251) ). Hence we need the relation between v, u and V,U in the vicinity 
of H only. More precisely we need only the function f (V) in the vicinity of 
V = when H is an infalling light cone since we have already the function 
U{u). Near V = the function V(f) can be approximated hnearly v — vh = 
nV + O(V^) where k is an irrelevant constant which depends on the shell's 
speed at V = 0. Hence in the region where H is an infalling light cone the 
modes take the forms 

C = ^ = ^ (3.36) 

We have used the fact that one may coordinatize the whole space with either 
f or V and used the continuity at the star's surface. After reflection at r = 0, 
H is an outgoing light cone and the in modes read 

= - , = ^= (3.37) 

V Attuj V Anu 

where K{= k,B4M, see eq. (I3.25P ) is a positive constant. Note that the 
approximation 

U = -54Me-"/^*^ = BU (3.38) 

is rapidly excellent since the corrections are of the order e""/^-^. Therefore 
the asymptotic form eq. (I3.37P is also rapidly valid. In this asymptotic regime, 
the Hawking radiation is independent of K because K can be reabsorbed into 
a redefinition of the origin of u. For convenience we choose K = AM and 
drop the phase — e"*"^^^. 

The rest is copying out the results of Chapter 2. We repeat the salient 
facts. Once more because the Bogoljubov transformation does not mix [/'ness 
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and V^'ness, each sector may be handled independently. Concerning the V 
modes, at this point, nothing need be said, the modes being of the form e"*'^'' 
throughout. The in-vacuum stays equal to the out-vacuum i.e. the vacuum 
defined by single particle states in R. 

The physics under present scrutiny is encoded in the overlap between 
the u part of the scattered in-modes given in eq. fl3.37p and the out-modes 



^out 



-iXu 



duo 



^ ,^out I n o'l 



out* 



(3.39) 



f+OO ^ 

ax. = I du^t*{u){-iduW:{U{u)) 



— oo 
+00 



(3xu = / du^T{'^){-^duW:\l^{u)) (3.40) 

J — CO 

Since U{u) is correctly approximated by BU only when the exponential cor- 
rections e^"/^^'^ are negligible, there is no universal form for axw and (3xu} 
when the Doppler shift is small i.e. when u; ~ A, hence there is no steady 
flux in this early stage. On the contrary, when the Doppler shift becomes 
exponential, the corrections fade out, and one does find a universal behavior 
for axuj and (3xuj when u; >> A. This is due to the finite interval in u in 
which these overlaps acquire their value. Hence, for >> A, one can replace 
the expressions for axui and [3xuj by the following expressions, thereby making 
contact with the Rindler-Minkowski Bogoljubov transformation, 

ax. = due^'-e^-'^^-''^ 



2tt ^Juj J- 



'-r(-^4MA) ( - 



'HMX 

^tt2MX 



271 \ U \U , 

fix. = e-'^'al^ (3.41) 

where we have used eqs (12.281) . Hence one finds a steady thermal flux at the 
Hawking temperature = l/SvrM = I/Ph since, for all u » A, one has 

ax. 
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We note that the transitory regime can nevertheless be handled analytically 
in some specific cases, see for instance Appendix |D] and reference [HB], and one 
explicitly verifies that there is an early stage without significant fiux which, 
after a certain w-time, tends exponentially fast to the asymptotic behavior. 
One can then verify that modes with / > do not contribute significantly 
to Hawking radiation. This is seen by comparing the Hawking temperature 
with the height of the centrifugal barrier in eq. (13.301) . 

Let us recall how a fiux arises. As in section 2.4, for a given uj and A, 
the integrals eq. (13.411) arise from the regions around the saddle u* at (see 
eq. ^Mj) 

Re(M*(cj, A)) = 4Mlncj/A (3.43) 



With Gerlach [3H], we call this saddle time the "resonance" time. The width 
of the significant region around u* is independent of uj and given by ~ \JaM/\ 
thereby justifying the replacement of eq. (13.401) by eq. (I3.4ip for >> A. This 
independence of the width allows for the derivation of an exact asymptotic 
formula for the rate i.e. for the fiux. First pick an interval ui < u < ui + Am 
with ui sufficiently large so that the correction term in eq. (13.251) can be 



neglected and Au >> y4M/A. We wish to calculate the flux in R due to 
quanta of frequency A. In this interval the frequencies uj which contribute to 
< a\a\ > (i.e. the number of quanta of frequency A) satisfy 

Whence the rate is given, following eqs (12.551 12.541) . 



lim < a\ax > - — = lim - — / dujiPx^l'^ 

4M /■Ae(-i+^")/4« 



= lim ^ , a ^ 

Au^oo Au JAe«i/4A^ UJ e'^H^ — 12tc 

= {e^"^ -1)-H/2tt (3.45) 

Since at fixed r, one has Au = At, < a\ax > / Au indeed represents a particle 
fiux per unit time. In this way the 6{0) which would come up by applying 
naively eq. (I2.80p is replaced by Aujlix. (This point has been the subject 
of some misunderstanding in the literature. There is no infra-red divergence 
in the number fiux at finite A! To define the total number fiux one should 
take into account the cutoff provided at small A by the potential barrier.) 
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Similarly the total energy flux is 



< T^u >= r^{e^-' - ir = ^Jh' (3.46) 

JO 2,71 iZ 

These results overestimate the flux since they do not take into account 
the suppression of the emission of low frequency quanta A < jSjj^ due to the 
s-wave repulsive barrier. 

Unruh [91] has given a nice mnemonic device to represent these results. 
Namely the steady state of Hawking evaporation (in the early stages when 
the decrease in M due to the radiation is neglected) is simulated by a vac- 
uum state - U(nruh) vacuum. This state is a cross between Schwarzschild 
vacuum (also called Boulware vacuum fl3j) and Kruskal vacuum (also called 
Hartle-Hawking vacuum [13]). Outside the star, U- vacuum is vacuum of 
Schwarzschild v modes (i.e. e'"^^"" / y^TWj) but the Kruskal vacuum of U 
modes (i.e. e~^^^ / \/ Atiuj). This is simply because in the collapsing situation 
the V modes are e"*"^^ and the u modes e^^^ (up to phases and irrelevant con- 
stants). The matching conditions on the surface of the star have universally 
related U to f/, eq. 03.381) . when the shell is close to the horizon H . This is 
the crux of Unruh's beautiful isomorphism. 

One must nevertheless bear in mind the following restriction of the Unruh 
mnemonic. The reference to the Kruskal character of the u part of the in 
modes is after their reflection from r = 0. There is an absolute frame of ref- 
erence in this problem which is given by the movement of the star's surface. 
This precludes the possibility of making arbitrary boosts (in Kruskal coor- 
dinates eq. (13. 6p ) so as to change the character of cj, A mixing in eq. fl3.43p . 
Such boosts correspond to translations in t, so that an event at t = with 
u; = A at resonance in the boosted frame could result say in uj = Ae*^/"^^^. 
The Unruh isomorphism is carried out in a flxed frame e.g. the star at rest 
until t = to. The Gerlach resonance condition uj/\ = e^* 1'^^ which describes 
the Doppler shift at and around u* cannot be boosted away. This restriction 
will play an important role upon confronting the consequences of the very 
high frequencies uj entering into the Bogoljubov coeficients eq. (13.411) as seen 
from eq. (13.431) . 

A contemporary and elegant derivation of the concept of U- vacuum is due 
to Hawking [46j and Damour and Ruffini [27j. Their method is the black hole 
analog of the technique illustrated in eqs. (12. 4 H 12.421) . It is of some interest 
to sketch their method since it makes contact with the tunneling methods of 
Section 1.2. 
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In Eddington Finkelstein coordinates eq. fl3.8l) the d'Alembertian near 



the horizon oustide the star is —dx{-^dx + 2dy) where x 
mode equation is 



d , X d 



2i\]x\{x)e 



-iXv 







dx'2Mdx 

having solutions Xx{^) = const, (f-modes) and the M-modes 



2M. The 



(3.47) 



A 



e{x) 



[X 



AiMX 



9{-x 



-X 



,4iMA 



(3.48) 



One checks that for x > 0, e'^^'x-'^^'^^^i 



*'^") is indeed a u mode in R. 



Completeness in the complete space spanned by the Eddington Finkelstein 
coordinates (collapsing star) requires the term in 6{—x) as well. We now 
require that in this complete space the modes be positive Kruskal frequency. 
Since at v fixed, dx is pure light like (it is in fact —2du) Xx{^) must enjoy 
upper half analyticity in x. The normed modes are thus 



Xx 



V2sinh/?j^A 



\4iAfA 



9{x) 



+ 



V2sinh/5HA 



9{-x) 



—X 



AiMX 



V47rA 



A > 



(3.49) 



as in eq. (I2.4ip . The modes xx are simply the rewriting of the in- modes 
which diagonalize the Bogoljubov transformation (13.391) . Indeed eq. (13.491) 
shows that the modes (p\ = e^^^^xx (— oo < A < +oo) (hereafter called 
Unruh modes) can be written in terms of the Schwarzschild modes y?^"* as 



^x = a^ifT' + Px^T'^* A>0 
<^„A = l3x^T* + ax^l^'^ A>0 



(3.50) 



where we have introduced the modes 



'iXv I 



-X] 



-iiMX 



-X] 



(3.51) 



To make contact with Section 1.2, we note the similarity between the 
differential equation for Xx{.x): [xd^ — 2i'iM\)xx{x) = and eq. (II. 7p . The 
analogy goes quite far. For example we take the Fourier transform of Xx{x) (= 
^a(p)) we get {dpp — 2iX)C^x{p) = 0. The modes ^a(p) having been obtained 
by integration over all x are complete and in fact can serve as a complete set 
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of in-basis states [71]. Their Fourier transform then gives them as a hnear 
combination of out states. Each of these latter hve on one side or the other 
of H, the horizon hne x = 0. 

This is then precisely the same mathematical mechanism used in eq. (11. 150 
to go from the in-state (proportional to 9{u)) to the linear combination of 
states defined in terms of the conjugate momentum to u(= v), a combination 
of 6{v) and 6{—v). The classical exponential approach in a given region of 
space-time to a horizon is translated in each case to a quantum formalism of 
this type. However, in the black hole case, the transcription of the modes to a 
manifest tunneling is ambiguous. We may write the operator ixdx as [11^—^^] 
where U = ^{iAd^ + A'^x) and ^ = -^{iAd^ - A'^x) with [H,^] = -i. 
But the constant A is arbitrary so tunneling in ^ is not only non local in x 
but occurs at arbitrary scales. Nevertheless, the above highlights that the 
collapse has produced pairs living on both sides of x and that these pairs 
are the realization of vacuum fluctuation which have "tunneled" into reality. 
The pair formation is illustrated in Figs (13.51) and (13. 6p in a Penrose diagram 
and in Eddington Finkelstein coordinates. To this end we take into account 
the fact when extrapolated backward in time they bounce off r = 0. Of 
course the 0{U) piece is never seen at finite Schwarzschild time, but only its 
backward reflected piece. As the wave progresses from T_ , this partner of 
the Hawking photon never manages to reflect, but, at t — oo, simply crowds 
into the horizon line extended into the origin. 
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Fig. 13.51 Tiie classical trajectory (stationnary phase) of a pair produced in 
the geometry of a collapsing black hole represented on a Penrose diagram. 
Only the member of the pair which reaches is on shell, the partner falls 
into the singularity. 
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Fig. 13.61 The same as in the previous figure (3.5) but in Eddington 
Finkelstein coordinates. 
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Fig. 13.71 Pair formation drawn in u, v coordinates. 



To compare with tlie accelerated mirror of Section 2.5 we also have drawn 
in Fig. fl3.7p the two pieces of ^\{U) in the extended u^v system. 

From these pictures and the above analysis it is seen that the conversion 
of fluctuations into particles, unlike in Chapter 1 (Fig. (11.21) ). does not occur 
over a well defined space-time domain. This is due to the absence of a scale, 
and it is reflected in the arbitrariness of the parameter A which appeared in 
the definition of 11 and ^ in the above paragraph (we note in passing that the 
introduction of mass does not help since the frequencies of the modes within 
the star rapidly become large compared to any known mass). Nevertheless, 
though tunneling seems inappropriate as a tool in black hole physics (as it 
appears in the present formulation) one should not lose sight of its concep- 
tual content. In chapter 1, particle production could have been understood 
either in terms of tunneling of modes or of their backward scattering in time 
according to the gauge choice. It is therefore useful to define a generalized 
tunneling concept, the conversion of a mode from a domain of virtuality (i.e. 
as a vacuum fluctuation) to a domain of realization (i.e. measurable as a 
quantum of excitation, in a counter for example). The process, in general, is 
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caused by the imposition of an external field. In black hole evaporation the 
process induced by the gravitational field of the collapsing star is summarized 
succintly: virtuality on X~ induces reality on X"*". The precise formulation of 
this aspect of the problem is the subject of Section I3.5[ 

We also note here that we have come upon a rather unexpected fragility 
of the theoretical foundations which have been set to work thus far. Is it 
reasonable to expect modes to remain free and unmixed in their infinite 
excursion from X~ to X+? Indeed, eq. f l3.43p tells us that after a time of the 
order of Am = MlnM the frequencies lo of the vacuum fluctuations converted 
upon resonance into Hawking quanta exceed the Planck frequency. More of 
this in Section 1X71 

3.3 Renormalized Energy Momentum Tensor 
in Unruh Vacuum. 

One of the approaches used to address the back reaction to black hole evap- 
oration consists in solving the semi-classical equations 

G^, = 87c{T^,) (3.52) 

wherein the mean value of the energy momentum tensor is taken as source of 
Einstein's equations. The fundamental assumption is that the source terms 
are (T^^) i.e. that the average solution is the solution of the equations using 
the average source. When the fluctuations are large, or when one wishes to 
evaluate the importance and the consequences of the fluctuations, one has to 
analyze other matrix elements than the mean appearing in eq. (13.521) . This 
is why we shall study in Section 13.51 non diagonal matrix elements of the 
operator T^^. 

The mean energy momentum tensor which appears on the r.h.s. of 
eq. (13.521) is formally inflnite and must be regularized and renormalized be- 
fore starting to solve eq. (13.521) . In this section we deal with this part only, 
in the next section we shall treat the gravitational back reaction in this semi- 
classical approximation. 

Until this time, the complete four dimensional computation of (T^^,) has 
not yet been carried out. What has been achieved is a good approximation to 
(T^jy) in the static spherical symmetric case [IH]. In the evaporating situation 
no analytical expression for the mean energy momentum tensor has been 
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obtained. However we recall that most of Hawking radiation is due to s- 
waves and that, for the s- waves, the problem simplifies a lot if one drops the 
residual centrifugal barrier. (Indeed one ends up with a conformally invariant 
two dimensional theory.) Hence one expects that the essential properties of 
the four dimensional theory are recovered from the effective two dimensional 
one. 

We shall therefore continue with this truncated theory and calculate an- 
alytically (Tp^v) in the Unruh vacuum following ref.[30]. Since the modes in 
the two dimensional model are rescaled by a factor of -\/47rr as compared 
with the original field (see eq. (I3.29p ). the fluxes should be rescaled by 47rr^ 

W") = (3-53) 

This equation applies for fi,!/ = U,V and we shall see that Bianchi indentities 
imply (Tgg) = 0. 

The renormalization of the energy momentum tensor has already been 
carried out, in two dimensions, in the case of a flat background geometry in 
Section [231 and we recall the result, eq. ( 12.861) (extended to the V modes as 
well) 

{Tuu)f,-{Tuu)o = -^f'/'dlf-'l' 

{Tvv)f,-{Tvv)o = ^9"'dlg-"' (3.54) 

Here {Tujj)q means Minskowski vacuum expectation value and U means in- 
ertial Minkowski coordinate. ( ) fg means the average with respect to the 
vacuum defined by the modes e~^^^^^^ / \/ Atiuj and g-^'^sC^) I^At:uj. An im- 
portant property of eqs (13.541) is that they can be inverted so as to express 
(Tff) and {Tgg) in terms of the inverse functions U{f) and V{g): all one does 
is flip the sign and replace /, g by their inverse functions, as is required by 
the reciprocity of these relations. 

In a curved background we shall generalize the subtraction (13.541) by 
subtracting from {T^^{x)) the value {I{x)\T^^{x)\I{x)) calculated from the 
inertial modes at x [65], i.e. those modes which most resemble Minkowski 
modes at x. This is taken on the principle that fiat space is a solution, i.e. 
that the expectation value of T^i, in Minkowski vacuum is zero. The natural 
generalization is to pose that in the local vacuum \I{x)) of local inertial 
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modes, {I{x)\T^^{x)\I{x))ren = 0. For an arbitrary state one then postulates 
that {Tf^i^{x))ren/4:TTr'^ is the gravitational source where 

{T^,yix))ren = {T^y{x)) - {I {x)\T ^^{x)\I {x)) (3.55) 

Furthermore, one must implement this difference of infinities with a regu- 
larization scheme which as in the previous chapters is taken to be the split 
point method e.g. 

{Tuu) = \midudu'{^{U)<i){U')) (3.56) 

To compute {T^u)ren we first need to construct inertial coordinates about 
X. To this end, we express the spherically symmetric geometry, in a "confor- 
mal" gauge {guv = C, guu = 9vv = everywhere) 

ds"^ = -C{U, V)dUdV + r'^iU, V)dn^ (3.57) 

The residual reparametrization invariance U U'{U), V —>■ V'{V) is fixed by 
requiring that the state of the field, in which one wants to compute {T^u)ren 
is vacuum with respect to the modes e"*'^'^, e"""^^. In this way the conformal 
factor C{U, V) encodes the vacuum state of the (p field as in the fiat exemple 
eq. I^M- 

The inertial set U, V based on the U, V set is 

Jua C{Uq,Vq) 

V-V.^f §^ ,V' (3.58) 

Jvo C{Uo, Vq) 



whereupon 



with U, V in eq. fl3.59p functions of U, V obtained by inversion of eqs fl3.58p . 
The coordinates U, V are affine parameters along the radial light-like geodesies 
U = Uq and V = Vq which pass through x = {Uq, Vq). Hence they are the 
inertial set at x which we shall use to define the local vacuum \I{x)) since 
they constitute the most inertial parametrization of the light like geodesies 
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U = Uq, V = Vq {in particular the Christoffel symbols vanish at x in coordi- 
nates f/, V^). One may now apply eqs fl3.54l) in terms of the function U{U). 
Lest the reader have no misunderstanding of the approximate status of the 
following computation, it comes from mode analysis, not from the form of 
the metric. We consider only the contribution of the s-waves and in addition 
drop the residual barrier so as to make the effective field theory conformal 
in the f/, V sector. 

One identifies the Minkowski coordinates of eqs (13.541) which defined the 
subtraction, with the coordinates ?7, V of eq. (13.581) and similarly one identi- 
fies the functions /, which defined the state of the field, with the functions 
U{U),V{y) whereupon 

{Tuu)ren = ~ TTT^^^^ ^IjC ^^"^ 
iZTT 

{Tvv)ren = - ^C^l^ dlC'^!^ (3.60) 

where we have used the reciprocity mentioned after eqs (I3.54p . 

Moreover, in addition to eqs (13.601) . the trace of the energy momentum 
tensor tr T^^ = ATjjv/C = in^cj)^ no longer vanishes even though classically 
it vanishes for a massless scalar field. This is because C is a function of both 
U and V contrary to / in eqs (I3.54p . Indeed, requiring that the renormalized 
energy momentum tensor be conserved 

TVy = ^Tuuy + (C-%y),a= (3.61) 

and the similar equation for T^.y and requiring also that and Tyv be 
given by eqs (I3.60p imposes a trace given by 

{T)ren = "^9^^ (Tuv) ren = TT^R (3.62) 

z47r 

where R = QlnC = 4C~^9(/(9y InC is the two dimensional curvature. Equa- 
tion (I3.6ip is the energy conservation in two dimensions but we emphasize 
that it is also valid in four dimensions if one divides the two dimensional 
fluxes by Anr^ and if one sets Tqq to zero. Then our {Tf^u)ren, divided by 
47rr^, is a legitimate source for the Einstein equations in four dimensions. 

One sees that the result {T)ren = when = is erroneous because one 
has not taken into account the change in the local vacuum upon changing 
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the point x = {Uq,Vq) which leads to the no n- vanishing of the first term 
of the r.h.s. of eq. fl3.6ip . More refined regularizations leads to the result 
directly [30] [10]. One can also interpret the origin of a non vanishing trace 
by remarking that the presence of a curvature within the Compton length of 
a vacuum fiuctuation causes the existence of a term 0{R/w?') which must be 
subtracted to keep the one loop approximation to the effective gravitation 
action finite. 

Note also that the quantum trace given in eq. (13.621) is a pure geometric 
quantity, the same for all state^. Indeed the difference in energy between 
two states (of which eqs (13.541) is a particular example) is given by conserved 
/STuu and ATyv-, hence with AT[/y = 0. 

The above procedure of subtracting {I\T^y\I) can in principle be applied 
in the original four dimentional theory and be generalized to fields with mass 
(see ref. [61]). A number of alternative procedures have been devised for renor- 
malizing the energy momentum tensor (see ref. [10] for a review). Happily 
they are all equivalent: indeed Wald [98] has shown that if the renormalized 
energy momentum tensor satisfies some simple and natural conditions it is 
completely fixed (see nevertheless ref. |17j). 

We now apply eqs (13.601) to our problem. As a warm-up let us start 
with Boulware vacuum. This is the case one would have if the star were 
eternally static at some fixed radius (= Rq) greater that its Schwarzschild 
radius (= 2M) and we shall have in mind the case Rq — 2M << 2M. 

In the Schwarzschild region, the modes e~^^^/ VAttX and c'^^'^/VAttX de- 
fine the Boulware (B) vacuum (e.g. usual vacuum at r = oo) . The con- 
formal factor C is (1 — 2M/r) in the Schwarzschild coordinates u,v de- 
fined in eq. (13.51) . Since C{u,v) = C{y — u) we see immediately from eqs 
(I3.60p that {Tuu)b = {Tvv)b thereby implying no fiux, {Trt)B = 0, as it 
should be for a static situation. Moreover since lim^^oo drC (r) = we 
have {Tuu)b = {Tyv)B = at r = oo. This is as it should be because the 
Schwarzschild metric is asymptotically fiat and because the Schwarzschild 
modes are identical to the usual Minkowski modes for large r. 

It is also very easy to obtain (T^m) b in the vicinity of r = 2M. Indeed, 
near r = 2M, 

C(r) ~ -e'-*/^^^ = _e(''-")/4^^ (3.63) 

This is true even in the full four dimensional theory wherein the trace anomaly is a 
quadratic form of the Riemann 4-tensor. 
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so that from eqs fl3.6UI) 



lim (T„„)b = {T,v)b = -(l/127r)(l/64M2) = -(7r/12)r| (3.64) 

r— >2M 

i.e. minus the asymptotic Hawking flux eq. (13.461) . The origin of this negative 
energy density hes on the Rindler character of the Schwarzschild geometry 
expressed in the Schwarzschild coordinates u,v in eq. (13.631) . Indeed, the 
Minkowski metric expressed in Rindler coordinates (eq. (I2.19P ) reads 

ds"^ = -dUdV = -e"(^-")rfMciw. (3.65) 

Thus the specification of being in Boulware vacuum becomes, near r = 2M, 
equivalent to being in Rindler vacuum in Minkowski space (see eqs. fl2.44[ 

IZI09D). 

It is to be observed that (Tuu) b being finite on the future horizon m — > cxo 
{U — > 0) leads to singular values of (Tmj) b on H: {Tjju) b = {du/dU)^{Tuu) b = 
{4:M/UY{Tuu)b- Similarly {Tvv)b blows up on the past horizon (V = 0) 
which exists in the complete Schwarzshild space. If this situation would 
pertain to the collapsing case (i.e. to the true physical state of affairs) one 
would arrive at a catastrophic situation in that as [/ — 0, the values of (Tuu) 
would tend to — oo; and it is {Tuu) which is close to that measured by an 
inertial observer near the horizon (since ds^ = —e^^dUdV on the horizon, 
see eq. (13. 7p ). The accomodation to this singular behavior is a remarkable 
feature of black hole evaporation. 

Indeed, in the collapsing case, in the Unruh vacuum, whereas the v- 
modes remain e'"^^"" / VAttX in the outer Schwarzschild region the u-modes 
rapidly behave hke e~^'^^ / y/Anuj, see eqs. (I3.37[ I3.38P and the discussion 
after eq. (13.441) . Without calculation we see from eqs (13.601) that near the 
horizon {Tuu)u vanishes (where the subscript U refers to Unruh vacuum). 
Indeed space is regular near the horizon and the inertial modes differ only 
slightly from the Kruskal modes. Therefore {Tuu{f = 2M))^ is finite and 
{Tuu)u = iU /AMY{Tuu)u vanishes quadratically on the horizon f/ = 0. The 
former singularity is thus obliterated. The quadratic vanishing is necessary 
for having no singularity for a free falling observer since the combined effect 
of the gravitational and Doppler shifts already encountered in eq. (13.261) is 
always present for any inertial trajectory crossing the future horizon. In refer- 
ing once more to the isomorphism between this situation and the Minkowski- 
Rindler case, we see that the u part of the state in the U-vacuum behaves 
like a regular Minkowski state on the horizon f/ = 0. 
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One can now either proceed by calculation to find {Tuu)u, using eq. 03.521) 
to go from Boulware to Unruh vacuum (with g{v) = v; f{u) = U{u) = 
— 4Me~'"/'^^^ see eqs. (13. 6[ 13.381) ). or easier yet (and perhaps more physical) 
appeal to the equation of conservation (13.611) . Indeed {Tuu)u differs from 
{Tuu)b by a function of u only. But this function of u must be a constant 
because we are in a steady state characterized by a rate eqs. (13. 43^ I3.44p . 
Whence 

{Tuuir))u = {Tuu{r))B - {T^u{2M))b (3.66) 

i.e. the constant is fixed by {Tuu)u = at r = 2M. Hence, from eq. (I3.64p . 
the constant fiux T\ = Tuu — T^v is given, in Unruh vacuum, by {Tuu{r = 
oo))u = (7r/12)T^ which is the thermal fiux at the Hawking temperature as 
in the mode analysis of Section 13.21 

The behavior of (T^^) for finite r is obtained from eqs (I3.60p with C = 
(1 - 2M/r) and dr* = {dr /dr*)dr: 



{Tvv)u — {Tvv)b — {Tuu)b 

TT .,48M^ 32M3 
12 " 



TkTM^ - (3.67) 



and 



,^ , 7r_,, 4M 12M, 

{Tuu)u = -Tl(l )2(1 + + (3.68) 



wherein the quadratic vanishing at r = 2M and the asymptotic behavior are 
displayed. 

All this is without back reaction i.e. without taking into account the 
decrease of M in time. But the energy conservation (the Einstein equations 
at large r) dictates that there is a necessary backreaction wherein 

(T.)u--^J^r (3.69, 



dt ' 12 V87rM 

We recall that the fiux should be integrated over the sphere and that 
our T^y has been rescaled by 47rr^, therefore eq. (13.691) is the usual four 
dimensional expression for the mass loss. Upon correcting for the existence 
of the potential barrier and adding the contribution of the higher angular 
momentum modes / > 1 one finds 

dM 1 ^ ^ 
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where ^ depends on the spin of the radiated field [HH]- Then, if one assumes 
that at later times the rate of evaporation is given by the same equation, the 
decay time for complete evaporation is Tdecay = in Planckian units. 

In eq. fl3.69p the sphere over which one calculates can be situated at any 
value of r since the fiux is conserved. But the interpretation of the integral 
changes with r. At large r one has {T[)u = {Tuu)u- This is a traditional 
positive energy outfiow. Instead, near the horizon {T^)u = —{T^v)b since 
{Tuu)u vanishes there. How these two properties contrive to modify the 
metric and describe a black hole with a slowly varying mass parameter is the 
subject of the next section. 

3.4 The Semi-Classical Back Reaction 

Up to this point in our primer we have concentrated on the response of 
matter to a fixed background geometry, that of the collapsing star. The 
dynamical response to this time dependent geometry is the emission of energy 
to infinity accompanied by an accommodation of the vacuum near the black 
hole's horizon. Indeed at the horizon the mean fiux is carried by (T^v) only. 
As emphasized at the end of the previous section, the rate of change of mass 
is {Trt)\ fixed r where the value of r is arbitrary. Hence near the horizon the 
description of the mass lost by the star is given in terms of a halo of negative 
energy which accumulates around the horizon. 

The above discussion then suggests the shrinking of the area of the hori- 
zon. But how does a significant reduction of the area, hence a large change 
in the metric, affect T^jy, and how in turn do these affect the metric? The 
complete answer to this question (i.e. wherein one takes into account the full 
quantum properties of the operator T^^) is the subject of present research 
and is far from resolution -very far indeed in that its ultimate elucidation 
may well entail (or lead to) the quantum theory of gravity (see Section 13.71 
in this regard). 

Until the present time the only quantitative treatment available is in 
the context of the semi-classical theory i.e. gravity is classical and Einstein's 
equations are driven by the mean value of T^,^ (see eq. fl3.52p ). Since the latter 
is a function of gfj_u and its derivatives, one has highly non trivial equations 
to solve. In what follows we shall report on what is known concerning this 
enterprise. We first give a semi quantitative description of what happens. It 
is then followed by more rigorous considerations and the presentation of the 
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properties of the evaporating geometry obtained by numerical integration of 
a simplified model. 

What are the features that one wishes to display? The first concerns the 
validity of eq. fl3.69p at later times. Is the mass loss at time t determined by 
M{t), the mass at that time, through Hawking's equation? In other words is 
the system Markoffian? One finds that the answer is yes, and that it is due 
to the continued negative energy flux across the horizon. This results in an 
outside metric which is determined by M{t) and not by history. The second 
question concerns the distribution of this negative energy. Is it a halo, or is 
it distributed homogeneously within the collapsing star so as to give rise to a 
sort of "effective" star of mass M{t)l One finds the former. There is a halo 
which accumulates in a region enclosed between the surface of the star and 
the horizon, which we now explain on qualitative grounds. 

In order to describe the geometry near the horizon it is necessary to choose 
a coordinate system. A convenient choice is the Eddington-Finkelstein set 
(r, t>) introduced in eq. (13. 8p . The advantages of these coordinates are: 

1) they cover both sides of the future horizon; 

2) in the case of a purely ingoing light like flux, the metric 

ds^ = -{l- '^^^)dv^ + 2dvdr + r^dn' (3.71) 
r 

where the only parameter is 

M{v) = J'' dvT^^ (3.72) 

is an exact solution of Einstein's equations (the Vaidya solution). We remind 
the reader that we are working with energy momentum tensor rescaled by 
Anr^, see eq. (I3.53p . In the evaporating situation near the horizon the flux is 
entirely carried by T^y hence 13.711 is valid in the vicinity of the horizon and 
can be used to describe the physics there. 

We are interested in the behavior of outgoing light rays because the modes 
which give rise to Hawking radiation flow along these light rays. Let us 
first recall what values of r a light ray, reflected from r = 0, visits without 
back reaction (i.e. with M = Mo). Photons which emerge from the star at 
r = 2Mq + e (with e > 0) proceed to infinity after adhering to the horizon 
r = 2Mq for a while (a v lapse of order Mq ln(e/2Mo)). Those which emerge 
at r = 2Mq — e also adhere to the horizon for a while before falling into the 
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singularity. These trajectories describe the locus 

u = const = v-2r* = v-2r- AMq In |r/2Mo - 1| (3.73) 

for r near the horizon. The main point is the existence of a horizon at 
r = 2Mq which separates the out's from the in's, those which escape from 
those which are trapped and crash into the singularity. 

Now suppose that the mass decreases due to the absorbtion of negative 
energy: (T^y) ~ — Then the locus which separates the out's from the 
in's is expected to shrink. Thus outgoing light rays feel a slowly diminishing 
gravitational potential and some of those which previously were trapped may 
now escape after having been sucked in for a while. The locus where they 
cease to fall in and start to increase in radius is called the apparent horizon. 
To describe this quantitatively we consider the equation for outgoing light 
rays when the mass is varying, in the metric 13.711 

dr _lr - 2M{v) 
dv 2 r 

The apparent horizon is the locus where dr/dv = 0, i.e. 

rah = 2M{v) (3.75) 

And it shrinks according to eq. f l3.72p i.e. dvah/dv ~ —1/rlf^. However 
the apparent horizon does not separate the geodesies which will ultimately 
fall into the singularity from the escaping ones since it is not an outgoing 
geodesic (for an evaporating black hole, when is negative, it is time like 
|47j). It is nevertheless very close to the event horizon (the light ray which 
does separate the two classes) . An estimate of the radius of the event horizon 
is the inflexion point of the outgoing geodesies d'^r/d'^v = 0: 

reh{v) = 2M{v) + SMdM/dv = Vahiv) - 0(1/M) (3.76) 

Therefore all the degrees of freedom which are in the whole region r < 
2M{v) — 0{1/M) remain, in this semi-classical treatment, inaccessible to 
the outside observer. This inaccessibility implies that the loss of mass does 
not come from the evaporation of degrees of freedom from inside to outside. 
Rather, it is due to the accumulation of negative energy. 

The slow rate of change of Vah suggests that the processes which occur 
in this geometry are the same as in the case with no backreaction and with 
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(3.74) 



M equal to M{v). For instance an important property of the outgoing light 
rays in the geometry 13.711 is that they are given, in good approximation, by 
the usual formula, eq. (13.731) . with M replaced by M{v) 

v-2r- 4:M{v) ln(r - 2M{v)) = u (3.77) 

provided 

M » r - 2M(w) » 1/M (3.78) 

where the restrictions come from the rate of change: dM/dv ~ 1/M^. Hence 
when this is valid the out-going geodesies are a scaled replica of what happens 
without back reaction with M replaced by M{v). 

Based on the above one may conjecture that the radiation emitted at 
time V is also a scaled replica of the radiation emitted in the absence of 
back reaction with M = M{v). To understand why the evaporation is only 
controlled by M{v) (and therefore why the past history of the black hole 
plays no role) and to have a full appreciation of the processes which occur 
around the apparent horizon (where eq. (13.770 is not valid) one must resort 
to a more detailed analysis. This is done in the following paragraphs. The 
upshot is that all that has been discussed qualitatively is correct to 0(1/M). 

[ Following Bardeen^ and York[_96j who followed up the ideas of Hajicek 
and Israel [12] we shall use the metric 

ds^ = -e^^{l - 2m{v, r) /r)dv^ + 2e^dvdr + r'^dVt^ (3.79) 

which describes a general spherically symmetric space-time. In these coordi- 
nates Einstein's equations are 



rpv 

V 

Trr/r (3.80) 

Note that ip is defined only up to the addition of an arbitrary function of v 
corresponding to a reparametrization of the v coordinate. For simplicity we 
shall suppose that the right hand side of (13.801) is given by the two dimensional 
renormalized energy momentum tensor discussed in section 13.31 (see eq. 13.601 
and l3.62] ). However the proof is general provided M >> 1 . The r.h.s. of l3.80l 



dm 

dv 
dm 

dr 
dr 
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can be taken to be the full 4 dimensional renormalized energy momentum 
tensor. 

We shall proceed in three steps following [62]: first we shall suppose 
that the renormalized energy momentum tensor resembles the renormalized 
tensor in the absence of back reaction. We shall then show that under this 
hypothesis the metric coefficients in (13.791) are slowly varying functions of r 
and V. Finally we shall solve adiabatically the Klein Gordon equation in this 
slowly varying metric and show that the renormalized tensor indeed possesses 
the properties supposed at the outset thereby proving that the calculation is 
consistent. 

Our first task is to obtain estimates for T^^ both far and near the black 
hole. We begin with the former. We suppose that when r is equal to a few 
times 2m (say r = 0(6m)) there is only an outgoing fiux Tuu{r >> 2m) = 
Lh{u) where Lh is the luminosity of the black hole. This is justified since in 
the absence of back reaction the other components of T^i, decrease as large 
inverse powers of r. For instance the trace anomaly, in our 2 dimensional 
problem decreases as M/r^ in the static Schwarzschild geometry. We shall 
also suppose that Lh is small (i.e. LhM'^ = 0(1)) and varies slowly. Hence 
when r > 0(6M) an outgoing Vaidya metric is an exact solution of Einstein's 
equations 

ds^ = -(1 - 2M{u) /r)du^ - 2dudr + r^dil^ 

M{u) = rdu' Lh{u') (3.81) 

The change of coordinates from (I3.80p to fl3.8ip is obtained by writing the 
equation for infalling radial null geodesies in the metric (13.811) as 

2dr 

Fdv = du+ — — (3.82) 

1 - 2M{u)/r ^ ' 

where F is an integration factor. Upon using (13.821) to change coordinates 
from the set (m, r) to (f , r) one finds that = F and m{r, v) = M{u). Hence 
when r > 0(6M) the r.h.s. of (I3.80p is given by 

T\ir>0i6M)) = -cHh 
-Tl{r>0{6M)) = 2LH/{l-2M/r) 
Trr{r > 0{6M))/r = ALh /r{l - 2M /rf (3.83) 

We now estimate T^^ near the horizon (i.e. near Tah = 2m{rah,v)) by 
assuming that the energy momentum tensor measured by an inertial observer 
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falling across the horizon is finite and of order Lh- Near r = one may 
neglect g^^ = — e^'^(l — 2m/r) and use the reparametrization invariance of v 
to choose ip{rah, v) = whereupon the metric becomes ds"^ ~ 2dvdr + r'^dVL^. 
Hence near rah-, r and v behave like inertial light like coordinates (ie. the 
proper time of an inertial infalling observer near the apparent horizon is 
cv + c~^r where c is a constant which depends on the precise trajectory of 
the observer). That the energy momentum is of order Lfj near the horizon 
is reexpressed as (for the components and Trr) 

T:{r^rah) = 0{Lh) 

T„{r-raH) = 0{Lh) (3.84) 

Where we have used the inertial character of the set f , r near Tah- In addition 
we shall determine T^(r ~ Tah) by making appeal to the conservation of 
energy 

T\r + = (3.85) 

where = 0{Lh,v)- Integrating the conservation equation from r = 2m 
to r = 0{6m) yields near the horizon in terms of its value where (13.831) 
is valid. Putting everything together, near the horizon we have 

-LHe^ + 0{mLH,,) 
0{Lh) 

0{LH/r) (3.86) 

As announced all metric coefficients vary slowly if is small and varies 
slowly. Integrating the equation for yields ~ r^'^ for all r > rah{v)- 
Hence ijj can safely be neglected up to distances r = 0{e^^^^). From now on 
we suppose for simphcity of the algebra that = 0. 

In order to calculate the modes and {T^i,)ren we must first investigate the 
outgoing radial nul geodesies in the metric (13.791) with ip = 0. As mentioned 
in \'S.77\ away from the horizon these geodesies are a scaled replica of the 
geodesies in the absence of back reaction. Near the horizon their structure 
is complicated by the distinction that has to be made between the apparent 
and event horizon. We recall (equation I3.75P that the apparent horizon is 
the locus where outgoing geodesies obey dr/dv = 0, therefore solution of 
Tahiv) = 2m{v,rah{v)). 



dm 

dv 
dm 

dr 

dil) 

dr 
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The event horizon r^hiv) is the last hght ray which reaches T+. It satisfies 
the equation of outgoing nul geodesies 



dreh 1 Teh - ^m^reh, v) 
dv 2 Tfj 

Setting reh{v) = Vahiv) + A(t;) one can rewrite 13.871 in the form A(t>) = 
2{rah{v) + ^){r-ah),v + ^,v) + 2m{v,reh{v)) - Vahiv). Solving recursively 
one obtains an asymptotic expansion for A the first term of which is A = 
rah{v){rah),v - -1/M (see equation [S^SD. 

To obtain the trajectory of the outgoing nul geodesies we change variables 
to the set (f , x = r — reh{v)) i.e. x is the comoving distance from the event 
horizon. In these coordinates the metric eq. (I3.79p becomes (using 13.871) 

ds^ = _ M^,^efe + x)x ^^2 ^ 2rft;da; + r^dn^ (3.88) 

This metric resembles the Edington-Finkelstein metric in the absence of back 
reaction in a crucial way. To wit gvy{x,v) vanishes on an outgoing null 
geodesic, the event horizon x = 0. Using this form for the metric it is now 
straightforward to obtain the outgoing null geodesies. Indeed when x « r^h 
the equation for radial outgoing nul geodesies can be solved exactly to yield 
an exponential approach to the horizon of the form v — 2\\ix = f{u) where 

' , 2m(v, reh(v)) 

dv ^ ; , 3.89 

This motivates the following ansatz for the outgoing radial null geodesies 
valid in all space time outside the collapsing star 

~ ~ 2\nx + 6= / —— + D (3.90) 



reh{v) J 2m{u' 

with D a constant of integration. This should be compared with the solution 
in the absence of backreaction given in eq. fl3.73p . The function of u on 
the r.h.s. of (13.901) has been written as du' /m{u') for dimensional reasons. 
The quantity 5 is of order 0{Lh[Mx + x^)/M^) for all x (this is shown by 
substitution of (13.901) into the equation for radial nul geodesies and integrat- 
ing the equation for 5 along the geodesies u = const). The function rh{u) is 
determined by requiring that the variable u in equation (I3.90p be the same 
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as in the Vaidya metric equation fl3.81l) . The difference M{u) —m{u) is then 
of order 0{MLh)- This is found by using fl3.90p to change coordinates from 
the set (v, r) to the set (u, r) at the radius r = 0(6M) where f l3.8ip is vahd. 

Equation (13.901) is sufficient to prove our hypothesis, to wit that the flux 
emitted is 0(M~^) and that the energy momentum tensor is regular at the 
horizon. Indeed, in our model, the flux emitted is given by eq. fl3.54l) : 

(T,,(m, J+)) = {l/127r){dU/duY/^dl{dU/du)-^/^ (3.91) 

where the derivatives are taken at fixed v. The variable U, defined in eq. 
(I3.13p . labels the outgoing geodesies as measured by an inertial observer 
inside the star. The jacobian du/dU is calculated by remarking that at the 
surface of the star the derivative at fixed v is dU / dr\v=vstar = —2. Hence 
differentiating 13.901 yields 

du/dU\y = {dr/dU)\y=^star{du/dr%=ystar 

= {l/2)Am{u)/x\,=,,tar = -^m{u)KU-Ueh) (3.92) 

Hence {Tuu{u,I^)) is equal to (7r/12)T^(u) where 

Tniu) ^ -4-- = + OiLj,)) (3.93) 

87im[u) SttM{u) 

is the Hawking temperature at time u when the residual mass is M{u). 

The calculation of (T„„(f , x))ren every place (and not only on X"*") is a 
slight generalization of the above calculation. The renormalized energy mo- 
mentum tensor is given by 

{Tuu)ren = {dU / duf {I / l2-K){dU / duf/"^ dl{dU / dU)-^'^ (3.94) 

where (j {u, v) is the inertial coordinate introduced in equation 13.581 Since 
U{u, v) is an affine parameter along radial nul geodesies v = const one obtains 
that U{u,v) = x{u,v). Hence 

dU/dU\y = dx/ du\vdu/ dU\y (3.95) 

which upon differentiating eq. I3.91l and using [3^921 is found to be finite on the 
horizon. Hence as in section [331 {Tuu)ren vanishes quadratically at the event 
horizon lA = Ueh- Thus the mean (T^j,(r, t>)) is the one computed without 
backreaction, in the collapsing geometry of a star whose mass is M{v) and 
this is valid as long as dyM{v) << 1.] 
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These conclusions have been verified numerically in a model [75] wherein 
the sources of the Einstein equations are a classical infalling fiux of spherical 
light-like dust (i.e. falling along v = constant) and the energy momentum 
tensor given once more by eqs. (13.601 and 13.621) properly rescaled see eq. 
fl3.53p . Since this quantum energy momentum tensor is a local function of 
derivatives of the conformal factor C, one is lead to a new set of differential 
equations of the second order which describe dynamically the evolution of 
space-time. The main advantage of this numerical integration is to provide 
the whole geometry from the distortion of Minkowski space-time due to the 
infalling shell up to the complete evaporation of the collapsed object. 

There is little point here in going through the algebraic formulation for 
this model in view of the general considerations which have just been set 
forth. We do wish however to point out one interesting feature of this work. 
In view of the T^,y used, the system of coordinates is that of eq. 13.571 It is rel- 
evant to remark that there has been an independent mathematical approach 
to the black hole problem inspired by some elements of string theory called 
the dilatonic black hole in 1+1 dimensions. It resembles in form the above- 
mentioned calculation. Indeed the dilatonic scalar field 0(f/, V) is played 
by lnr{U, V) [19] [82] [81] [77] . The original hope was that the dilatonic black 
hole model could be developed into a complete quantum mechanical theory. 
Unfortunately these hopes have not been realized and at the present time 
our knowledge is restricted to the semiclassical theory. Since we are now 
possessed of a more physical semiclassical theory we shall not present this 
work in this review. 

We now present the numerical results in a series of figures. 
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Fig. 13.81 The outgoing nul geodesies (u = const ) in the geometry of an 
evaporating black hole, depicted in r, v coordinates. 

Figure fl3.8p sliows what liappens in terms of the coordinates r and v (see 
eq. 13.711 and 13.791) wherein the axis v is drawn horizontally and the axis r 
vertically. The units are Planckian and Mq is equal to eight. The classical 
trajectory of infall is not designated but it is a packet centered around the 
line f = 25. The contour lines which are shown are equally spaced u = const 
outgoing light rays (that is equally spaced on X"*", here represented essentially 
by r = 20). The event horizon appears as a thick line because it is where 
all these outgoing geodesies accumulate. The geodesies which fall into the 
singularity are not represented. They would lie in the white zone which lies 
beneath the event horizon. All of these geodesies emanate from the point 
of reflection r = (that is well within the star). The first ones continue to 
increase in r on their voyage to whereas the later ones first increase in 
r, then decrease in the evaporating geometry, then increase once more. The 
locus where they reexpand is the apparent horizon. It is interesting to remark 
that the paths of photons in the white zone which begin their trajectory at 
r = and fall back towards r = resemble those of a closed Robertson 
Walker universe. Indeed the maximum radius encountered by the outgoing 
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light rays in this zone is 2Mo and then the radius diminishes according to the 
evaporation rate. At the apparent horizon rah = 2M{v) one finds the throat 
which connects the interior region to the external one. Thus at the end of 
the evaporation, one has two macroscopic (smooth i.e. wherein the mean 
geometry is far from the Planckian regime) regions connected by a throat of 
Planckian dimensions. This is the situation that precedes the splitting of the 
geometry into two disjoined regions (universes) with a change of topology 
[48j . From the semi-classical scenario this seems unavoidable and confirms 
that information is forever lost to the outside universe. 

In this theory there is a singular space like line r = Tq, = 0(1) (which 
is not represented in figure 13.81) . This line is singular in that the dynamical 
equations loose meaning at this radius. When the apparent horizon reaches 
Ta (i.e. when the residual mass of the hole is ra/2), one has to stop the 
numerical integration. Thus the semiclassical model cannot describe or give 
any hint about the endpoint of the black hole evaporation. This is as it 
should be since at the end of the evaporation, the mean curvature reaches 
the Planckian domain where the semiclassical treatment has no justification 
whatever. 

Two other figures of interest are Figs (13. 9p and (13.101) . These represent 
the same geometry in other coordinates. 
In Fig dSlD 
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Fig. 13.91 The r = const lines in the geometry of an evaporating black hole, 
depicted in U,V coordinates. 

we present the geometry in tlie W, V coordinates defined in eqs. (13.131) . These 
are the inertial coordinates in the Minkowski region inside the spherical in- 
falhng shell. We have represented the contour lines of r = constant since 
they offer a convenient visualization of the evaporation process. Indeed, in- 
side the shell, where there is no matter one has Minkowski space time in 
which r = constant are always time like straight lines. Outside the shell, 
in the absence of back reaction, one would have Schwarzschild space with 
mass Mq and the apparent horizon coincides with the static event horizon at 
r = 2Mo which separates time like from space like r = constant lines. 

As we have explained, in the presence of back reaction, the evaporation 
process is accompanied by the shrinking of the apparent horizon rah{v). This 
horizon is the exterior boundary of the trapped region wherein r = constant 
lines are spacelike. The inner boundary of the trapped region (the other 
solution of dr^r = at fixed U) lies within the shell and is space like. [Indeed 
because the flux T^y is positive here, due to classical matter falling into the 
black hole, the apparent horizon is space like. In the evaporating situation 
T^^ is negative and the associated apparent horizon is time like.] In this 
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figure we have also drawn the singular space-like line r = which meets the 
apparent horizon rah{v) at V ~ 52. 

In the U, V coordinates, the whole evaporation period is contained in a 
tiny U lapse (—14 < U < —12.8). This is why we have presented in Fig 
f l3.10p the evaporating geometry in the m, v coordinates. 



event horizon (r^J 




V 



Fig. 13.101 The r = const lines in the geometry of an evaporating black 
hole, depicted in u, v coordinates. 



These later are the inertial ones at r = oo where they are normalized by 
V — u = 2r. One sees the dramatic effect of the exponential jacobian eq. 
fl3.92p relating utoU which blows up the region between Tah and r^h- In the 
M, V coordinates, the apparent horizon appears as an almost static line [where 
static is defined as follows : if a mirror is put along it, two inf ailing light rays 
separated by Af (i.e. by At at fixed r) will be reflected with the same Au 
(i.e. by the same At).] The property of staticity is obviously satisfied by the 
r = constant lines in a static geometry but strictly speaking no longer in the 
evaporating situation where, as we have just described, r = constant passes 
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from space-likeness back to time-likeness. We also note that the singular line 
at r = Tq, is not present in this figure since it is beyond the last u-line (i.e. 
the event horizon reh{v)) at which the apparent horizon meets this singular 
line. 

In Fig (13.101) we have drawn an extra line which sits outside the trapped 
region. This line designates the locus where (T„u(r, u)) reaches half of its 
asymptotic value ((T„„(cxd, m)))). Being well outside the trapped region, this 
proves that the fiux is concerned with the external geometry only and is char- 
acterized by the time dependent mass M{u) and not by the whole interior 
geometry. This shows also that the infalling matter is not at all affected by 
the evaporation process since the infalling matter is in the causally inacces- 
sible past of the places where the mean fluxes build up. 

The semiclassical theory which we have presented in this section is a math- 
ematically consistent and well understood theory which predicts that black 
holes evaporate following the law dM/dt = —^M~'^. However the validity of 
the semiclassical theory even when the curvature is far from the Planck scale 
is far from obvious, owing to the fact that very small distance scales are in- 
voked in order to derive the solution. Indeed the jacobian dlA/du = e~"/'^^(") 
used in obtaining the flux at inflnity in eq. l3.91l becomes exponentially small, 
ie one makes appeal to the structure of the vacuum inside the star on expo- 
nentially small scales in order to derive the radiation at later times. Another 
related problem is that the distance between the event horizon and the appar- 
ent horizon is 0(1/ M{v)) which for macroscopic black holes is much smaller 
than the Planck length. [This distance is an invariant: it is the maximum 
proper time to go from one horizon to the other. This is seen by writing the 
metric eq. (13.791) near the horizon in the form ds"^ = 2drdv + r'^dQ'^ since we 
can drop the term in (r — 2M)/2Mdv^]. We shall discuss in more detail this 
problematic aspect of the theory of black holes in section 13.71 

3.5 From Vacuum Fluctuations to Hawking 
Radiation 

In the previous sections we have described the mean value of the energy mo- 
mentum tensor and its effect on the background geometry. In this section 
we consider fluctuations. More precisely we describe the fleld conflgurations 
which evolve into a particular Hawking photon using the weak value formal- 
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ism. This will provide us with the explicit history of the creation of the 
particle out of vacuum. It will also provide us with the matrix elements of 
T^i, which play the role of the mean when calculating backreaction effects to 
S matrix elements wherein the final state contains this specific pair. 

As in the case of the accelerating mirror, we follow the history of the 
vacuum fluctuation associated with the emission of this photon. What follows 
is a resume of Sections 12.51 12.6.31 and the forthcoming paragraphs. 

Consider a created particle in a wave packet emitted at retarded time mq, 
of frequency A (where A = 0{M~^)) with width Am = 0(A~^). From the 
Gerlach resonance condition u*{uj,\) = 4M \n{uj / X) , the frequency uj of the 
mode comprising the fluctuation within the star that is converted into this 
Hawking photon is = 

0(^e«o/4M)_ This fluctuation is set up on X and is 
composed of three parts. Let f = f be the backward reflected light cone of 
the horizon U = (see Fig. 13. 2p . Then, as in Section [275| the fluctuation, 
represented as a packet, straddles v = vh- It is spread out on either side of 
V = vh with a spread Af = 0{uj~^) = 0(AMe~'"°/^*^). The part with v > vh 
has positive energy density. The fluctuation with v < vh also has a positive 
energy hump, as well as an oscillating broader distribution of energy. This 
latter is net negative and the total energy of all these contributions vanishes 
as behooves a vacuum fluctuation. 

As in Section [231 the two different contributions, 9{±{v — vh)) possess 
very different future destinies. The piece for t> — < 0, reflects off r = 0, 
travels through the star and gets converted into a Hawking photon. How- 
ever unlike the mirror it is only at large radius that space gets flat and that 
one gets a true on mass shell quantum. Once the photon gets out of the 
star, the oscillating net negative piece of the fluctuation becomes negligible 
because of the time dependent Doppler shift [dlA/duY encountered in con- 
verting Tuu to Tuu- The emerging photon is a net lump of positive energy, a 
propagating outgoing photon. On the contrary, the piece of the fluctuation 
for V — Vh > ^ cannot get to the horizon in finite Schwarzschild time. It 
approaches the horizon exponentially near the center of the star and there it 
sits carrying net positive energy of 0(M~^e"'"''^^). However, as in the case 
of the accelerating mirror, the average energy carried by all fluctuations is 
zero. Thus an observation of the absence of a photon emitted aX u = uq is 
associated with an "anti-partner" fluctuation carrying negative energy near 
the center of the star. So on the average the fluctuations carry no energy. 
The problems raised by the gravitational back reaction to the exponentially 
large fluctuations near the horizon will be discussed in Section 13.71 
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The detailed evaluation of the weak values of T^^j for the effective two 
dimensional model proceeds exactly as in the accelerated mirror problem. 
Indeed the two dimensional part of the black hole geometry depicted in the 
M, V coordinate system (see Fig. 13.41) is almost identical to the accelerated 
mirror problem (see Fig. 12. 3p . The role of the mirror is played by the center 
of the star r = (see the reflection condition eq. fl3.35p . The only difference is 
that the conformal factor of the metric is trivial in the mirror problem {ds^ = 
—dudv) whereas it is nontrivial in the black hole problem. This introduces 
some complications when computing the renormalized energy momentum 
tensor which we now address. But it does not affect the fluctuating part 
of T^y since we are considering the simplified s-waves which are conformally 
invariant. 

In Section [331 the mean energy momentum tensor of the truncated model 
was computed. In order to calculate any matrix element (and not only the 
mean) of the renormalized energy momentum operator we need a slightly 
more general formalism wherein the renormalized energy momentum opera- 
tor is written as 

T^Ax)ren = T^u{x) - {I {x)\T^,{x)\I {x)) I (3.96) 

Here T^^{x) is the bare energy momentum operator, |/(a;)) is the inertial 
vacuum, I is the identity operator, (J(x)|T^,^(x)|J(x))is the expectation value 
of T^i, in the inertial vacuum which is conserved upon including the trace 
anomaly. Note that in this expression the trace anomaly is entirely included 
in the second term so it is state independent as required. A convenient 
reexpression of eq. (13.961) which isolates the mean of T^^, from its fluctuating 
part is 

Tfj^ii ren -T^ii ■ ~\- {^(^in\Tp,y\^in) (-^(''') I -^/tti' | -^(i^) ) ) 

-T^u- ~l~ (Ojri I -^jui^ I Om)ren (3.97) 

where {Qin\T^u\^in) is the expectation value of T^^ in the Heisenberg vacuum 
|0i„) and :T^j,i, : is the energy momentum operator normal ordered with respect 
to this Heisenberg vacuum. 

We are now in position to recopy the results of Section [2751 We first recall 
(see eqs fl3.49|3.50|3.51l) that to each out Schwarzschild mode 
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there corresponds a partner mode 

= 751 ((^'""""(^ - - (3,99) 

where we have set for simphcity the constants which appear in eq. (13.371) to 
K = 4M, B = 1, K = 1. The Unruh modes which are positive frequency on 
are given by 

<^A = a^^T' + PxvT'^* A>0 
<^^A = PxipT'* + ax^T'^ A>0 (3.100) 

To these modes are associated the operators a^"*, a^"*^, ax- 

The Heisenberg state |0m) is that anihilated by the dx operators. It can 
be expressed in terms of out states as 

io-) = ^ne"^ ' " Kut) (3.101) 

thereby exhibiting the correlations between the produced Hawking quanta 
(a^"*^) and the partners (a^"*^^). Following the development of eqs (12.1021) 
et seq. one shows that to each produced hawking photon in a wave packet 

^lj,= dX -i,xvT (3-102) 

JO 

there corresponds a partner in the packet N^^ Jo°dX 'jix (Px/o^x) v^a"*^ where 
Ni is a normalization factor given by = J dX\'jix{l3x/<yx)\'^- Note that the 
wave function of the partner does not have the same mode decomposition as 

We can now compute the weak value of the energy momentum tensor, i.e. 
the energy correlated to the creation of a photon in mode ijji. 

^ ^^^'^ " (o.„|n|oj ^^-'"^^ 



As in Section [231 the projector 11 selects the final state wherein one photon 
is produced in the mode ipi. 

U = Ipartners ® dX-fixal\Oout){Oout\ d^l*xO'\ (3.104) 

Jo Jo 
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A calculation similar to that leading to eq. fl2.105p yields for Tyy = dyCpdyCp 
(a similar expression obtains for T^u) 

vv/w ^ 



(Oouti -Tyy. \Oin) 



(Oin I -^Dt) I Ojn)ren (3.105) 



The second and third terms of eq. fl3.105p are background (they are indepen- 
dent of ipi)- The third term was the subject of Section [373] and the second is 
the difference between Unruh and Boulware vacuum. Therefore in computing 
these weak values the background is that of Boulware vacuum. This is the 
precise analogue of Section 12.51 where the background was Rindler vacuum. 

The first term (hereafter referred to as (T^i/)^J describes the energy mo- 
mentum of the vacuum fluctuation which will become the Hawking photon 
ipi and its partner. From eq. (13.1051) it is apparent that the energy of this 
vacuum fluctuation vanishes. Indeed the annihilation of the vacuum by the 
total energy operator dv -.Tyy : |0i„) = implies that the integral of the 
flrst two terms on the r.h.s. of eq. f l3.105p vanish as in eq. f l2.116p . 

It is convenient to rewrite (Tyy)^^ in terms of out modes to obtain 

iJ.^dX ^a{( 3l/al)dycpr) ilo^^dX ' il^dy^T^ 
S^dX H^mjal) 



{J-vv}i,i = 2&{Vh-V) roo^x I., 19/^2/ 2\ + 



2e{v-VH)- 



j^dx IWiPl/al) 

(3.106) 



In this form it is apparent that the 6{v — vh) piece is real and positive whereas 
the 9{vh — v) piece is not. It is complex and oscillates in such a way that 
the total energy dvTyy vanishes. 

Only the piece 6{vh — v) is reflected in flnite Schwarzschild time and 
reaches X^. There it takes the form 

^ ""^'"^ - Io^dX\U'{f3l/al) ^'-'"'^ 

It carries the Schwarzschild energy Aq of the post selected photon. Indeed 
one has 

du{T^u)^, = ^TSSTTH i^r^/ — ^ - Ao 3.108 

Jo "A I'jixrlfJx/axr 
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where Aq is the energy of the produced hawking photon. 

We mention that one may in similar fashion postselect the presence of any 
number Hawking photons. In this case one will find a decomposition of the 
weak value very similar to eq. fl3.105p . There will be a term corresponding 
to the Boulware energy (the second line of eq. (13.1051) ) and a sum of terms 
of the form {T^^)^., one for each post selected photon. Thus the background 
contribution can be naturally separated from the fluctuating contributions. 

As in the electric field (see eq. (I1.59P ). the physical relevance of the 
imaginary part of (Tuu)^, can be seen by modifying slightly the background 
geometry: g^y = g^^ + 5g^u- Then the change in the probability of finding 
on X+ the Hawking photon selected by H (eq. fl3.104p ) is given by 

Pg,.+&9,. = Pg,. (l - / d'x^6g,,2lm [(T^.)^]) (3.109) 

3.6 Thermodynamics of Black Holes 

In order not to disrupt the continuity of the text we shall introduce the no- 
tion of black hole entropy by giving a thermodynamic interpretation to the 
evaporation phenomenon which has been the subject of this review until this 
point. The argument pursued will be heuristic and based on analogy. Sub- 
sequently a more rigorous derivation based on a true equilibrium situation, 
the eternal black hole will be presented. 

A convenient pedagogical crutch to start with is the idealization used in 
the preceding chapters, the 1+1 dimensional problem with unit transmission 
coefficient. We emphasize at the outset that the notion (and the value) of 
the entropy of the black hole itself which is deduced from this model, is by 
no means contingent on the idealization used to get it. 

The central result of the idealized model is that at a certain time a ther- 
mal flux emerges as the collapsing star approaches its horizon. Then the 
mean energy emitted during Av (or At at fixed r) is given in eq. (I3.46P : 
< Tuu > At = j^PjfAt and it is distributed thermally (with = SnM. 
The probability for the simultaneous occurrence of quanta in a state whose 
occupation numbers for frequency Ui are rii is 

P|„j = Z-ie-^^S"''^^ (3.110) 

We remind the reader that the density matrix eq. fl3.110p is a consequence 
of the Bogoljubov transformation eq. (13.391) . The normalizing factor Z is 
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likened to a partition function of a one dimensional gas (c.f. eq. fl2.45l) ) 
whose volume is At, i.e. 



InZ = - 




(3.111) 



and the mean number of quanta with frequency uji is 

in.) = ^ = (e^-' - 1)-^ (3.112) 

The fact that the photons are all outgoing does no injustice to the application 
of usual statistical mechanics and thermodynamics defined by the canonical 
ensemble. 

The conditions under which evaporation is taking place are isothermal to 
a very good approximation. By that we mean there is a At sufficiently small 
so that even when a macroscopic number of photons evaporate in this time 
interval, jSn does not change. For example for M ~ 1 solar mass (= 10^^ 
Planck masses) during the time interval that 10^° photons are emitted one 
has AM = O{10^yM) and APh/Ph = 0{AM/M) = O{10-^^). 

One may liken this quantum evaporation to the irreversible process in 
which a large mass of liquid evaporates into a tiny amount of vapor in time 
At. The whole system is insulated from its surroundings (i.e. Efotai = 
constant) with a little vapor occupying a space above the liquid. Then 
increase the volume a little bit. To make the illustration more cogent we 
neglect terms proportional to the chemical potential (after all it doesn't cost 
anything to make a vacuum fluctuation). Then the only increase in entropy 
is due to the increase of volume of the vapor phase. The liquid loses a "lit- 
tle" energy AE and the vapor increases in energy AE where "little" means 
AE / Etotai << 1 and we envision that almost all the energy is localized in 
the liquid. Then the temperature change is negligible {AT /T = 0{AE/E)). 
One says that the liquid is a "reservoir". The "reservoir limit" for the change 
in its entropy is AS* = AE/T. [In thermodynamics one introduces a reservoir 
to convert entropy considerations of a system to its free energy, the relevant 
function to describe isothermal processes undergone by a system in contact 
with the reservoir]. 

From the above interpretation we are led to ascribe to the black hole 
an entropy. In point of fact it was this remarkable insight of Bekenstein 
[8] that must have incited Hawking to attribute a temperature to the black 
hole and hence evaporation. During the time At wherein a mass AMbh 
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(= -/\Mgas = -J2i{ni)uJi) is evaporated ( with \AMbh\/M « 1), the 
process is effectively isothermal and the black hole acts as a reservoir. Its 
change in entropy is thus 

ASbh = Ph^Mbh = StiMAMbh ^ A(yl/4) (3.113) 

where we have introduced the area of the horizon surface {A = 47r(2M)^). 

A noteworthy feature of black hole evaporation as compared to the liquid 
vapor analogy is that no time is required to get the "vapor" into an equilib- 
rium state after it evaporates. To complete this discussion, one may calculate 
the increase in entropy occasioned by the evaporation. Using eq. (13.1101) one 
has 

AStotal = ASBH + ASgas=/3H{AMBH+AMgas) + ^^At 

- ^A<^pA* ,3.114) 

Recall that the pressure is given by din Z/dV and in our case AV is the one 
dimensional volume At. Thus the entropy increase is pAV as for the analog 
liquid- vapor example (with zero chemical potential). 

But how is that? We started with one state and we are now calculating 
the increase in the total number of states. The answer of course is that we 
are describing the gas by a density matrix and have completely forgotten 
the correlation to the degrees of freedom left within the black hole. It is 
therefore the extra trace over these degrees of freedom which is responsible 
for the increase of entropy AStotai and this is a perfectly legitimate description 
for the outside observer. 

It is thus possible to develop a phenomenological thermodynamics for the 
outside observer which ascribes to the black hole an entropy change which 
is (3h{M)AM. We may then calculate the entropy of the black hole by 
integrating this change from to M as to obtain 

i-M 

Sbh{M) - Sbh{0) = (3H{M')dM' = 'n{2Mf = A/A (3.115) 

Jo 

One may ask why Sbh{M) is intrinsic to the black hole i.e. equal to logQ{M), 
where Q is the number of degrees of freedom of the black hole. After all the 
black hole was used as a reservoir only. The answer is that the black hole 
even if it is a reservoir for the radiation is at the same time in very good 
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approximation a microcanonical ensemble unto itself. This is because M 
changes so slowly {dM/dv ~ — 1/M^). Thus /3h = dMlogVt{M) is a valid 
microcanonical expression for how the density of states varies with the mass. 
In usual statistical mechanics, one turns the argument the other way around 
since one starts with the density of states. The (temperature)"^ is then 
introduced as the logarithmic derivative of VL. 

A possible constant of integration in eq. (13.1151) is a subject of much 
debate [23j. If there is a remnant at the end of evaporation then Sbh{0) is 
its entropy. Since we have as yet no way to think about this remnant it is the 
better part of valor not to commit oneself as to its value. If the remnant leaks 
away then this leakage should be accompanied by further entropy increase, in 
a model dependent way and in the sense of a density matrix. The situation is 
analogous to that sketched in Fig. (12.41) for a decelerating mirror. In principle 
one could measure the correlations between these late leakage photons and 
the earlier radiated ones. Then one would be able to check whether or not 
there is a pure state with no change in entropy. Thus the former entropy 
increase would have resulted from neglecting the correlations. If otherwise 
the remnant sinks through a singularity these degrees of freedom get lost 
and the entropy increase is given by eq. (13.1141) . i.e. we are stuck with the 
density matrix description and quantum mechanics applied to this problem 
has become non unitarity as initially suggested by Hawking [16]. Nothing 
could be more interesting-or exasperating. For further discussion, see Section 
1X71 

We mention that when one takes into account the transmission coefficient, 
much of what has been said still survives. Eq. (I3.110p must be changed 
in that each factor e"^*^' must be multiplied by a transmission factor F,, 
and one must include all the angular momentum modes. The value of Z 
changes, but thermodynamics is retained in a modified sense since for N » 
1 (where = J2ini), relative energy fiuctuations are still 0{1/^/N) because 
the distribution remains Poisson, and the process is still isothermal. All 
that is modified is the numerical value of AStotai- In particular eq. (13.1141) 
is retained since —d\iaZ/dj3\/3^ = AMgas for the modified ensemble as well. 
The entropy ascribed to the black hole and its physical interpretation as 
described after eq. (13.1141) . is independent of the mechanical details. And 
one still has AStotai = pAt > 0. 

The idea of black hole entropy in terms of a reservoir is reinforced by 
the analysis of the equilibrium situation which is characteristic of the eternal 
black hole. This idea is once more a seminal discovery of Hawking [IS]. We 
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shall follow Hawking's idea, but, unlike him, we shall also deduce the black 
hole entropy itself. 

In a closed box of volume V, it has been shown by Hawking that for a 
sufficiently large energy, the energy becomes partitioned into a black hole 
surrounded by radiation in thermal equilibrium. We shall see below that 
there is a limit wherein essentially all the energy is in the black hole whose 
radius is nevertheless much smaller than the radius of the box. Then, the only 
role of the radiation is to furnish a temperature. In the limit envisaged, this 
latter is the Hawking temperature, Th, which value is then used to compute 
the entropy of the hole. 

The relevant geometry which describes the static situation characteristic 
of the eternal black hole is the full Schwarzschild quadrant R (r > 2M) of 
Kruskal space. This geometry is depicted in Fig. (13.11) . In the case where 
almost all the mass is in the black hole the geometry in R approximates to 
that of empty Schwarzschild space. This space has both past and future 
horizons. For the case envisioned of a black hole formed from a collapsing 
shell the past horizon is not present (see Fig. (13.21) ). Nevertheless, the use 
of the eternal black hole geometry is a correct mathematical idealization to 
be understood as follows. Outgoing photons that issue from the black hole, 
in the course of its formation at late stages are Kruskal in character (i.e. the 
modes used to describe Unruh vacuum). Because the black hole is enclosed 
in a finite volume, these modes get refiected off the walls and come back 
as Kruskal modes. So after some time a stationary state gets established 
in which both the incoming and outgoing modes (or pieces thereof) are of 
Kruskal character. Note that all angular momenta participate in thermal 
equilibrium since considerations of reversibility render useless any appeal to 
the smallness of the transmission coefficient for the higher angular momentum 
waves. In this equilibrium case, the hole is surrounded by a gas of energy 
density proportional to T^. So the s-wave truncation makes no sense in this 
case. The state is then set up in terms of the modes which traverse the whole 
space. But the part of the space which is physically relevant is restricted to 
that part of R which is bordered by the surface of infalling matter. 

Every stationary state of the photon gas in the eternal geometry except 
the Hartle-Hawking vacuum gives rise to a singular energy momentum tensor 
on at least one of the horizons. The Hartle-Hawking vacuum |33] is that 
constructed from the quanta of Kruskal modes. It is the analog of Minkowski 
vacuum in Minkowski space and R is the analog of the right Rindler quadrant. 
One can understand the origin of the theorem on singular energy density 
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near the horizon by reference to eq. fl3.64l) et seq. where it was shown that 
in Boulware vacuum < T(/^ > is singular. To undo this singularity clearly 
requires a very special condition (which, for the collapsing case, was shown 
to be {Tuuir ~ 2M)) = 0[(r - 2Mf]). 

In order to prove this point, to see the thermal character of Hartle- 
Hawking vacuum in R, and in fact to construct the corresponding density 
matrix, it is useful to make the euclidian continuation of the Schwarzschild 
geometry t — >■ it (a construction which is possible owing to the staticity of 
the geometry). Forgetting angles the line element is (1 — 2M / r){d{it)Y + 
(1 — 2M/r)~^dr'^. This has the form of a cigar which terminates at r = 2M 
and extends as a cylinder of radius 4M out to infinity. To see this note that 
near r = 2M, the line element takes the form 

ds^ p^dO^ + dp^ (3.116) 

where p = 2^2M{r - 2M) ; 9 = it /AM. Equation fl3.11fip is the line element 
in the neighborhood of r = 2M written in local euclidean polar coordinates. 
No conical singularity at that point requires periodicity of 9 equal to 2n. On 
the other hand as r — > oo the metric goes over to the metric of a cylinder 
(= {AMyd9'^ + dr"^), the periodicity in 9 remains equal to 27r. 

Thus regular functions defined on this space are periodic in it with period 
equal to 2tt{4M) = (3h- For large r, Schwarzschild t coincides with proper 
time, hence is temperature at large r. And Green's functions defined as 
Hartle-Hawking expectation values are regular. A discussion of these Green's 
functions, in general, is given in ref. [37]. For the purpose at hand it suffices 
to note that near the horizon the geometry is regular (eq. (13.1161) ). In this 
region Hartle-Hawking Green's functions are related to Boulware Green's 
functions in the same way that in flat space Minkowski Green's functions are 
related to Rindler Green's functions. Near the horizon the analogy is exact 
with the acceleration replaced by p~^. 

Having displayed the properties at equilibrium, let us now construct the 
necessary conditions to derive the black hole entropy. Our condition of neg- 
ligible energy in the gas is 

M » VM-^ thus M » V^l^ (3.117) 

On the other hand we require that the volume occupied by the gas be much 
greater than that occupied by the hole in order to validate the estimate of 
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its energy (i.e. Egas = VT"^ as in absence of gravity) 

V » (3.118) 

Equations 03.1 17p and 03.1181) are compatible if >> i.e. M » 1 . 
In that case the total energy at equilibrium is in very good approximation 

E = M + V(3h'^ (3.119) 

We now use the fact that this equilibrium configuration should be derivable 
from the variational principle of entropy SStotai\Ey = 0. Indeed, by attribut- 
ing an entropy to the black hole (Sbh) and by taking the variation with the 
energy repartition {SM = —5 Egas) , one finds 

= '^^(^-/5h) = (3.120) 

whence the equality of the temperatures gives 

dSsH 



dM 



SnM (3.121) 



In eq. ( 13.1201] . we have used conventional canonical thermodynamics for the 
radiation (i.e. dEgas = TudSgas) and ascribed to it the temperature albeit 
that the total system is microcanonical in which one phase (the black hole) 
acts as a reservoir. The above considerations are quite general and the use 
of the photon gas was by way of illustration. Equation O3.120p will be true 
for any model of matter. In this sense eq. 03.12ip is a very powerful result. 
Whatever happens in the ultimate destiny of black hole physics one would 
be loath to give it up. 

A critique (see ref. [76]) of the notion of black hole entropy follows from a 
closer inspection of the distribution of energy density in the Hartle-Hawking 
vacuum [35] . Vacuum polarization effects in < T"^ > ( which lead to negative 
energy density near the horizon) prevent a clean split of total entropy into 
its matter and black hole components. Therefore the value Sbh = ^4/4 is to 
be considered as valid only in the reservoir limit as discussed above. This 
does not mean that Stotai does not exist but rather that in the general case it 
cannot indeed be unambiguously divided into two terms (matter and black 
hole) because gravitation is a long range force. 
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That the entropy is proportional to the area of the horizon has been based 
on the fact that the Hawking temperature is proportional to M~^. This in 
turn is related to the imaginary period of the relevant Green's functions of 
quantum fields in Schwarzschild space. A transfer of dimensions has occurred 
from mass~^ to length (or time since c = 1) owing to the existence of h. 
Thus it would seem that the identification of entropy with area is essentially 
quantum in character. But quite surprisingly, a purely classical development 
already foretells a good bit. On one hand, Hawking and collaborators have 
shown that in the classical evolution of matter-gravitation configurations, 
the area of the horizon always increases. [This review is not the place to 
prove this important classical topological theorem. The reader is referred to 
the important monograph of Ellis and Hawking [Sj.] This can be related 
to more general results that the total entropy of black holes and matter 
always increases (see ref. jp]). On the other hand, Bardeen, Carter and 
Hawking [7j have used a tool called the Killing identity to investigate the 
role of the horizon's area as an entropy, in purely classical terms. The point 
of departure is an exact geometrical identity satisfied by Killing vectors, upon 
which one grafts Einstein's equations to relate the curvature, which appears 
in the identity, to the energy momentum tensor. The version we shall present 
below is applicable to static spherically symmetric systems (See also ref. [23] 
for a similar derivation in the framework of the hamiltonian formalism). In 
addition we refer to efforts [SH] , [IB] which identify the black hole entropy as 
an action integral for the pure gravitational sector. The value Sbh = A/4 is 
then obtained from the classical Einstein action taking due care of boundary 
terms. The bearing of this result on the characterization of the degrees of 
freedom locked within the black hole (but coupled to the external world) 
remains a subject of debate. 

Spherically symmetric static systems can be described by the line element 

ds^ = -e^^^^^de + e2^('^)dr2 + r^dfi^ (3122) 

We assume that at large r there is no matter present, so Birkhoff's theo- 
rem applies and one has 

e-2A =1 r ^ 00 3.123 

r 

Further we fix completely the coordinates by imposing (/)(r = 00) = 0. M 
is the total Keplerian mass measured from infinity. It is the sum of the 
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black hole and matter mass. Indeed, the equation Rq — R/2 = SttTq can be 
integrated exactly to yield 

g-2A = (1 — 

m(r) = M- drAnr^pir) (3.124) 

The horizon of the black hole is the radius r = 2m b at which e~^^ van- 
ishes. Regularity of the geometry then implies that in the vicinity of the 
horizon one has 

r,« r — 2mn 9, ,r,r — 2mn 
g-2A _ B . _ ^2 £ (3.125) 

2mB 2mB 

The value of k differs from unity owing to the presence of matter. It plays 
an important conceptual role in what follows. 

The relevant geometrical (Killing) identity in this approach is 

i?OrV+^ = -(e'^-VV)' , (3.126) 



a consequence of the staticity of the geometry . One then uses Einstein's 
equation Rq = S7i{Tq — |T). We shall suppose that the matter is a perfect 
fluid = (p + p)^o^'ji -~ P^'j[ (introducing a difference between radial and 
tangential pressure does not change the final answer), to give 

i.(e<^-A^2^) = Arcip + 3p)e'^+^ (3.127) 
ar ar 

whereupon integration from r = 2m b to r = 00 yields 

M-kmB = ATx (p + 3p)rV+^dr (3.128) 

J2mB 

(This is an alternative way of writing M as compared to M = mB+An pr^dr; 
see eq. (Kim .) 

The term kmB will be expressed henceforward as a term proportional to 
the black hole area {A = 47r(2mB)^) 

kmB = -^A ; K = k/AmB (3.129) 

477 

The constant k, called the surface gravity, has very important physical sig- 
nificance. It is the gravitational acceleration at radius r measured at infinity. 
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for example from the tension in a string attached to a test mass located at r 
[70] . To see this, suppose the test mass, /i, initially at rest, is dropped by a 
static observer at r. It picks up kinetic energy T (as measured by a static ob- 
server at r + Sr) in an infinitesimal distance 6r (ie. proper distance 61 = e^Sr) 
equal to fi{d(j)/dr)Sr. When measured at r = cxd, this energy is redshifted to 
fie'^(f)'6r. Equating it to a work term F^6l one finds = fie'^'^cp'. For r 
at the horizon this gives F^/ ^ = k. It is noteworthy that in the absence of 
matter (with e^'^ = e~^^ = (1 — 2ms/r)) one has F^o = ixmB/r"^- Consistency 
with eq. fl3.127p when p = p = is an essential point. Were Newton's law 
for the force other than r~^, the repercussions would be serious indeed. 

The reason why we have belabored the physical interpretation of k is 
that in the euclidean continuation of the metric near the horizon, eq. (13.1251) 
gives rise to the polar coordinate representation of fiat space which for 
Schwarzschild space is given by eq. (13.1161) with 6 = it /{AM). The only 
difference when k ^ 1 is 6 = k{it) / {AniB) = K{it). When transcribed into 
quantum mechanics this period in proper euclidean time is transformed into 
the inverse Hawking temperature (3h = STfrriB/k = 27r/K. 

To see that the area has to do with an entropy one compares two static 
solutions by varying external parameters p, p, tub, etc in such manner as to 
be consistent with Einstein's equations. The steps require a bit of algebra 
and is relegated to a parenthesis. We first quote the result 



6M = —6A + / {fd6S.matter + flddKmatter) + PBkA5{2mB) (3. 130) 



Here T and p are the local temperature and chemical potential, scaled cor- 
rectly by the red shift (Tolman scaling) T = e'^T{r), jl = e'^p{r). SSmatter 
and SNmatter are the local entropy and particle number of matter. There are 
no pressure terms at infinity since the space is asymptotically flat but we 
have retained the work term due to the pressure pb occurring near the hori- 
zon. This term doesn't appear in the original formula of Bardeen, Carter, 
Hawking who assumed that p and p vanish at the horizon. However it be- 
comes relevant for instance when we compare geometries of black holes in 
a box surrounded by radiation in a Hartle-Hawking state. If instead of the 
model considered here above one considers a black hole and surrounding 
matter enclosed in a finite volume having walls, the energy of neighboring 
configurations have to differ by a term like —p6V but also by a term (model 
dependent) taking into account the stress in the wall. 
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[The difference between two neighboring solutions are characterized by 
5M, 6mB, Sk, 6p, 6p, 5(j) and 5A, quantities that cannot be all independent 
because of the Einstein field equations. Starting from eq. fl3.128p we obtain 

SM = —S{kA) + 47r(5 / (p + 3pye^+^dr (3.131) 

4:71 J2ms 

The variation of this last integral is the main task of the calculation. This is 
most easy done by splitting it as a sum of two parts using —R = Sn^—p + Sp) 



471 {p + 3pye'^+^dr = / -fle^+\^dr 

POO 

+87r / pe'^^^'^dr (3.132) 

J 2mB 

The first term, hereafter denoted /i, on the right hand side is the gravitational 
action. Its variation must take into account boundary terms which arise both 
because tub varies and by integration by parts to yield 

51^ = e'^+^r^n2mJmB + ( lim - lim )e'''-\\6(j)' + (j)'6(j) - (- + (p')SA) 

POO 

+ (G°50 + G';M)e<^+Vcir (3.133) 

J 2mg 

To prepare the evaluation of the variation of the second term, let us recall 
the thermodynamic relation 

6E = -pSV + TSS + pSN (3. 134) 

which in terms of local densities on 3-surfaces t = constant becomes 

6{p4:7!-e^r^dr) = -p6{4:7ie^r'^dr) + T6{s4:7ce^r^dr) + p5{n47ie^r'^dr) 

= -pA7Te^r^5Adr + Td6S + pd5N (3.135) 

Acordingly we obtain 

5l2 = -87rpe'^^''r%mJi2mB) 

POO 

+2 / e'^+^47rpr^6^dr 

POO r- -| 

+2/ e'^ Td6S + pd6N) - p4nr^e^6Adr (3.136) 
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The variation of kA results from the variation of both and and A: 

+e<^-\2(50' + - 6A)hms (3.137) 

Putting all together and using the Einstein equations Gq = 8np, = Snp 
and the expression of i?/2 = l/r^-e-^^{^+(j)"+(l)'^-(f)'A'+2(p'/r-2A'/r) and 
the asymptotic behavior of the metric components e*^^^ = 1, Scj)' = SM/r"^, 
5(f) = —6A = —6M/r and the condition defining the radial coordinate at the 
horizon e~'^^{6A + A' 6(27713)) = we obtain 

roo _ 

6M = kdmB+ {Td6S + fid6N) -4ne^+^pr^\2mJ{2mBX3.138) 

J2mg 

where T = e'^T and fl = e'^fi are the Tolman temperature and chemical 
potential. On the horizon, staticity implies Rqo = 0, ie. p + p = (the same 
can be deduced from the regularity of \/—c^2mB = kAm^), so the last term 
can be translated into a term giving the work (measured at infinity) due to 
the pressure at the horizon when the volume of the exterior of the black hole 
varies.] 

In eq. fl3.130p it is difficult to resist setting the black hole entropy equal 
to A/4 since k/2'k is the periodicity of the euclidianized time hence the tem- 
perature. How is it that the classical theory anticipates quantum mechanics? 
To this end, we cite the efforts of Gibbons and Hawking and Brown and York 
referred to above, based as they are on the classical euclidean gravitational 
action with imaginary time of period (3. According to the tenets of quantum 
statistical theory, one converts action to entropy universally through divi- 
sion by h. It is only through quantum mechanics that entropy acquires an 
absolute sense. The term ndA in eq. fl3.130p is derivable from the difference 
between the classical actions of gravity which arises owing to different mat- 
ter configurations. From the above identification of h/2t{ with Tjj/h, one 
rewrites KdA/Sn as TndSBH hence dSsH = dA/ATi. Thus it seems that the 
classical Killing indentity "derivation" of entropy has no thermodynamical 
interpretation without quantum mechanics (See ref. [33]). 
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3.7 Problems and Perspectives 



As elegant as is the semi-classical theory of black hole radiation, it is fraught 
with severe conceptual problems. 

The most essential is concerned with the consequences of the exponential 
increase of the energy densities of the vacuum fluctuations inside the star 
and about the horizon which are converted into Hawking photons [51], [53] . 
This is due to the Doppler shift, eq. (13.431) 

00 = Xe^l^^ (3.139) 

as the surface of the star approaches the horizon {u oo). Therefore the 
Planck scale uj = 0(1) is reached very early in the history of the evaporation 
for a characteristic value of A of 0{M~^)) after a time u = O(MlnM). This 
should be compared to the lifetime of the black hole 0{M^). Thus after the 
emission of a few photons {An = O(lnM)), Hawking radiation is concerned 
with the conversion of "transplanckian" vacuum fluctuations (of frequencies 
u) into "cisplanckian" photons (of frequencies A). The point to emphasize 
is the mixing of radically different energy scales which is occasioned by the 
exponential growth of the Doppler shift. This is in contradistinction to more 
conventional physics wherein different scales remain separate (e.g. atomic 
versus nuclear). 

We recall that the observation of a photon near X"*" implies the existence 
of a particular localized vacuum fluctuation all the way back on X^. The 
energies involved in this vacuum fluctuation are of the order of uo. Thus there 
is an implicit assumption of inflnite mean free path of modes. This is implied 
ab initio through the use of free fleld theory (00 = 0). Whilst there is nothing 
wrong with that insofar as elementary particle interactions are concerned 
(their scale being of 0(Gev or Tev) and presumably being asymptotically 
free), it is most probably incorrect since the gravitational interactions really 
get in the way. In fact, if anything, the existence of the Planckian mass 
scale would inevitably tend to stronger forces at higher energies (in local 
fleld theory at least, but maybe not for string theory). Be that as it may, 
if one uses the Newtonian law of gravitation as a guide, two spherically 
symmetric shells of mass m at radius r will have a gravitational interaction 
energy which exceeds the mass m, for m = 0(1) and r = 0(1). From the 
above considerations, it is seen that these scales are attained at the threshold 
of entry into the transplanckian region {u = 0(1), u = O(MlnM)) as the 
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fluctuation approaches the center of the star. The assumption of free field 
theory is thus a priori completely inadequate. 

Whatever is the accommodation of the transplanckian fluctuations within 
the star, there is another vexing short distance problem which comes about 
from the fluctuations of the position of the apparent horizon due to the 
emission or non emission of a Hawking photon. In the semi-classical theory 
the change in radius of the apparent horizon is due to {T^v)ren only, since 
{Tuu)ren = ou the horlzou. However the post-selected presence (absence) 
of a photon around uq is contiguous with the contraction (expansion) of the 
apparent horizon with respect to the mean evolution. More precisely, when 
the Hawking photon of frequency A = (measured on X+) is within 
of the horizon (i.e. r — 2M at fixed v is less than M~^) the classical geometric 
concept of the apparent horizon (the locus where dvr\u = 0) loses meaning 
since its radius fluctuates by 0{M~^). So the approximation of a free field 
in a fixed background seems to break down near the horizon. Thus how can 
we be sure that the result of the semiclassical theory wherein there is no 
coupling between u and v modes corresponds to the true physics near the 
horizon ? Might there not arise crucial correlations which are absent in the 
semi-classical theory {{TuuTm) — {Tuu){Tvv) 7^ 0) ? Then putting the whole 
blame on (T„„) for the change of the geometry near the horizon could be 
quite misleading. The mechanism of evaporation might then be closer to 
pair creation near the horizon. 

Does this mean that everything contained in the semi classical approxi- 
mation is irrelevant j54j ? Probably not. For one thing it is very unlikely that 
Hawking radiation does not occur. One must distinguish between the theory 
of vacuum fluctuations based on free field theory (given in Section 13.51) and 
Hawking radiation as a mean theory (see Section [?!^ . In the first, expecta- 
tion values of (T^^,) in the in-vacuum wash out the fluctuations and provide 
a dynamical origin of the Hawking flux. So one can refer to the derivation 
of Hawking radiation in a strong or weak sense, i.e. as derived from free 
field theory or as some effective theory giving rise to a similar (T^u)- In 
support of the existence of such a effective theory, one needs but appeal to 
the regularity of (T^^) on the horizon(s). Indeed both in the collapsing and 
eternal situations (Sections 13.31 and 13.61) it has been emphasized that one 
can envisage Hawking radiation as a response to an incipient singularity of 
the Boulware vacuum at the horizon(s) in such a manner as to erase that 
singularity. (To see the connection between the eternal and collapsing case 
imagine punching a small hole in the surface of a recipient containing the 
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"eternal" black hole. It then becomes "ephemeral". Radiation leaks out 
and this is neither more nor less than Hawking radiation albeit with a small 
transmission coefficient. In fact the potential barrier which stops s- waves 
with energy smaller than 0(1/M) and almost all higher angular momen- 
tum modes plays already the role of the small hole since one has almost a 
thermal equilibrium behind the barrier [20]). Such a general consideration 
might well be inherent in the complete theory wherein one accounts for the 
gravitational quantum back reaction. This latter could give rise to violent 
fluctuations at Planckian distances which nevertheless leave a regular mean 
to drive a classical background geometry at larger scales. 

How might we then envisage the fluctuations ? Almost certainly one 
would expect that the reduction of the problem to free s-wave modes is in- 
correct. Rather modes will interact mixing angular momenta out to distances 
where the notion of a free field is legitimate -sufficiently outside the horizon 
so that the Doppler shift does not entail transplanckian frequencies when 
the Hawking photon is extrapolated backwards in time. Within the interior 
region one would expect a "soup". And it is to be noted that this soup is 
present in Minkowski space, since it is inside the star. 

It is interesting to speculate on how one might describe this situation. 
After all, even if there is an ungainly soup, it nevertheless fluctuates in a 
systematic way if one imposes spherical symmetry and translational symme- 
try in time. The fluctuations can then be sorted out according to angular 
momentum and frequency; perhaps they must be endowed with a lifetime 
as well. These things are quasi particles. If Hawking radiation exists then 
we know for sure that as a quasi particle passes through the horizon region 
it gets converted from "quasi" to the free field fluctuations which we have 
treated in this review. So it should be possible to come to grips with this 
problem somewhere in the middle region. It is not impossible that one will 
be able to prove that the Hawking radiation develops out of this Planckian 
nether nether land and that the true transplanckian fluctuations are irrele- 
vant After all these latter don't seem to bother us very much here and 
now so there is some hope that they are no nuisance there either. They may 
even result, as suggested above, in a Planckian spread of the region between 
the apparent and real horizon since they incorporate gravity within them. 

■^In a recent numerical calculation in a model which is analogous to the black hole, 
Unruh [?] has shown that a severe modification of the dispersion relation uj{k) for k 
greater than some threshold value in no way affects the thermal spectrum of Hawking 
emission. 
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Another very disturbing problem arises in the semi classical theory. This 
is the so called unitarity issue [IH]. In usual evaporation of an isolated 
system into vacuum, correlations are always present. At the early stages 
an evaporated molecule gives a kick back to the unevaporated mass caus- 
ing correlations between what has and what has not evaporated. At later 
stages these correlations get transferred to correlations among the evaporated 
molecules [OS]- If one were in a pure quantum state, it is these correlations 
which encode its purity. Of course one says evaporation is accompanied by 
entropy increase, hence increase in the number of states. But we attribute 
this to the coarse graining that is implicit in the definition of entropy. 

In black hole evaporation, in the semi classical theory, the correlations are 
in nature similar to that encountered in usual evaporation in its early stages. 
Reference to Section 13.51 shows that each Hawking photon leaves behind a 
"partner" , a field configuration in vacuum of a specified local character within 
the star. As the evaporation proceeds these correlating configurations build 
up. However in the semi classical theory these configurations never get out 
(they are shut up in the closed geometry which devellops during the evap- 
oration see Section 13.41) and the Hawking radiation contains no information 
on the quantum state of the star. One does not recover unitarity at r = oo. 

A few options (see for instance the review article of Preskill [78] ) seem 
available on how to confront this situation: 

• As suggested by the semi classical theory, these inside configurations are 
forever lost to the outside observer (they could end up in the singularity, 
or in an infinitely long lived remnant). Unitarity is truly violated as 
originally claimed Hawking |46j. 

• The semi classical theory fails at Planckian size black hole, evapora- 
tion stops and a finite long living remnant forms whereupon the cor- 
relations between the Hawking quanta, their partners and the star's 
matter are recovered at r = cxd to reconstruct purity. We recall that for 
the accelerating mirror (Section 12.51) the correlations to the partners 
are completely recovered upon decelerating the mirror. However these 
considerations cannot be applied directly for the black hole since, for 
purity, the degrees of freedom of the star should be recovered as well. 

• The semi-classical picture is all wrong in this regard and the correla- 
tions occur outside the horizon. All the information about the state of 
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the star leaks out in the radiation [JT] [52] [SI]- This option however ne- 
cessitates either a violation of causality or a fundamental revision of the 
concept of background geometry at scales large compared to the Planck 
scale. Some authors entertain the thought that different backgrounds 
are appropriate for different observers |89j. e.g. Schwarzschildian or free 
falling observer. The backgrounds would be "post-selected" by the ob- 
server. 

In the opinion of the authors it is futile to confront the unitarity issue 
without some clear ideas about the transplanckian issue i.e. quantum gravity. 
Indeed the unitary problem cannot be settled without a deep understanding 
of the dynamical origin of black hole radiation. 

We also wish to point out that in this quest the fact that the entropy of 
a black hole is proportional to its area may play a vital role. How is it that 
the entropy of a black hole is equal to the number of Planckian cells that 
are necessary to pave its surface? The horizon seems to block out the cells 
which lie deeper than a Planck length within the hole. Is this related to the 
expected scenario that emission will occur at the surface outside the apparent 
horizon i.e. where a quantum fluctuation begins to belie its presence? 

Such are the problems that one must face. Whether their solution will 
lead to the quantum theory of gravity or the inverse is a moot point. And 
this primer is certainly not the place to speculate any further on the question. 

No doubt there will turn up further stormy weather to stir up the already 
troubled waters that must be traversed on this journey to terra incognita. 
Nevertheless we wish the reader at least some fair weather. Good Luck and 
Bon Voyage. 
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Appendix A 

Bogoljubov Transformation. 



In superfluid helium at rest a macroscopic number of particles occupy the 
zero momentum state < agao >= and Nq/N is a finite fraction as well as 
N/V where V is the volume, the total number. The thermodynamic limit 
is — > oo, \^ — > oo with Nq/N and N/V fixed. The commutation relation 
[ao,ao"] = 1 is then a negligible consideration when considering operators 
containing oq and Oq as products that multiply the typical unperturbed states 
which make up the vacuum (ground state), these unperturbed states being 
eigenfunctions of < >. 

The hamiltonian of interacting particles is (with C,k = k'^/'2m) 

H = Ho + V 

k 

^ = v{ki,k2;k3,k4) dk^+k^M+kiOk^Ok^ak^ak^ 

kiMMM 

(A.l) 

The interaction potential has matrix elements of 0(1/V"^) since the un- 
perturbed states in a box are e*^**/ VV. 

Each J2k is of 0{V) so each of the terms in eq. flA.ip is 0{N), /x is a 
chemical potential put in for convenience so that one can allow A^ to fiuctuate 
albeit such that < AA^^ > / <N >^= 0{1/N). 

The unperturbed ground state has Nq = N, fi = 0, and E = 0. The idea 
of Bogoljubov [121 was to develop a perturbation theory in the small number 
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[{Nq/N) — 1]. So the technique is to keep terms in leading order in -y/iVo in eq. 
(1A.1I) where one counts = i/iVoe*'^, = i/iVoe"*'^ in accord with neglect 
of the commutator. The leading orders are then 0{Nq) and 0{Nq). One 
returns to terms of 0{y/l%) and 0(1) in a standard perturbative procedure 
as a subsequent step. 

Thus to 0{Nq), the perturbation V becomes 



2 . 



+t;(0, k; 0, A;)a^afc + 0; 0, k)aka^] (A. 2) 



The third and fourth terms of eq. (]A.2I) are standard Hartree Fock single 



particle energies and may be absorbed into i^o 
Ho = ^nk{Ek-^j) 

k 

Ek = efc + b(0,A;;0,A;) + t;(A;,0;0,A;)]Aro/2 (A.3) 

The result is a quadratic hamiltonian. This is diagonalized in the follow- 
ing (Bogoljubov) transformation 

bl = alal^ + plak (A.4) 

where the phase of a is e'^"^ and of /? is e"'"*'^ and 

Ek\akf3k\ = {\ak\^ + \f3k\'')\Vk\ (A.5) 

with \Vk\ = \v{k, —k;0,0)\. Canonical commutation relations for bk gives 
\ak\' - \(3k\' = I. 

Thus the unperturbed ground state is not an eigenstate of the total hamil- 
tonian H but is unstable. In particular, had we chosen this state as the initial 
state it would evolve to the true ground state through emission of k, —k pairs. 
We leave to the reader the pleasure to confirm the Nambu Goldstone theorem 
for this case 

[H- < n\H\n >]=J2 '^kbibk (A.6) 
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where limfc^o<^fc = C\k\ + 0{k'^). He may also show that h\.\Q) corresponds 
to the creation of a longitudinal density fluctuation in this approximation, 
i.e. a phonon. 
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Appendix B 

Functional Integral Technique 



We use standard field theoretical techniques. The uninformed reader will find 
an account in refs. [HE] [10] [ZS] . The in vacuum to out vacuum amplitude is 



P0 e^^^"^) = (out,0, |0,m) (B.l) 



where (/) is a complex scalar field, 5'(</)) is the action to go from initial to final 
configurations [ which we take to be a quadratic form: free field theory in the 
presence of an electromagnetic and/or gravitational field ] in time t (which 
tends to oo ) and the mass^ has a small negative imaginary part. The mass 
dependence of S is of the form / —m'^(j)(j)*d'^x (with ^/g = 1) so that 



dW 



dm? 



(i x(oMt,0, 10 |0,m) = - j GF{x,x)d X 
tiGp , (B.2) 



where Gf{x, x') being the Feynman propagator to go from x' to x. In terms 
of the heat kernel K one has 



/■oo . 

Gp{x,x')= ds e-'"" 'K{x,x';s) (B.3) 
Jo 



where 



d 

(n + i7r)K = i6(s)6(x-x') (B.4) 

OS 
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In eq. (1B.2P we have written dW/dm^ as a trace. Clearly this is formal 
and one has to watch one's step on the measure. The following steps are valid 
because the passage from x representation to u representation is unitary [73] . 
Putting it all together we have upon integrating eq. flB.2l) over 

W = - —e-'^^^iKis) (B.5) 

JO s 

We now specialize to the case of constant electric field in the gauge = 0, 
At = Ex. As in Chapter [1], we can label the modes by u and work in the u 
representation, whereupon 

W=-tJ2 —e-'"" ' / K^{u, u- s)du (B.6) 

ij Jo S J 

where K^{u,Uq, s) obeys 

[ud^ + \- dr]K^ = 8{t)6{u - u^) (B.7) 

and r = 2E s, and we used eq. f ll.7p with ie replaced by d/dr. The solution 
is 

K^{u, uo; t) = e{T)6{ue^ - uo)e^^^ (B.8) 
to give from eq. (IB. 51) 

i f°° ds e~*™^* 



2 ^ jQ s sinhE's 
= i^ln(l + e-"™'/^) (B.9) 

where the last equality is obtained by picking up the poles (with = w?—ie) 
recovering therefore the Schwinger formula eq. fll.44p since = ELT /2%. 

When the path parameter r is expressed in terms of proper time these 
poles on the imaginary axis are related to multiple excursions through the 
tunneling region. For example, Born approximation (WKB) for the tunneling 
amplitude corresponds to the pole at Es = in and represents one excursion 
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back and forth {ATproper = 27r/a in the movement of the wave packets of 
Section 2.3). 

One may perform similar tricks to put into evidence the instabihty of 
Schwarzschild vacuum in the complete space spanned by the Eddington 
Finkelstein coordinates. Here the Schwinger counting parameter s is once 
more proportional to proper time. The complete analysis is considerately 
more tedious and complicated than the above but the essential features are 
the same. We refer the reader to ref. [73] for details. 

In this case it turns out that the convenient variable that describes the 
effective motion of a packet is p, the momentum conjugate to the Eddington 
Finkelstein coordinate x = r — 2M at fixed v. There are three classes of 
paths which contribute to W, those in which initial and final momenta have 
the same sign and that in which it goes from negative to positive momen- 
tum. Upon taking the trace initial and final momenta are set equal, so this 
operation requires a careful limiting procedure. One finds that it is the third 
class that encodes the instability and that the time to execute this movement 
is Av = iSnM = i/Sn i-e. jSn is the imaginary time to go from p to —p and 
back. One finds 

T r°° 

ImWBH = —^ diu \n\l- e"^«"] 
47r JO 

(B.IO) 

precisely the one dimensional partition function. In this manner (3h^^ is 
indeed interpretable as the action for a Hawking photon to tunnel out into 
existence. 
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Appendix C 



Pre- and Post-Selection, Weak 
Measurements. 

Pre- and post-selection consists in specifying both the initial and the final 
state of a system (denoted by S in the sequel). Pre and post selection is not 
an unusual procedure in physics. For instance when dealing with transition 
amplitudes, scattering amphtudes, etc... one is performing pre and post 
selection. 

In the first part of this appendix we shall implement post-selection in a 
rather formal way by acting on the state with projection operators which se- 
lect the desired final state(s) following the treatment of [63]. This generalizes 
the approach of 

In the second part of this appendix we show how post- selection may be re- 
alized operationally following the rules of quantum mechanics by coupling to 
S an additional system in a metastable state (the "post selector" PS) which 
will make a transition only if the system is in the required final state(s). The 
weak value of an operator obtained in this manner changes as time goes by 
from an asymmetric form to an expectation value, thereby making contact 
with more familiar physics. This extended formalism finds important ap- 
plication when considering the physics of the accelerated detector since the 
accelerated detector itself plays the role of post selector. In this way one can 
study the EPR correlations between the state of a uniformly accelerated de- 
tector and the radiation field, thereby clarifying and generalizing the results 

of [M], m, m- 
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C.l Weak Values 



The approach developed by Aharonov et al.^ for studying pre- and post- 
selected ensembles consists in performing at an intermediate time a "weak 
measurement" on S. In essence one studies the first order effect of S (ie. 
the back reaction) onto an additional system taken by Aharonov et al. to be 
the measuring device. But the formalism is more general. Indeed when the 
first order (or weak-coupling) approximation is valid, the backreaction takes 
a simple and universal form governed by a c- number, the "weak value" of 
the operator which controls the interaction. 

The system to be studied is in the state \ipi) at time tj (or more generally 
is described by a density matrix pi). The unperturbed time evolution of this 
pre-selected state can be described by the following density matrix 

Psit) = Usit,UM){tlj,\UsiU,t) (C.l) 

where Us = exp{—iHst) is the time evolution operator for the system S. 
The post-selection at time tf consists in specifying that the system belongs 
to a certain subspace, TCg, of Hs- Then the probability to find the system 
in this subspace at time tf is 

Pno = rr5[n°p(t/)] = Trs[u'sUs{tf,UM){^,\Us{U,tf)] (C.2) 

where 11^ is the projection operator onto TY^ and Trs is the trace over the 
states of system S. In the special cases wherein the specification of the final 
state is to be in a pure state \ipf) (ie. II^ = \ipf){ipf\) then the probability is 
simply given by the overlap 

Pf = m\Us{tf,UM)\' (C.3) 

Following Aharonov et al. we introduce an additional system, called the 
"weak detector" (WD), coupled to S. The interaction hamiltonian between 
5* and WD is taken to be of the form Hs-woit) = ef{t)AsBwD where e is a 
coupling constant, f{t) is a c- number function. As and -Bvfd are hermitian 
operators acting on S and WD respectively. 

Then to first order in e (the coupling is weak), the evolution of the coupled 
system S and WD is given by 

pitf) = \^itf)){^itj)\ 
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where 



= [Usitf,ti)UwDitf,ti)-ie [ ' dt Us{tf,t)UwD{tf,t)f{t)AsBwDX 
Us{t,U)UwD{t,ti)]\iPi)\WD) (C.4) 

where Us and Uwd are the free evolution operators for S and WD and \ WD) 
is the initial state of WD. Upon post-selecting at t = tf that S belongs to 
the subspace Ti^. and tracing over the states of the system S, the reduced 
density matrix describing the WD is obtained 



pwDitf) = Trslu'^sPitf) 



(C.5) 



In the first order approximation in which we are working it takes a very 
simple form 

pwnitf) ~ Puo\^WD{tf)){^WD{tf)\ 

where 

|*Wd(^/)) = UwD{tf,ti) -ie dtUwD{tf,t)f{t)Asyjeak{t)BwDUwD{t,ti) 

Jti 



where is the probability to be in subspace 71% and 



\WD) 
(C.6) 



^Sweak 



it)- 



Trs 


U''sUs(tf,t)AsUs{t,UM){A\Us{t 




Trs 


U^sUs{tf,tiM){ilJi\Us{U,tf) 





(C.7) 



is a c-number called the weak value of A. If one specifies completely the final 
state, = then the result of Aharonov et al. obtains: 



A 



Sweak 



it) 



{rf\Us(ff.f)AsUs(f.f,m 

{^jjf\Us{tf,uM) 



The principal feature of the above formalism is its independence on the in- 
ternal structure of the WD. The first order backreaction of S onto WD is 
universal: it is always controlled by the c-number Asweakit), the "weak value 
of A" . Therefore if S is coupled to itself by an interaction hamiltonian, the 
backreaction will be controlled by the weak value of i^int in first order per- 
turbation theory. For instance the modification of the probability that the 
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final state belongs to Tig is given by the imaginary part of i^int weak- Indeed 

(1 -ijdt //int)p.(l -^jdt H,nt) 

weak ) (C.9) 



The weak value of A is complex. By performing a series of measurements 
on WD and by varying the coupling function f{t), the real and imaginary 
part of Asweak could in principle be determined. Here the word "measure- 
ment" must be understood in its usual quantum sense: the average over 
repeated realizations of the same situation. This means that the weak value 
of As should also be understood as an average. The fluctuations around 
Asweak djce eucoded in the second order terms of eq. flC.4|) which have been 
neglected in eq. flC.4l) . 

To illustrate the role of the real and imaginary parts of Asweak, we recall 
the example of Aharonov et al consisting of a weak detector which has one de- 
gree of freedom q, with a gaussian initial state < q\WD > = e~'^ ,— cxd < 
q < +00. The unperturbed hamiltonian of WD is taken to vanish (hence 

UwD{ii^t2) = 1) and the interaction hamiltonian is ifs'_vKD(i^) = ^Sit—to)pAs 
where p is the momentum conjugate to q. Then after the post-selection the 
state of the WD is given to first order by 



<q\WDitf)> = {l-tepAsweak{to))e-'^"^^^ 

-((7-eAs„eafe(*0))V2A2 



e 

-(g-eReAs„eafc(to))^/2A2 +jegImyls„e^fc(to)/A2 



(C.IO) 



The real part of Asweak induces a translation of the center of the gaussian, 
the imaginary part a change in the momentum. Their effect on the WD is 
therefore measurable. The validity of the first order approximation requires 

eAsyjeak/^ << 1- 

It is instructive to see how unitarity is realised in the above formalism. 
Take II^ to be a complete orthogonal set of projectors acting on the Hilbert 
space of S. Denote by Pj the probability that the final state of the system 
belong to the subspace spanned by 11^ and by A-'^^^^i^ the corresponding weak 
value of A. Then the mean value of As is 

mAsm=T.PAeak (C.ll) 
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Thus the mean backreaction if no post-selection is performed is the average 
over the post-selected backreactions (in the linear response approximation). 
Notice that the imaginary parts of the weak values necessarily cancel since 
the l.h.s. of eq. flC.ll|) is real. Equation flC.ll|) is the short cut used in the 
main text to obtain with minimum effort the weak values. 



C.2 Physical Implementation of Post Selec- 
tion 

Up to now the postselection has been implemented by projecting by hand 
the state of the system onto a certain subspace Ti^.. Such a projection may 
be realised operationally by introducing an additional quantum system, a 
"post-selector" (PS), coupled in such a way that it will make a transition if 
and only if the system S is in the required final state. Then by considering 
only that subspace of the final states in which PS has made the transition, 
a post-selected state is specified. This quantum description of the post- 
selection is similar in spirit to the measurement theory developed in ref . |95] : 
by introducing explicitly the measuring device in the hamiltonian the collapse 
of the wave function ceases to be a necessary concomitant of measurement 
theory. As we have mentioned, this formalism is the basis for a general 
treatment of the energy density correlated to transitions of an accelerated 
detector. 

We shall consider the very simple model of a PS having two states, ini- 
tially in the ground state, and coupled to the system by an interaction of the 
form 

Hs-ps = Xgma^Qs + aQl) (C.12) 

where A is a coupling constant, g(t) a time dependent function, the oper- 
ator that induce transitions from the ground state to the exited state of the 
PS, Qs an operator acting on the system S. The postselection is performed 
dX t = tf and consists in finding the PS in the exited state. 

For simplicity we shall work to second order in A (although in principle 
the interaction of PS with 5* need not be weak). The wave function of the 
combined system S + WD + PS is in interaction representation to order e 
and order A^ 
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1-iJdt {Hs-woit) + Hs-Ps{t)) 

-Ijdtj dt'T [Hs-Ps{t)Hs-Ps{t')] - j dt j dt'T [Hs-wD{t)Hs-Ps{t')] 

+ '2! dtj dt' j dt"T[Hs-Ps{t)Hs-Ps{t')Hs-wD{t")'^ |0p5) 

(C.13) 

where \Qps) is the ground state of PS and T is the time ordering operator. 
The probabihty of being in the excited state at t = t/- at order is 



excited 



/ dtg{t)Q\ / dt'g{t')Qs\iJ,) (C.14) 



Upon imposing that the PS be in its excited state ai t = tj the resulting 
wave function is, to order e and A^, 



-I / dt\g{t)Qs{t) 



-JdtJ dt'r[ef{t)As{t)BwD{t)X9it')Qsit')] 



|V^,)|iyZ})at|0p5) 

(C.15) 



Making a density matrix out of the state ( jC.151) . tracing over the states of 
S and PS yields the reduced density matrix |\E'TyD)(^vyD| of WD to order e 
where 



\^wd) 



l-iej dtof{to)BwD{to)AZe:k{to, 



\WD) (C.16) 



and 



^e.eued(^ . ^ I dtg{t)Qlit) I dt' g{t')r [As{t,)Qs{t')] 1^,) 
{A\Idtg{t)Ql{t)Jdt'g{t')Qs{t'm 

Note how the weak value of As results from the quantum mechanical inter- 
ference of the two terms in eq. (IC.lSp . 

There are several important cases when the time ordering in eq. (1C.17P 
simplifies. If g{t) is non vanishing only after t = to then Aw takes a typical 
(for a weak value) asymmetric form 

j^e.cued(^ ^ ^ ^ I dtg^jit) J dt' g{t')Qs{t')As{toM) 
' midtg{t)QUt)Idt'g{t')Qs{t'm 
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If in addition g(t) = S{t — tj), tj > and Qs = H^, eq. (1C.7P is recovered 
using (J1%Y = n°. This is expected since in this case the post-selector has 
simply gotten correlated to the system in the subspace TCg 

If on the other hand g{t) is non vanishing only before t = to then the time 
ordering operator becomes trivial once more and eq. flC.17|) takes the form 

^excUed^^ ^ ^ {H I dtgmljt) Asjto) J dt' git')Q sji'M) 

midtg{t)Ql{t)Jdt'g{t')Qs{t'M) 

This is by construction the expectation value of As if the PS has made 
a transition. It is necessarily real contrary to eq. flC.181) when the weak 
measurement is performed before the 'collapse' induced by the post-selection. 

Finally, the weak value of As if the PS has not made a transition can also 
be computed. Once more the two cases discussed in eqs (jC.181) and (ICIOP 
are particularly simple: if g{t) is non vanishing only after t = to one finds 

Af:::r{to) = -— ^ — (imsm 

-^excited \ 

-A2Re[(^,| J dtg{t)Ql{t) J dt'g{t')Qsit')As{toM) 

(C.20) 

On the other hand if g{t) is non vanishing only before t = to one finds 

J- -l^excited \ 

-iA^Re 



mAsito) I dt I dt'Tgit)QUt)git')Qsit'M^)\j 

(C.21) 



These are related to the mean value of As and to eq. (]C.17|) through the 
unitary relation eq. (1C.11I) : if g{t) is non vanishing only after t = to 

PexcUedAZfal'ito) + (1 - PexcUed) AfZ^t^' = (^.l^slV'.) (C.22) 

if g{t) is non vanishing only before t = to 

p A excited/ J \ i /i JD \ Adeexcited 

J^excited^Sweak V^O) + U " J^excited) ^Sweak " 
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{^,\As\^^) + X^^^l J dtg{t)Ql{t)As{to) j dt' g{t')Q\{t')\^,) - 
-\^Re{ij,\As{t,) JdtJ dt'Tg{t)Q\{t)g{t')Q\{t')\^,) (C.23) 

where the right hand side is the average value of As before eq. (]C.22p and 
after eq. (1C.23I) the detector has interacted with S. 
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Appendix D 



S-wave Hawking Radiation for 
a General Collapsing Spherical 
Star Without Back Reaction 



We shall consider a collapsing sphere of matter with a well defined surface. 
Exterior to the surface the geometry is Schwarzschild, parametrized by a 
fixed mass M. This part of the space shall be coordinatized by advanced 
Eddington-Finkelstein coordinates {v,r,6,(f) defined in eq. (13. 8p . Interior to 
the star we shall use a set (T, X, 6, (p) in terms of which the length interval 
for a general dynamic spherically symmetric space is 



with dfl^ = d^^ + sin^ Odip"^. This system is used in the numerical calculation 
of (T^i/) that was reported in Section [3731 (with U,V = T ^ X). In that case 
it was used to cover the whole space so as to account for the back reaction 
outside the star as well. Since we are not taking into account the back 
reaction in the present case, we use eq. llD.ll) in the inside region only. The 
curvature stemming from the metric components a, b is driven by classical 
sources. 

The coordinate X is a radial coordinate, but areas of spheres are Airb"^. 
Thus at a given point of coordinates (T, X) on the star's surface one has the 
identification 



ds^ = -a\T,X) [dT^ - dX^] + 62(T,X) ciO' 



(D.l) 



r = 



6(r,x) 



(D.2) 
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Wc shall choose the origin of X so that r = coincides with X = 0. 

Hawking radiation follows from the combined gravitational and Doppler 
shifts that a photon experiences on its voyage from X^to X+. It begins on 

as a packet (superposition) of modes e~'^'^'" and end up on X+ as a packet 
of modes e*'^-''^"^ . We consider a particular radial packet that emanates from 
a point on X~ denoted hy v = u. The function f{u) varies according to the 
value of u so we will call it u{u) . The shift of frequency on X+ is 

uj{u)\.j-+ = uj{du/du) . (D.3) 

One may deduce this, for example, by the conservation of the number of 
oscillations in a given segment of wavefront (= urdv = urdu = Lj\j+du) 
There are three important points which the ray visits. These are: 

1. Pj = {X,T) where the ray penetrates the star. This labels the sphere 
which is the intersection between the incoming lightcone v = u and the 
star's surface. 

2. Pc where reflection on the axis X — occurs. 

3. Pq = (X,T), where the ray leaves the star. 

One has 

f - X = T + X (DA) 

since the coordinate time interval (= T — T) is equal to the coordinate 
distance toward in the star {= X + X). 

We now must specify the surface of the trajectory. In any given model this 

specification is, of course, correlated to the metric coefficients a, h along the 
trajectory, but there is no need to keep track of this in the general analysis. 
We parametrize the trajectory of the surface in each of the systems according 
to 

V = Vs(r) , 

X = Es{T) . (D.5) 

Then the first intersection point Pj is given in terms of the ray u by the 
solution of 
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u = Vs{r) 

X = Ss(T) (D.6) 
The second intersection point is obtained from eq. flD.4p 

f -Ss(r)=T + Hs(f) (D.7) 

The wave then propagates out to X+ according to e*'^^ where U = T — X. 
The exterior coordinate u is given by 

u = Vs{f) - 2r\f) (D.8) 

where f = 6(T, Ss(T)). Inverting this last relation gives T{f). One then 
computes X{f) from eq. (lD.2p hence lJ{f). Together with eq. OD.Sp this chain 
then gives the required relation between u and U . This matter is the sought 
function U = u{u). In this way it is seen that the phase of the outgoing wave 
gives a sort of " X-ray picture " of the interior of the star. 

To perform these various inversions is an arduous task. At the end of the 
section we shall present the results of a calculation wherein the collapsing 
star consists of dust. As announced in Section 13.41 almost all of Hawking 
radiation occurs in a small interval where the point Pq = {X,T) is near the 
horizon, Ti. As in Section 13.21 we shall linearize the equation of motion of 
the star's trajectory for such points. Similarly, the point P/ = (X, T) is near 
the extension of 7i into a past lightcone (the last null rays that are reflected 
into a future lightcone expanding up to the asymptotically flat infinity). See 
Fig. (13. 2b ). We shall denote the intersection of the star's surface with Ti 
as O = (Xo,To) and the intersection with the backward extension of H as 
/ = {Xj, Tj). The linearized forms of the trajectory of the surface of the star 
for points near O and / are 

v-vo = ko{r-2M) , 

V — vj = ki{r — rj) , (D.9) 
which in interior coordinates will be written 

X-Xo = MT-To) , 

X~Xj = Pj{T-Ti) . (D.IO) 
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The relation between jSo and ko (or /?/ and fc/) is obtained by equating 
the expressions of proper time intervals (= ds^) on the star's surface at the 
intersection points O and /, when expressed in both coordinate systems. 
Thus 

(1 _ —) dv^ - 2drdv = o'ldT^ - dX^] (D.ll) 
r 

with all differentials taken along the trajectory of the star's surface. In such 
differentials we also have 

dr = b'dX + bdT . (D.12) 

Dots are derivatives with respect to T and primes derivatives with respect 
to X. Using eqs.dES [PlOl [PlTl [Dl2|) gives at X = Xq 



[{l-^)ej-2kr][bj + bM' = aj{l-(3]) (D.13) 

-2ko[bo + b'oPo? = ao{l~l3'o) (D.14) 

Here qq, bo, {aj, bj) denotes the values of the metric components a and b at 
points O and / respectively. 

Let us now track the various chain of variables by following the ray back- 
wards in time. Near O we have on the suface point {X,T) the value of u 
given by 

f - 2M 

u ^ vq - AM - AMln (D.15) 

2M ^ ' 

where 

f = Hs(r) = 2M + boif - To) + b'o'^^if - To) 

so that 

f-2M={bo + b'o Po){T - To) (D.16) 
And from eq. (lD.7p we have 

f-To = T-Tj+ (Ss(T) - Ss(T,)) + (Ss(f ) - Ss(To)) (D.17) 

to give after linearization 

{T-To) = ^^(T-Ti) (D.18) 
J- — Po 
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Equations (ID. 151 ID. 161 ID.lSp then permit one to express u in terms of (T — 
Tj). We must go one step further back and express (T — Tj) in terms of 
V — vj. For this we use 



Vi 



ki 



ri 



b'j{X-Xj)+bi{T-Ti) = {b'jf3i + bj){T-Ti) . (D.19) 



Substituting for (T — T/) the value {v — vj) of eq. ( ID. 191) and proceeding the 
chainwise through eqs. (ID. 181 ID.16t ID.15|) and using the relationships eqs. 
( Im3|) and (iDlil) yields the final result 



u 



vo - AM -AM In 



'-^)k]-2kj ao 



kiy/—2ko 



ai 



1 


f/3o 




-l3o 


1 


-Pi 




f/3/ 



Vl 



2M 



(D.20) 



for the coordinate u of the outgoing ray in the exterior space which began 
on X_ at the point v = u. 

Equation (ID. 201) has a physical interpretation when one writes the total shift 
as a product of three shifts (denoted Di, D2 and D^).'T\ie factor Di is the 
shift produced in the voyage from 1~ to /, D2 from / to O and from O 
to X+. The identifications are 



D3 



k, 
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It is of course D3 that is the important shift that gives rise to the steady 
state Hawking radiation. 

To illustrate this discussion and to exhibit how the Hawking flux reaches 
its asymptotic value we consider a star consisting a cloud of dust in parabolic 
collapse. The trajectory of the surface of the star is 
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so that V = at r = 2M. The interior metric is a Roberston- Walker one, 
given by 



ds' = (^) ' + dX^] + X' ' dn^ (D.22) 

where M is the mass of the star, and X = the equation of motion of its 
surface. (The value of fixes the density in the star). 

The figures flPTTl) and flDl2l) show the outgoing flux ( times 2Atc {2M^)) 
measured at infinity ($00) and by a free falling observer on the surface of 
the star The fluxes are calculated using the two dimensional energy 

momentum tensor discussed in Sect ion Their analytic expressions (where 
we have choosen X^, = lOM) are: 

^ 1 18p^/^ + 27p - 21pV^ - 15 

247r4M2p''(8 + 2V/2 + 26p+12p3/2 + 2p2 ^ ' 

^ = ^-^(pP + p-2Vp + 21n(l + yp) 

-2[p + ln(p-l)] (D.24) 
^ 1 -(4p3/^ + 24p + 39pV^ + 15) 

247r 4M2p^(8 + 24pi/2 + 26p + 12p3/2 + 2p2 ^ ' ^ 

^(l-p3/2) (D.26) 



2M 3 

where the parameter p = r/2M represents the radial Eddington-Flinkenstein 
coordinate of the surface of the star and the proper time measured along 
it, too being the minkowskian time at infinity. The Hawking flux attains 
its asymptotic value at infinity when the last term in eq. flD.24p becomes 
dominant, ie. when the approximation eq. flD.20p is valid. 
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Fig. ID. 21 Outgoing Hux seen by an observer on the surface of the star as a 
function of his proper time. 
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In connexion with the remark at the end of Section 1.3 concerning the 
thermal distribution of pairs produced in an external field and in connex- 
ion with the analysis of the mean energy emitted by a uniformly acceler- 
ated detector we would like to mention the work of Nikishov[Ni] (see also 
Myrhvold[My]) who anticipated the Unruh effect by an analysis of the pho- 
tons emitted by accelerated electrons. 

As mentioned in footnote 10 (pages 153 and 154), upon taking into ac- 
count recoil effects the properties of the emitted ffuxes are drastically mod- 
ified. One can show that the correlations which encode the weak values 
are also modified. The result is that due to recoil effects the weak value 
Eq. (2.149) vanishes on X~ . Upon taking into account quantum gravita- 
tional effects in the black hole problem (see the end of Section 3.5), such 
recoil effects may play an important role and could for instance modify the 
properties of the weak values. Indeed we emphasize that the peculiar prop- 
erties of l(rvv)ipi (Eq. (3.106)) result from the (unjustified) assumptions of 
a free field theory evolving in a given classical background geometry. We 
hope that the study of the weak value (Tyy)^. can be used to investigate the 
validity of both assumptions. 

In order to understand the role of the transplanckian frequencies that are 
involved in the emergence of Hawking quanta, Unruh showed in a recent pa- 
per [93] (see footnote 13, page 201), through a numerical analysis that Hawk- 
ing radiation is unaffected by a truncation of the free field spectrum at the 
Planck scale. We have investigated and extended his result in Ref. [BMPS95] 
where we show analytically how the appeal to transplanckian frequencies can 
be avoided whilst retaining the thermal spectrum of emitted particles. 
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